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Abstract This review paper outlines background information and covers recent ad- 
vances made via the analysis of spectra and images of prominence plasma and the 
increased sophistication of non-LTE {i.e., when there is a departure from Local Ther- 
modynamic Equilibrium) radiative transfer models. We first describe the spectral in- 
version techniques that have been used to infer the plasma parameters important for 
the general properties of the prominence plasma in both its cool core and the hotter 
prominence-corona transition region. We also review studies devoted to the observa- 
tion of bulk motions of the prominence plasma and to the determination of prominence 
mass. However, a simple inversion of spectroscopic data usually fails when the lines be- 
come optically thick at certain wavelengths. Therefore, complex non-LTE models be- 
come necessary. We thus present the basics of non-LTE radiative transfer theory and 
the associated multi-level radiative transfer problems. The main results of one- and 
two-dimensional models of the prominences and their fine-structures are presented. 
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We then discuss the energy balance in various prominence models. Finally, we outline 
the outstanding observational and theoretical questions, and the directions for future 
progress in our understanding of solar prominences. 

Keywords Solar Prominences • Spectroscopy • Radiative Transfer ■ Diagnostics • 
Modelhng 
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1 Introduction 

Solar prominences are large magnetic structures confining a cool (temperature T < 
10* K) and dense (electron density 10^ — 10^^ cm~^) plasma in the hot solar corona. 
Typically, the prominence plasma is 100 times cooler and denser than its coronal sur- 
roundings. This situation raises important questions about the origin of the prominence 
plasma, and the energy and force equilibria which allow it to remain in the corona for 
a relatively long time. 

Many important characteristics of solar prominences, along with a presentation of 
some of the mathematical and physical tools necessary to un derstand the observation s, 
are given in the monograph by Einar Tandberg-Hanssen ( Tandberg-HanssenI Il995l l . 



Because of the impressive amount of exquisite observations obtained since that book 
was published some fifteen years ago, and the ever increasing power of numerical models 
and simulations developed in order to explain these observations, we felt that it was 
time to present an updated vision of the current status of our knowledge in the field 
of solar prominence physics. 

Early ground-based observations {e.g., in the Ha line of hydrogen at 6564 A) have 
led to a distinction being made between a prominence, which appears in emission when 
it is seen above the solar limb, and a filament when it is seen against the solar disk 
in absorption. This fundamental difference in the appearance of a prominence is easily 
explained in terms of opacity and emissivity of the plasma. In this review, we will 
often use both terms interchangeably as they refer to the same object. The prominence 
radiation in Hq allows us to view the densest parts of the structure in detail. The 
Hof line is routinely observed from the ground as well as from s pace (Fig. [J), as it is 
currently done by the S OT instrument (So lar Optical Telescope, iTsuneta et ahllioOSl ) 
on the Hinode satellite ( Kosuei et al.l[2007l l. 



A further distinction is usually made between quiescent and active prominences. 
For a discussion on the classification between different types of prominences and fil- 
aments, we refer the reader to Sects. 2.2.1 and 5 of the companion review paper 



Fig. 1 A quiescent solar pro minence observed above the limb in Ha with the Hinode/SOT 
NFI filter jHeinzel et al.ll200Sf) . 



( Mackav et alJEoiol . hereafter Paper II). Quiescent prominences are relatively stable 
features lying mostly outside active regions. Their lifetimes range from a few days up 
to several months (several solar rotations). Their dimensions are in the range of a 
few 10** — 10'°^ km in length, and a few 10'^ — 10'* km thick (with fine structures of 
the order of 10^ - 10^ km). Their heights are of the order of 10"* - 10^ km. Figure [2] 
shows several prominences observed with the Extreme ultraviolet Imaging Telescope 
(EIT, Pclaboudinicr c et al. 1995) on board the Solar Heliospheric Observatory (SOHO, 
IPomingQ et al.il 19951 ) on April 21, 2005 in the 304 A channel (where emission is pre- 
dominantly from He II). The size of the Earth is indicated for comparison. While they 
are globally stable, high-resolution movies from Ha ground-based and space-based ob- 
servations do show that even quiescent prominen ces exhibit rapid variations in the ir 
fine structure on time-scales of a few minutes (e.g. . iLiii et al.l[2007l : lBerger et al.ll2008l V 

Active prominences are dynamical features typically occurring in the vicinity of 
active regions and are usually short-lived (their lifetime is smaller than the lifetime 
of the associated active region). Their altitude is often smaller than that of quiescent 



Fig. 2 Quiescent and eruptive prominences observed by SOHO/EIT on 21 April 2005 in the 
304 A channel where emission comes predominantly from He II. 

prominences. Both quiescent and active prominences form an integral part of the solar 
corona. The conditions for their stability are inherently related to the associated mag- 
netic configuration (see Paper II, Sect. 2). Prominences and filaments may undergo 
large-scale instabilities which will disrupt their equilibria and lead to eruptions. These 
eruptions are often associated with flares and Coronal Mass Ejections (CMEs). 

In this review, we are primarily interested in the determination of the prominence 
plasma parameters which play a role in the force and energy equilibria, while Pa- 
per II describes recent progress in four areas of prominence research: their magnetic 
structure, the dynamics of prominence plasmas (formation and flows), the dynamics 
of magneto-hydrodynamic (MHD) waves in prominences, and the formation and large- 
scale patterns of the filament channels in which prominences are located. Both papers 
focus primarily on non-eruptive prominences. 

We outline general goals of the prominence spectral diagnostics and non-LTE mod- 
elling. Non-LTE stands for departures from Local Thermodynamic Equilibrium, or 
LTE. Particular approaches and techniques, together with the most important results, 
are detailed in subsequent sections. The ultimate goal of the prominence spectral diag- 
nostics is the determination of all important physical (thermodynamic) quantities and 
their variations in space and time. Some of these quantities (as known in the eighties) 
are listed in Table [T] while others will be mentioned in following sections. Primarily 
we are interested in those quantities which play a role in establishing the prominence 
equilibria, namely the momentum and energy equilibrium. The other main physical 
quantity, the magnetic field, is presented in Paper II. 
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Table 1 Main prominence plasma parameters in the cool core and in the PCTR (prominence- 
to-corona transition region); Electron temperature, microturbulent v elocity, electron densi ty, 
gas pressure, hydrogen ionization ratio, flow velocities. Adapted from lEngvold et al.l lll99Cri . 



Physical Parameter Core PCTR 



r(K) 


4300 - 10000 


10* - 10^ 


it (km s-i) 


3-20 


30 


He (cm~3) 


109 - 10" 


10^^ - 10** 


p (dyn cm"'^) 


~ 0.02 - 1 


~ 0.2 


N(R+)/N(R°) 
V (km s-l) 


0.2 - 0.9 




~ 5 


~ 10 



Two basic approaches are used to obtain the information on the plasma parame- 
ters. The first one is a direct inversion of spectral data, which is possible under some 
specific conditions. Namely, in the case of both cool as well as hot, optically thin plas- 
mas (plasmas from which all emitted photons leave freely without interaction) , one can 
obtain spatially averaged values of some parameters, as we discuss later. For hotter op- 
tically thin plasmas, no radiative transfer is needed, but in some cases one has to solve 
the statistical equilibrium equations for multilevel atoms and ions in order to com- 
pute the line emissivities. Sect. [5] presents spectral inversion techniques for prominence 
plasma parameters in the cool part. The inferred values of the electron temperature, 
non-thermal velocities, electron density, ionization degree, pressure, and abundances, 
are given. Sect. [3] follows the same approach and focuses on the plasma of the PCTR 
(prominence-to-corona transition region). The mechanisms of line emission from an 
optically thin plasma are first described, and then the inferred values for the electron 
temperature, electron density, gas pressure, and non-thermal velocities are presented. 
The small-scale structure of the PCTR is also discussed. Sects. [3]and[S]review measure- 
ments of velocities and observations of various types of bulk motions in prominences, 
and the determination of the mass of the prominence plasma. 

The second approach requires the solution of the radiative transfer problem, using 
the coupled equations of radiative transfer, statistical equilibrium, and other constraint 
equations. This is generally called non-LTE radiative transfer, or non-LTE modelling. 
It is applied to multilevel atoms or ions, and is necessary in the case of optically thick 
plasmas. We devote a substantial part of this paper to prominence non-LTE modelling. 
Contrary to direct inversion techniques, the non-LTE modelling represents a forward 
method. This means that starting from a given prominence model (spatial distribution 
of temperature, pressure, gas density), one evaluates the excitation and ionization 
balance for given species, determines the opacities and emissivities, and finally solves 
the transfer equation along the line-of-sight to get the emergent synthetic spectrum. 
The latter is then compared to the observed spectrum. In this way, one can adjust 
the initial model in order to get an optimum agreement with the observations. This 
procedure is iterative and the final models are called semi-empirical models, provided 
that they are at least partially data-driven. In Sects. [5]and[71 we present the basics of the 
radiative transfer theory and of multilevel non-LTE problems relevant to prominences. 
Then, Sects. [8] and [9] present the basic results from one-dimensional (ID) and two- 
dimensional (2D) non-LTE radiative transfer computations, with particular emphasis 
on the hydrogen and helium spectra of model prominences. 

The radiation which prominences or filaments emit, and which we observe in the 
form of monochromatic images or in the form of spectra, has actually a two-fold impor- 
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tance for the prominence physics. On one hand, it provides us with the diagnostics of 
the prominence structure as we explained above. On the other hand, the radiation field 
inside the prominence plays a crucial role in the global energy budget of prominence 
structures via the radiative losses or gains, and has also to be considered when evalu- 
ating self-consistently the prominence internal structure (momentum equilibria). The 
prominence densities are generally low enough so that the radiation field determines 
to a large extent the plasma ionization state which enters the equation of state of the 
plasma - the relation between the gas pressure and density depends on the tempera- 
ture and on the degree of ionization of hydrogen and helium. Therefore, the non-LTE 
physics must be integrated into the global radiation magneto- hydrodynamical (RMHD) 
modelling in order to obtain consistent models. Such theoretical models which are not 
data-driven are usually called ab initio models. They should play a crucial role in our 
understanding of the prominence physics. Considerations on energy balance issues are 
given in Sect. llOl 

Finally, Sect. [11] summarises open issues related to the material reviewed here, and 
we make some suggestions on possible future directions to progress in our understanding 
of solar prominences, both observationally and theoretically. 

In the following, we describe various aspects of the prominence spectral diagnostics 
and non-LTE modelling in detail, and summarize the most important results. Our dis- 
cussion is focused mainly on optical, ultraviolet (UV), and extreme ultraviolet (EUV) 
diagnostics, with some extension to selected infrared (IR) lines. However, we don't 
include the complex issue of the prominence polarimetry. The Stokes polarimetry is 
currently used to diagnose the prominence/filament magnetic fields via the Hanle effect, 
and the results are briefly discussed in Paper II (Sect. 2.1.4). 



2 Spectral Inversion for the Cool Prominence Plasma 

In this section, the focus is on relatively direct measurements of temperature, density, 
pressure, and other quantities (such as ionization degree and abundances) from the 
direct inversion of Gaussian spectral lines. Methods involving extensive modelling for 
more complex c ases are described in Sect. [8l P revious reviews can be found in e.g., 
iHiravamal (|l990l ) and iTandberg-HanssenI (|l995h . 



2.1 Temperature and Non- Thermal Velocities 

The simplest way of deriving the temperature is to analyse the profile of the line. In 
this section, we limit ourselves to optically thin plasmas emitting spectral lines which 
are Gaussian-shaped. For a Gaussian-shaped line, the Doppler width is given by: 

AX,^'-l('-^Y\ (1) 
c \ m J 

where Aq is the rest wavelength position of the line, T is the ion temperature, and m is 
the mass of the ion. In the literature, we may find different values for T. In some cases 
the temperature of maximum ionization fraction is taken. In other cases, T corresponds 
to the temperature of formation of the line. 
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Gaussian-shaped spectral lines from optically thin plasmas have a width that gen- 
erally exceeds the thermal broadening. Unresolved motions can be responsible for ad- 
ditional broadening. If a line is broader than expected from the thermal component, 
an additional non-thermal term ^ is included in ([T}; 

Aq f2kT nV''^ 

a\d = — — + r . (2 

c \ m J 

The parameter ^ is the non-thermal velocity (NTV) component and may have different 
interpretations. Generally, it describes an averaged value along the line-of-sight (LOS) 
of the unresolved motions, due to waves or turbulence. Equation ^ also assumes 
that the non-thermal motions have a Maxwellian distribution. Uncertainties associated 
with the line width measurement and the determination of T affect the result for ^. 
This is particularly true for cool lin es produced by neutrals and singly ionized ions 
( Chae et al.lll99AlBaudin et al.ll2007h . for which it is difficult to establish the formation 
temperature. Estimating this quantity as a function of temperature is important for 
understanding the distribution of energy in the different layers of the prominence- 
corona interface. Important implications in the energy budget may be deduced (see 
Sect. HI). 

The thermal velocity and the NTV may be deduced by a line fitting technique 
and by using the measured Full Width at Half Maximum (FWHMmcas). Under the 
assumption that the instrumental profile is Gaussian (FWHMiiist)i the true solar line 
width (FWHM sun) is given by: 

FWHMsun^ = FWHMmcas^ " FWHMi„st^ • (3) 

An immediate technique for separating the temperature from ^ consists in measuring as 
m any lines from d ifferent atoms and io ns of different atomic ma sses as possible (Fig. 1 
of lHiravamalll990h . Results obtained bv lStellmacher et al.l (120031 ') with Ca II 8542 A and 
He I 10830 A are shown in Fig.O they derived a temperature in the range 8000 — 9000 K 
and ^ between 3 and 8 km s^^ in two quiescent prominences. It is difficult to assess 
whether the dispersion of observational values results from uncertainties or intrinsic 
variations inside each prominence. It should be pointed out that the technique assumes 
that the lines are optically thin. Otherwise, temperature and ^ are overestimated. 



2.2 Electron Density 

As is clear in the Hvar model (see Table [T] and lEngvold et al]|l990l ). and as we will see 
below, measured values of the electron density vary greatly. Some of these variations 
are, no doubt, due to differences between the various techniques that are used, but 
there are also variations among individual prominences, and between parts of the same 
prominence. 

A primary ground-based technique for measuring the electron density is to use the 
Stark effect. The Stark effect broadens lines in such a way that the de nsity can be esti- 
mated, e.g., from the number of the highest resolvable Balmer line fsee llnglis and Tellerl 
However, the broadening can also be due to electric field and instrumental broad- 
ening, so this technique essentially provides an up per limit on rte . The metho d was 
applied on numerous occasions by Hirayama fe.o. . lHiravamalll97ll . Il972l . Il985l ). who 
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0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 0.28 



AAq (Ca*8542 A) 

Fig. 3 Relation between observed Doppler widths of He I 10830 A and Ca II 8542 A for two 
prominences (dots and crosses). The calculated curves for kinetic temperatures of 8000 K and 
9000 K are shown along with the values of the NTVs (values in the squares in km s~^). From 
IStellmacher et al.l ll20oij'l . 



obtained upper limits on the electron density between lO"'^" and 10^"'^ cm~^ in quiescent 
prominences. 

Lines in the mid-infrared have also been used as a tool for electron density di- 
agnostics. Based on the observations of the H I line at 808.3 cm~^ (12.37 fim) in a 
quiescent prominence. | Zirker (^1985) obtained an upper limit to ric of 1.9 x 10^*^ cm^'^. 
IChang and Demind (| 19981 ') measured Stark broadening in the spectra of a quiescent 
and two time-frames of an active prominence. They found upper limits to Uc of 
2 X 10"'^'^ cm"'^ in the former, and 7.4 x 10^° and 3.5 x 10^" cm from the two spectra 
in the active region prominence. 

The Hanle effect, which is a depolarization due to the magnetic field, can also be 
used to calculate no by utilizing the fact that there is additional coll isional depolar- 
ization t hat in creases with density in some lines and not in others. iBommier et al.l 
1 1986al lbl. \l99^ ) used the Hanle effect to measure no in a set of combined 24 promi- 
nences, finding values between 1 x 10^ and 6.4 x 10^" cm^'^. 

Landman developed a technique for estimating no from ratios of resonance lines 
of Na, Mg, Sr"*", and Ba"*", which draws on some assumption on the intensity of the 
La line. His revised average v alue for a comb ination of quiescent and active region 
pro minences is 9 x 1 ^*^ cm ~^ ( Landmanll 19861 ). with a range of about a factor of two. 

iKoutchmv et al. 1 (| 19831 ) determined an average electron density limit of > 5.7 x 
10^ cm~'^ for a prominence observed during an eclipse, using a diagno stic based on the 
ratio o f the hydrogen H/3 intensity to that of the white-light continuum. Istellmacher et alJ 
(|l98d ) used Thompson scattering measurements during the same eclipse to estimate 
no in a highly activated/erupting prominence as > 2 x 10^ cm^"^. These calculations 
provide only lower limits because they assume a homogeneous source; a filling factor 
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less than one would l ead to a higher rie. A mor e complex diagnostics of this type was 
recently developed bv ljeicic and Heinzell (|2009l ) who derived similar density values. 

There have been a number of electron density estimates based on emission measure 
methods. The emission measure yields the amount of plasma emitting the observed 
radiation along the LOS, and is described in Sect. 13.2.11 Such methods require the 
modelling of the source vol u me, e specially along the LOS, and of the filling factor. 
For instance, ISastian et al determined the emission measure of a quiescent 

prominence from microwave data. By estimating the extent of the source along the LOS 
and assuming a filling factor between 0.03 and 0.3, they derive d an el ectron density in 
the range 1 x lO^'' — 3 x 10^^ cm~^. Similarly. iHarrison et al.| 1 1993h used microwave 
observations of an erupting pro minence t o calc ulate rie = lO^'' cm~^, assuming a filling 
factor of 0.33 and T = 6000 K. IlI et al.l (|l998l ) provided an emission measure estimate 
using Ha line fits combined with a cloud model. Using an estimate of the prominence 
thi ckness, they calculat ed no = 1.8 x 10^° cm~'^. 

iKucera et al.l ( 19931 ') observed what they thought to be second order plasma radi- 
ation emitted at 1.4 GHz from an erupting prominence. From this assumption they 
calculated a lower limit to the density of 6 x 10^ cm^'^. 

Thus we see that we have measurements of prominence electron density from 10^ 
to 10"'^^ cm~'^. Other (indirect) methods, such as non-LTE modelling (described in 
Sect. [8]), or prominence seismology (Paper H, Sect. 4.4), yield similar electron density 
values. The aim of on-going work is to shorten the two-orders-of-magnitude range of 
values, inasmuch as the large diversity of prominences allows for such an achievement. 



2.3 Ionization Degree 

The electron density is not the only density to be determined. If one wants to derive 
e.g., the total mass of the prominence or the gas pressure, one has to measure the 
density of neutral hydrogen (along with minor atoms and ions such as He). This raises 
the issue of the ionization degree of the hydrogen plasma, that we define as the ratio 
of the proton density (roughly equated here to the electron density no) to the neutral 
hydrogen density (roughly equated here to the population density of the ground level, 
ni). Actually, there are not many studies devoted to this issue. They require a simul- 
taneous derivation of the no and ni densities. We have discussed electron densities in 
the section above. Neutral hydrogen densities are addressed in Sect. (5] 



The earlier measurements from Skylab/A TM (Apollo Telescope Mount, Tousey 



1981 



:o pe 

1977 ) provided ric/ni values between 1 and 2 l Orrall and Schmahllll980l : iKanno et al 



From the ratio of La to Ca II K intensities m easured by the LPS P (Labora- 



toire de Physique Stellaire et Plan e taire) instrument ( Bonnet et al.lfl"978l ) on OSO-8 



(Orbiting Solar Observatorv). Ivial (|l982al ) deri ved a value of 3.4 (within a factor 3), 
definitely higher than the 0.09 value derived bv iLandmanl ( 19841 ) which led to suspi- 



ciously high pressures of the order of 6 dyn cm . The Hvar reference atmosphere gives 
iV(H+)/iV(H°) between 0.2 and 0.9 (Table[I]). 

The ionization degree of the plasma plays an important role in the momentum 
balance. Its importance with respect to the da mping of MH D waves in prominences is 
emphasized in Sect. 4.3.3 of Paper II (see also lOliveij[2009l ). 

At this stage, we should emphasize that only a comparison between the full set of 
profiles (including hydrogen) and the predictions of non-LTE modelling can provide an 
adequate answer for this important parameter. 
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2.4 Gas Pressure 



Early Rr ound-based effor ts to measure prominence gas pressure are discussed exten- 
sively bv lHiravamal(|l990l ). A number of early values measured with ATM, OSO-4, and 
OSO-6, yielded values between 0.01 and 0.08 dyn cm~^ for quiescent prominences, 
although iKieldseth-Moe et al.l 1 19791 ) found pressure values between 0.07 and 0.22 dyn 
cm~^ in a more active prominence. 

One often used pressure diagnostic i s the ratio of H/3 to Ca II 8542 A, originally 
suggested by iHeaslev and Milkev ll 19781). The method has yielded a range of values 
from low 0.01-0.04 dvn cm~^ (|Heaslev and Milkcv"l978':'Stellmacher and Wichr"l997'. 
[2OOQ) to high values of 0.3—0.5 dyn cm~^ ( Bcridli ri ct al. 1988). Gouttcbrozc and HcinzeJ 
(I2OO2I) revised models of the relev ant line emissions. They applied the results to the data 
of IStellmacher and Wiehil (|2000l ) and found the pressure to be above 0.1 dyn cm 
Although the values vary greatly from promine nce to prominence, the ratio tends to 
be fa irly uniform f or a given prominence ( e .g., de Boer et al.l 1 19981 : lOammasch et all 
I2OO3I ). For instance IStellmacher and Wiehii (|200ol ) report variations in the pressure of 
about 30% over a given pro minence. 

As discussed in Sect, [m iLandmanl ('1986') used ground-based spectra to calculate 
the electron density. He then derived a pressure of 1.4 dyn cm^^, a value which, like the 
density values (whether electron or neutral hydrogen), is higher than that calculated 
using the H/? to Ca II ratio discussed above. 



2.5 Abundances 

The precise measurement of abundances could be a precious source of information on 
the origin of prominence material, and hence on the prominence formation mechanisms. 
Establishing whether the prominence plasma more closely resembles the photospheric 
or chromospheric plasma tells us something about its formation. The different models 
for the formation of prominences invoke plasma condensation (where chromospheric 
abundances are expected) or siphon and emerging flux models, where prominence ma- 
terial originating from the photosphere is transported into the corona (leading to pho- 
tospheric abundances). We refer the reader to Sect. 3 of Paper II for an in-depth 
discussion of the different models of prominence formation. 

However, very few studies aimed at measuring quiescent prominence abundances 
exist in the literature, probably because of the difficulty in treating this partially-ionized 
en vironment . 

iMariskal l|l980l ) derived absolute element abundances for C, N, O, and Si in differ- 



ent solar regions, including a prominence, taking as reference a Si IV line from ATM 
measurements (which was then used to scale the abundances to hydrogen). His aver- 
aged values in the prominence indicated the N abundance to be less, C as abundant 
as, and O more abundant than in the quiet Sun. However, there are important uncer- 
tainties in these measurements, primarily arising from the atomic physics calculations 
(uncertainties of at least 20%). 

The Ne/Mg abundance ratio in quiescent prominences was obtained from Skylab 



by ISpicer et al. (|l998l ). Their results, obtained under the assumption that the two 
lines used are formed a t the same tem perature, supported a photospheric origin for 
the prominence plasma. ISheelevI (|l995l ) arrived to a similar conclusion by finding Ne- 



enriched plasma (that is, close to the photospheric value) in emerging flux regions and 
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990 1000 1010 1020 1030 

Wavelength (Angstroms) 

Fig. 4 A spectrum obtained by SOHO/SUMER in S VI, LjS, and O VI. A portion of the disk, 
a filament, and the filament extension above the limb as a promi nence are visible. A, B , and 
C point to three different parts of the prominence. Adapted from lSchmieder et al.l (119991) . 



close to the footpoints of a filament. I Widing et al. derived the same ratio for an 

erupting prominence and found a value intermediate between that for the photosphere 
and corona. This was probably due to the heating and mixing of the prominence plasma 
with the surrounding coronal environment. 

These analyses are difficult to perform, because of the uncertainties in the absolute 
UV line photometry coupled with the uncertainties in many atomic parameters. 

Let us note that no n-LTE mod elling techniques can also be used to study abun- 
dances. For instance, He aslev and M ilkcv ( 197^) inferred a helium abundance of 0.100± 
0.025 by number with respect to hydrogen. 



3 Spectral Inversion for the Prominence-Corona Transition Region 

In this section, we describe plasma diagnostic techniques that apply to most of the 
EUV (100 - 1000 A) and UV (1000 - 4000 A) plasma emission, i.e., coming from 
the chromosphere (temperature of 6 x 10"^ — 10* K) to the corona (T > 10*^ K), with 
particular emphasis on the PCTR (T > 10* K). 

Since 1995, the two most important instruments for UV-EUV spectroscopy of 
prominences have been SUMER (Solar Ultraviolet Measurement of Emitted Radia- 
tion) and CDS (Coronal Diagnostic Spectrometer) aboard SOHO. Each instrument 
has its own advantages, and both are therefore complementary and extremely useful 
to study the PCTR. SUMER wavebands include a large number of chromospheric and 
transition region (TR) lines. In addition, it has a high spectral resolution. Figure |4] 
shows an example of a spectrum obtained by SUMER in S VI, L/3, and O VI. The 
solar disk, the filament on the disk, and the prominence above the limb, are visible 
on the spectrum. CDS consists of two types of spectrometers: the Normal Incidence 
Spectrometer (NIS) and the Grazing Incidence Spectrometer (CIS). CDS wavebands 
are richer in density s ensiti ve lines at coronal temperatures than those of SUMER. 
IPatsourakos and Viaj l|2002ll reviewed the great contribution of these instruments to 
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Table 2 Main cha racte ristics of current s pace- based spectrometers . References: 1 - 
IWilhelm et"all iHQQSI) : 2 - iHarrison et al.l jlQQSlj : 3 - ICuIhane et al.l l l2007f) . 



Instrument / Obser 
vatory 


Wavelength range 
(A, 1st order) 


Spectral pixel 
size 


Spatial pixel size Ref. 
(arcsec/pixel) 


SOHO/SUMER 


660 - 


- 1610 


42 - 45 mA 


1 


1 








pixel" ^ 






SOHO/CDS 


150 - 


- 800 


70 mA pixel" 1 


1.68 


2 








(NISI), 117 mA 












pixel- 1 (NIS2), 












400 mA for 12 












pixels (GIS) 






Hinodc/EIS 


170 


- 210 and 


22 mA pixel -1 


1 


3 




250 - 


- 290 









solar prominence physics, along with other SOHO instruments. These two spectrome- 
ters were joined in 2007 by EIS (Extreme-ultraviolet Imaging Spectrometer) on Hin- 
ode. The possibility of obtaining combined observations with these three instruments 
provides a rare opportunity for detailed spectroscopy of the solar atmosphere. This 
has led to a numb er of joint observation campaigns dedicated to prominences {e.g., 
iHeinzel et al.|[20o3 ) as well as other targets. 

We summarise the main characteristics (spatial and spectral resolution, wavelength 
coverage) of current space-based spectrometers which significantly contribute to solar 
prominence science in Table [2l The CDS spatial and spectral resolutions were affected 
by the temporary loss of the SOHO satellite in 1998. Post-recovery values are different 
(in some cases the performance was degraded), and these can be found on the CDS 
websiteQ. Post-recovery CDS data show distorted line profiles, but standard routines 
have been made available in the CDS software package to recover the original line 

profiles. 

Let us also mention the Ultraviolet Coronagraph Spectrometer (UVCS. lKohl et"al] 

Il995l ). the other SOHO spectrometer. UVCS can be used to study erupted prominences 
high in the corona {R > l.S-R©). Even though the spatial and spectral resolutions are 
lower than those of the low corona instruments, UVCS data are unique in their kind. 
More details on the di agnostic t echni q ues presented in this section can be found 
Feldman et al.1 (Il977l ): iMariskal (|l992l '): iMason and Monsignori Fossil (Il994l ): and 



IPhillips et alj These techniques apply to optically thin plasmas with electron 



densities rio < 10^^ cm~^. In such environments, the radiation which is emitted by 
the source completely escapes the source itself without being absorbed. This is true 
for most of the lines emitted in the PCTR, although some exceptions apply. There 
are many atomic processes involved in the formation of the observed radiation. These 
include excitation and de-excitation of atomic levels, and ionization and recombination 
for each ion. However, not all these processes act at the same rate, and in most cases 
of interest we can draw a simplified picture. Electron temperature, electron density, 
plasma velocities and other plasma parameters can be derived by measuring the ob- 
served spectral lines and applying the appropriate inversion technique. Table [3] lists 
the main optically thin lines used in recent works to diagnose the PCTR. This table 
also lists the line wavelength, the formation temperature of the line and the SOHO 
spectrometers that detect the line. 



^ |http : //solar .bnsc .rl . ac.uk/softMare/ugulde/ugTilde . shtml] 



13 



Table 3 Main lines detected by SOHO instruments that can be used to diagnose the plasma 
conditions in the PCTR: emitting ion, wavelength, formation temperature, and spectrometers 
detecting each line. 



Ion 


A (A) 


logT(K) 


Instrument 


Si II 


1259 


4.1 


SUMER 


C II 


1037 


4.4 


SUMER, UVCS 


Si III 


1206 


4.7 


SUMER 


C III 


977, 1175 


4.8 


SUMER 


Si IV 


1402 


4.9 


SUMER 


O III 


525 


4.9 


SUMER, CDS 


N V 


1242 


5.3 


SUMER, UVCS 


O IV 


554 


5.3 


CDS 


O V 


1218 


5.4 


SUMER 


O VI 


1037 


5.6 


SUMER, UVCS 



The UV-EUV prominence emission also contains various continua emitted from 
the cool part of the structure (mainly hydrogen). Specific diagnostic techniques can 
be applied to derive, for example, the plasma temperature from the Lyman continuum 
(see Sect. 1511) . 

3.1 Line Emission from an Optically Thin Plasma 

The optically thin emission discussed in this section is assumed to come from a plasma 
in ionization equilibrium: the population of a given ionization state of an element is 
constant in time. All the lines listed in Table [3] are optically thin in prominences. 
Optically thick lines are discussed in Sects. |S] and 

The emissivity of the transition (power per unit volume per unit solid angle) due 
to the spontaneous decay of electrons from the ion's excited level j to lower level i 
{bound-bound transition) is given by; 

Pi^ji) = 4^;^.. ^3^j» [erg cm s sr ] , (4) 

where nj represents the number density of the level j, and Aji is the spontaneous 
radiative decay rate from level j to level i. 

Under the PCTR density conditions, the ionization and recombination processes 
are slower than the coUisional excitation in populating excited levels. In this case the 
ionization state of an element can be treated independently from the excitation balance. 
The number density rij is obtained by solving the statistical equilibrium equations for 
a number of low-lying levels and including all the relevant excitation and de-excitation 
mechanisms. 

Further simplification arises under the coronal approximation. This assumes that 
only the ground level (<?) and the excited level (j) are responsible for the emitted radia- 
tion, and it is generally valid for resonance lines. In this case the statistical equilibrium 
equation is reduced to: 

TLcTLgC gj = TljA.jg , (5) 

where Cgj is the coUisional excitation rate and no the electron number density. 

In this case the problem of calculating the emission is reduced to considering ex- 
citation from the ground level to the excited one, balanced by spontaneous radiative 
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decay back to the ground level. The line intensity is then given by the integration of 
(|4]) in the emitting volume V of cross-sectional area A, along the LOS: 

H>']g) = ^ j ncUgCgjdV [erg cm s sr ] , (6) 

where ^ was also used. The population of the g level can be expressed in the form of: 

ng iVioniVcliVH 



■ne , (7) 



where Ug/Ni^^^ ~ 1 under the assumed conditions, Ni^^^/N^i is the ion density relative 
to the total number density of the element; N^i/Nj^ is the abundance Ab of the element 
with respect to hydrogen, and N-^/ric = 0.8 in the case of a fully ionised hydrogen- 
helium plasma. 

The coefficient Cgj is given by: 



where Tgj{T) is the thermally averaged collision strength (which is related to the 
electron excitation cross-section), ujg is the statistical weight of the groun d level, and 



a Ma xwellian velocity distribution for the colliding electrons is assumed (jPere et al.l 
Il997l l. 

Inserting (O and ([8| into ((G]), we obtain: 

where the element abundance is assumed constant in the emitting volume. We notice 
that the line intensity is proportional to n^. In ([9]), we have introduced the function 
G(T), called the contribution function, which contains all the atomic physics parameters 
and is peaked in temperature (see Fig. [S] for an example): 

^^^^ ^ 8.63 X io-^r,,(r) j^iVH^„,,,^^^^_ ^^^^ 

ujg N^i no ^jg^T^ ^jg 

Equation Q was derived under simple conditions. The coronal approximation does 
not always apply, so the statistical equilibrium equations should be solved for all low- 
lying level populations, including all the relevant excitation and de-excitation processes 
in a multi-level atom. For example, the proton collision excitation and de-excitation are 
important for those transitions where the energy between the levels involved is much 
smaller than kT. 

In addition, the calculation of Cgj is quite complex because it requires the knowl- 
edge of the electron excitation cross-secti on. The accuracy of this calculation depends 
mainly on the approximations used (see iMason and Monsignori Fossil 1 19941 . for a re- 
view) . 

If the radiative decay transition probability is small enough, other mechanisms 
may start to compete with it in de-exciting the level. This is the case for forbidden 
or intersystem transitions, which de-excite the level by electron collisions (for electron 
densities large enough to allow coUisional de-excitations to play an important role). 
The levels which are sensitive to this are called metastables, and their populations may 
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Fig. 5 Contribution functions for the lines belonging to the O III- VI ions listed in Table [3] 
(CHIANTI V. 5.1) 



become as important as that of the ground level. The intensities of their lines depend 
on the electron density, making them good candidates for electron density diagnostics 
(see Sec. imj) . 

The inversion technique aimed at deriving the plasma parameters from the spectral 
line intensities requires the knowledge of all the terms in ((SJ. This includes the ioniza- 
tion degree of the element under consideration Nio^/Nd in (O. This term is calculated 
by solving a set of equations invoking the ionization balance among the ionization 
stages from to z+1, where z is the stage of ionization of interest. For each stage 
k < z, all the relevant ionization and recombination processes are taken into account: 

NkQk = ^^(fe+i)«(fc+i) , (11) 

where q and a are the ionization and recombination rate coefficients, respectively. For 
an optically thin plasma, the dominant processes are: coUisional ionization, excitation 
follow ed by auto-ionization, and radiative and dielectronic recombinations (see lMariskal 
Il99l for more details on these processes). 

When dealing with spectroscopic diagnostic techniques to derive physical param- 
eters of the emitting plasma, it is often necessary to compare the measured emission 
with that expected under assumed conditions. There exist various atomic databases 
and spectral codes that allow the calc ulation of the line intensity from Some are 
freely available, e.g., CHIANTfl (|Dere et all 1 19971 ) and APEC/APED ()Smith et~al] 

^ |http : //mm . chianti . rl ■ ac . uk/ 1 
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l200lh . In the physics of the solar UV-EUV emission the most commonly used is CHI- 
ANTI. This database, built by an international collaboration (USA/Italy/UK), con- 
tains atomic data (energy levels, wavelengths, radiative transition probabilities and 
excitation data) for a large number of ions. It provides software for diagnostics analy- 
sis using the Interactive Data Language (IDL). This code is kept up-to-date with results 
of new atomic physics calculations and laboratory measurements. CHIANTI also in- 
cludes several ionization equilibrium calculations and elemental abundances. This is an 
important aspect that allows comparisons of results under different assumptions. The 
diagnostics analysis, in fact, depends on these assumptions. Note that CHIANTI does 
not take into account the effect of velocity fields. The investigator should also be aware 
that uncertainties arising from atomic physics calculation are usually at least 10%. 



3.2 Electron Temperature 

Knowing the thermal structure of prominences is an essential element to understand 
their energy balance. It is still not understood how a prominence can stay stable for 
days with a core plasma at about 8000 K surrounded by the 1 MK corona. There exist 
several spectroscopic techniques to derive the electron temperature from the observed 
optically thin emission. Here we present those which are most often used in the analysis 
of prominence observations: the emission measure and differential emission measure 
techniques, which both involve the electron density at a given temperature. 

3.2.1 Emission Measure 

The Emission Measu re (EM) , first introduced bv llvanov-Kholodnvi and Nikol'skiil (|l96lh 



andlPottasch 



lll963l) , is a quantity which yields the amount of plasma emissivity along 



PottaschI (l963h used a simple approach by noticing that for most lines, the 



the LOS. 

contribution function G{T) is a narrow, peaked function of temperature. This author 
assumed that G{T) can be approximated by a constant function, equal to 0.7 G(rmax), 
in an interval AT around the temperature Tmax where the function reaches its max- 
imum, and zero elsewhere. The temperature interval is such that ZiT/Tmax ~ 1 — 2. 
Under this condition, the optically thin line intensity can be rewritten as: 

/(A,<,) = i-A6{EM){G(r)) , (12) 

where (EM) is the column emission measure of the plasma at a temperature T in the 
interval AT, averaged along the LOS h: 

(EM) = / nldh [cn-T^] . (13) 
Jh 



In (|12|l . {G{T)) is the contribution function as defined bv lPottaschI ( 1963l l. The emission 
measure so defined can be estimated directly from the observed intensity of a single 
line. It should be remembered, however, that this is a crude estimation of the real 
emission measure which takes into account the emission of the spectral line over the 
total temperature interval. 

Em ission measure analys es on prominences have been performed since the time o f 
Skvlab. lMariska et all (|l979l ') derived the EM under the assumption of lPottasclil (|l964l 'l 



and found a steeper decrease of EM with height for hot material {T > 10^) than for 
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cooler areas. Further indications led the authors to conclude that their observations 
better fit a picture with cool threads having a thin TR (see Sect. 13.5)1 . However, while 
the core material remains constant with height in the prominence, the TR of the thread 
gets thinner. 

3.2.2 Differential Ermssion Measure 

When several lines formed at different temperatures are available, a more sophisticated 
technique may be used to derive the distribution in temperature of the plasma along 
the LOS, that is the differential emission measure, or DEM: 

DEM(T) = n2^ [cm-^K-l]. (14) 
In this case the total line intensity can be expressed by: 

= —Ah / G(r)DEM(r)dr . (15) 

The DEM contains information on the multi-temperature nature and the fine structure 
of the prominence along the LOS, and the dominant physical processes at work. To 
derive the DEM, (jlSp needs to be inverted using the observed line intensities, through 
the calculation of G{T) and assuming a constant elemental abundance. The wider 
the temperature range of the lines, the more the inversion procedure is constrained. 
However, this inversion is not an easy task and requires some assumptions concerning 
the solution. In addition, the solutions are sensitive to t he atomic physics calculations 
and their uncertainties. There exist different codes f e.g. . iMonsignori Fossi and Landinil 
ll99ll : lLandi and Landinilll997l ) available to solve the DEM problem. Note that one can 
find slightly different DEM definitions in the literature. For example, the DEM may be 
defined for a unit volume (~ n^AV /dT). This shows that some precaution is needed 
when comparing results from various sources. 

In the case of prominences, the DEM provides insight into the temperature gradient 
in the PCTR, which is one of the key parameters for the prominence modelling. A 
strong temperature gradient in the PCTR was identified decades ago in Skylab data. 
Skylab results indicated a lowe r pressure and thinner TR for prominences than for the 
chromosphere-corona TRre.g.. lvialll99d ). The radiative losses in the PCTR seem to be 
too high to be compensated by either waves or conduction (see e.g. . lAnzer and Heinzell 
[2OO8). However, these conclusions often rely on models in which, due to the lack of 
detailed information, important assumptions concerning the prominence fine structure 
and plasma conditions are made. 

Deriving the DEM from observations is not a trivial task. Furthermore, given a 
DEM, there is not a unique way to obtain the density and temperature distributions 
as functions of the path length h. It should be noted that the emission from promi- 
nences spans the whole UV-EUV temperature range, from 10'' to 10^ K. The DEM 
calculation does not provide any information on the filamentary origin of the promi- 
nence emission, so that a comparison with expected DEMs from different geometric 
models is needed. Based on observational evidence, there are two main models that 
could satisfy the observed thermal properties and relate them to the unresolved fine 
structure (see Fig. I10|) . The first model considers thre ads having a cold core and a 
thin TR interfacing with the coronal environment fe.o. . [chiuderi and Chiuderi-Dra"gol 
I1991I I. The observed EUV emission originates mostly from the thread's TR, and the 
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Fig. 6 DEMs f or a quiescent prominence from ISchmahl and Orralll lll986f) (left) and 
iParenti and Viall l l2007f) (right). 



DEM gives a direct information on the gradient of PCTR of the threads. The second 
model pictures prominences as made of isothermal threads , each of them having a dif- 
ferent temperature in the observed UV-EUV range (e.g.. [Poiogalll993 ). In this case, 
most of the EUV emission comes from the threads, and the DEM relates to the amount 
of threads at a given temperature. 

Quiescent prominences. The DEM derived from ISchmahl and Orralll | 19861 ) was ob- 
tained including the effect of H Lyman continuum absorption (Fig. [G]). They tested 
three geometric models (threads with hot sheaths and cool core; isothermal threads; 
threads with temperature gr adient along the longitudinal magnetic field) and did not 
arrive at a conclusive result. lEngvoldl (|l988l ) showed a steeper gradient of the DEM 
below 10^ K than at higher temperatures. He concluded that extra, unknown, heating 
mechanisms should play a role, in addition to thermal conduction, in order to balance 
the observed radiative losses. 

It is important to constrain the DEM calculation with as many lines as possible, 
distributed over a wide temperature interval. Since the early measurements mentioned 
above, the development of multi-wavelength instruments at high resolution has allowed 

us to obtain more reliable results. 

Examples of these results are those of IWiik et al.l (|l993h . who use d lines from the 
rocket instrument High Resolution Tele scope and Spectrograph (HRTS. bartoe and Bruecknerl 
I1975I ) over the range 4.3 < logT < 6.1: ICirigliano et al.l ||2004) who combined SUMER 



and CDS data to better constrain the results; and Parenti and Via] (|2007l . see our 



Fig. [B]) who used more than 60 SUMER li nes (lying abov e the H Lyman continuum 



head) from their prominence spectral atlas (jParenti et al. II2OO4, 2005b,). These resuhs 



show that the DEMs in prominences are very similar to those of the quiet Sun. The 
main differences are to be found in the amplitude of the gradients at low and high 
temperatures and on the position of the minimum between them. 

The minimum temperature at which the DEM can be derived is limited by the 
optical thickness of the plasma, because the total intensity of the line is modified. The 
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DEM is correctly derived for optically thin emission, while lines at low stages of ion- 
ization or n eutrals may be affected by opacity and by the hypothesis behind the DEM 
derivation ( Avre tt 20071 ). Even if the effect is less str ong than in th e quiet Sun, some 
absorp tion is still present ( Schmahl and OrralllllOsil ). For example. IParenti and Vial 



20071 ) put a limit to about 4.2 in logT. 

For te mperatures smaller than l og T ^ 5—5.5, the DEM de c reases r oughly as T"^ — 
(see iKucera and La ndi' '2006'. for discussion). l Engvoldl dlQSSl ): ICirighano et all 



( 20041 '): and lParenti and Vial (2007) found the minimum of the DEM between logT 



5.1 — 5.2, which is a smaller temperature th an what was found for the quiet Sun and 
for activated prominences (log T = 5.4 - 5.6, iKucera and Landil[20oi : IParenti and Viall 
[^007; Kucera and Landi 2008). 

Discrepancies among the results stimulate further investigation. Beside the uncer- 
tainties arising from the different methods used, the lines selected and the assumed 
theoretical atomic physics, such differences may have a physical origin. It may be that 
the differences are linked to t he dif ferent levels of activity in the prominence. Also, 
IChiuderi and Chiuderi- Dragol (|l99lh showed that the DEM is proportional to cos 6 



(the angle between the direction of the magnetic field and the gradient of tempera- 
ture, which enters in the thermal conduction term in the energy equation), so that 
the derived thickness of the PCTR change s with the angle of view of the observed 
prommence. The modelling of lHeinzel et al.l (|2005l ') is consistent with this picture. 



Another aspect worth mentioning is that the derived DEM generally includes that 
of the background and foreground emissions. Studies addressing this are mentioned in 
the next section. 

Activated and erupting prominences We report on a few examples of temperature mea- 
surements in activated and erupting prominences. The int ensification of activity in 
prominences involves higher velocity fields and temperatures. [Kucera and Landil (l2006l) 



studi ed a jet in an activated prominence on 17 April 2003 using TRACE i Handv et al.l 
Il999l '). SUMER, and Ha from the Kanzelhohe Solar Observatory. Their DEM analy- 
sis on different parts of the jet shows curves with similar profiles below logT = 5.4 
(the local minimum), but with varying gradients and coronal components for higher T. 
The low temperature part of the DEMs is following a power law function with index 
equal to —3, although this may be artificially steepened by a gap in the temperature 
coverage in the DEM. This value is steeper than that found for quiescent prominences 
and active regions. These DEMs are one of the few existing in the literature where a 
background emission has been removed from the data. This is a very important point 
(see the authors' discussion). It is difficult, generally, to have good data for the instru- 
mental stray light and real background emission. The authors are more uncertain on 
the background subtraction at coronal temperatures than in the low TR part. 

iKucera and Land j (|2008l ') also studied an erupting prominence with SUMER. Ap- 



plying again the background subtraction, they calculated the DEM in various parts of 
the prominence. They found the TR DEM to be similar to that of the non-erupting 
prominence, also with a minimum at logT = 5.4. In addition, they found significant 
hydrogen Lyman continuum absorption (at wavelengths below 911 A) which varies 
inside the structure. In particular, the effect was less important near the base of the 
prominence and in a kink region, which they estimated to be areas of mild heating. 

Interesting results come from attempts to measure erupte d filaments in the high 
corona obtained with UVCS data. Applying the EM technique. ICiaravella et al] (2000) 



found a flat EM distribution in the range 4.6 < logT < 5.5 at 1.7 Rq. In one case 
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ICiaravella et alJ (120031 ') found temperatures up to 1.6 x lO^K in the filament upper edge, 
through detecting Si XII and Mg X emission. This could be the result of prominence 
plasma being heated and ionized as it merges with the ambient corona. 



3.3 Electron Density and Gas Pressure 

Here we focus on relatively direct measurements of density and pressure in the PCTR. 
More exte nsive reviews of p revio us prominence dens i ty me asurements can be found in 
reviews bv lHiravamal (|l990l ') and iTandberg-HanssenI (l995l ). 



3.3.1 Electron Density from Line Ratios 

Lines formed through allowed transitions, such as those we have described in the pre- 
vious section, have an intensity which depends on the electron density squared 

In forbidden and intersystem transitions, the upper level is a metastable (m) and 
its radiative decay rate is very low (1 — 100 s~^). This means that if the electron density 
is high enough, the coUisional de-excitation can compete with spontaneous emission 
in depopulating the level. This happens at different densities, depending on the actual 
lines used in the analysis (see an example in Fig. (Tjl . Note that a high rate of collisions 
could also result in a level population of m comparable to that of the ground level. 

The properties of the metastable level populations can be used to diagnose electron 
density. Let us consider the general case of two lines from the same ion formed through 
the transitions j — >■ p and m — >■ fc . Assuming the same emitting volume, from (|6|) their 
ratio may be written as: 



/(A 

mk) l^m Cjyik ^jg 



(16) 



At low electron densities, only lines arising from the levels excited from the ground 
state g are visible. On the contrary, at higher electron densities, the lines originating 
from the de-population of the metastable level will be also visible. Looking at the term 
ng/Nion in ©, we see that for allowed transitions Ug oc no, while Um oc F{T,ne), 
because of the density dependence of the population of this level. This means that the 
intensity ratio of (|16p is a function of density: 

This variation in density dependence of lines formed from the same ion becomes 
a tool to estimate the averaged electron density of the emitting region: the ratio of 
line intensities from an allowed and forbidden transition is a function only of density 
(assuming a constant temperature in the emitting volume). It should be pointed out, 
however, that this ratio is an averaged value based on the emitting plasma along the 
LOS. 

Since 1995, electron density values have been calculated using SUMER. They have 
an advantage over some earlier measurements, in that the ratios consist of intensities 
of lines from a common ion, reducing possible errors related to abundances and tem- 
perature variations. In prominence studies, for example, a commonly used line ratio 
is the one between the 977 A and the 1174.9 A lines produced by C III (the line for- 
mation temperature is 6.3 x lO** K). The dependence of the theoretical line ratio on 



21 



C III 977.02/1 1 74.933 
1 00 [ 



o 



in 
c 

0) 



1 

^q7 ^qIC ^q11 ^q12 

-3 



Electron Density [cm ] 

Fig. 7 The C III density-sensitive ratio (CHIANTI v. 5.1). 



the electron density is shown in Fig. [T] The comparison of the measured ratio with the 
values shown in Fig. [7] gives a unique solution for the electron density in the density 
range 10*^ - lO" cm"^. 



3.3.2 Results of Electron Density Determinations 

As is clear below, measured densities vary greatly. Some of these variations are likely 
due to differences between the various techniques, but there are also true variations 
between prominences, and within prominences themselves. 

Space-based measurements of ne have been done using EUV observations. As a 
result, they are chiefly focused on emission from lines in the hotter portions of the 
prominence, i.e., the PCTR. Thus it would not be surprising if the density values were 
somewhat different than those measured using emission from Ha. The chief method 
has been to use line ratios, which have the advantage that they are independent of 
the filling factor, although they are based on a number of assu mptions wh i ch ar e 
described above. One such assumption is that, as pointed out bv IWiik et al] l|l993l l. 



the calculations assume ionization balance, which may not be correct in activated or 

erupting prominences. 

Th e electron density can also be determined via an emission measure method. IChad 
( 20031 ) derived density values from TRACE observations of a jet in an active region 



prominence. The emission measure was combined with different models of the source 
volume incorporating a number of threads with Gaussian density profiles about their 
central axes. With a Gaussian e-folding width of 1000 km and a slab of depth 3450 km, 
the on-axis electron density was about 10^'^ cm""^. 
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Electron density values for a number of space-based observations are listed in Ta- 
ble |4l along with the temperatures of formation of the lines used for the diagnostics. 
These temperatures are taken from the original papers, and different assumptions based 
on different ionization equilibrium calculations result in different values. Like the values 
derived from optical data, the values from EUV line ratios vary over a significant range. 
Most, but not all, values derived from EUV line ratios are in the range cm 
(the lowest values being mostly related to quiescent promine nces). The densi t y val - 
ues also vary over the prominence. In quiescent promi nences, Cirigliano et al.l ( 200' 



measured values from 1 x 10 ^ to 3.2 x 10^ cm ~^, and IParenti and Viall (|2007h from 
5.75 X 10® to 3.6 X lO'' cm'^. IWiik et all (|l997l l measured density values from 3 x 10^ 
to 3 X 10^^ cm~^ in an erupting prominence. This large range of values partly reflects 
the invalidity of most hypotheses made in the respective derivations (techniques of line 
ratios, EM, etc). 

3.3.3 Gas Pressure 

Pressure values for the PCTR are based on the density sensitive line ratios discussed 
above, and are also listed in Table [J] 

More recently, measurements have varied to a greater extent. As with the den- 
sity measureme n ts, th ese include variations within the prominences. For instance, 
I Cirigliano et all (|2004l ) found pressure valu es of 0.02 at the edge (0.06 dyn cm 
at prominence centre), and iMadiarska et al.l ( 19991 ) found values in the range 0.03 



2 



0.31 dyn cm The aim of future observations will be to reduce the one-order-of- 
magnitude range of values within each type of prominence. 



3.4 Small-Scale Prominence Motions and Non- Thermal Velocities 

Spectroscopic measurements of line position (Doppler shift) and width (Doppler width) 
are commonly used to study prominence motions. The PCTR lines are generally Gaus- 
sian, and are narro wer than those of the so lar TR; however, excess broadening has 
also been observed ( Tandberg-Hanssenll 19951 ) . The method to derive thermal and non- 



thermal broadening is described in Sect. 12. Il 

For quiescent plasma, additional broadening may have different origins: the ion and 
electron temperatures may not be in equilibrium, the ionization equilibrium may fail, 
and absorption effects may be present (Mar iska 1992). For a quiescent prominence, 
it is generally assumed that the first two effects do not occur, or can be neglected, 
mostly because of a high enough rate of collisions. The effect of plasma absorption may 
be overcome by a careful a nalysis of the data, with the identificatio n of modified line 
profiles and intensities (e.g. . [chae et al.lll998l : |Parenti and ViafcoOTi ) . Once those lines 
are excluded, measurements of widths larger than the thermal ones give information on 
the unresolved, non-thermal motions inside the prominences. However, careful attention 
should be given to the possible under-estimation of the instrumental profiles. 

A correlation of N TVs with temperature on the quiet Sun has been observed by 
several authors {e.g., IChae et al.l 1 199^ : IParenti and Viall |2007| ). A similar (but not 
identical) behaviour is also observed in prominences. The line width seems to in- 
crease, increasing the total brightness of the line (Wiik et al. 1993; de Boer et al. | j 



IWiik et al.lll999l ). even though this is not always the case jde Boer et al. 199S). All 



199S 



these observational aspects imply that the non-thermal widths may have an origin in 
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Table 4 Electron density and gas pressure measurements from space-based EUV observations. 
The temperatures shown in the table correspond to the formation te mperatures of the li nes 
u sed for the diagnostics in t he c orresponding papers. Refere nces: 1 — iNov es et al l lll972f ): 2 
- lOrrall and Schmahll (!l976l); 3 - |Kieldseth-Moe et al.l (119791); 4 - iMariska et al.l ||1979I'| T5 - 
Poland and Tandbcrg-Hansscn (1983); 6 - Widing et al. (198^); 7 - IWiik et al.l lll993ri; 8 - 
Wiik et al.lll997.): 9 - Madiarsk a et al.i ([1999.): 10 -.Chae L2003 ): 11 - ICirigliano et al.l i206i ): 
12 — iKucer^^id Land i (2006); 13 — IParentiand Viall ||2007|) . Pressure values marked with * 
are based on ricT values quoted in the paper assuming p = 2nckT, with k being the Boltzmann 
constant. 



Instrument 


Ions 


Prominence 


ric 

(10^ cm-3) 


T 
(K) 


P 

(dyn cm~^) 


Ref. 


OSO-4 and 
OSO-6 


C III 


quiescent 


1 


80000 


0.01 


1 


ATM 


C III 


9 quiescent 


0.7-2.7 


90000 


0.17-0.59* 


2 


ATM 


multi- 
species line 
ratios using 
Si III, C III, 
Si IV, O III 


3 quiescent 




56000 - 
67000 


0.04 - 0.08 
(quieter) 

0.07 - 0.22 

(activated) 


3 


ATM 


C III/S IV 
C III / EM 


quiescent 


9.4 - 40 
8.4 - 29 






4 


SMM/UVSP 


Si IV/O IV 


quiescent 


500, 300 
and below 


100000; 
180000 




5 


ATM 


O IV, Ne V, 
Mg VII 


erupting 


8 - 300 






6 


HRTS 


multi- 
species line 
ratios using 
O IV, N IV, 
Si IV, N V 


activated 
quiescent 


100 


100000 


0.28 - 2.8* 


7 


SUMER 


O IV line 

ratio 


erupting 


3-300 






8 


SUMER 


O IV 


quiescent 


1.3 - 14 


170000 - 
200000 


0.03 - 0.31 


9 


TRACE 


EM -1- 
modcl 


active 
region 


7-19 


250000 




10 


SUMER 


C III 


quiescent 


2 


70000 


0.04 (0.02 
edge, 0.06 
center) 


11 


SUMER 


O V 


activation 
in plage 


3 


250000 




12 


SUMER 


C III 


quiescent 


0.6,3.6 


70000 


0.012,0.07 


13 
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Fig. 8 Quiet Sun (left ) and quiescent prominence NTV (right) as function of temperature 
llParenti and Vialll2007h . 



the unresolved fine structure of the prominence (each thread having different plasma 
conditions and dynamics), and/or the presence of waves (Paper II, Sect. 4) and micro- 
turbulence. The dynamics of a prominence, including the non-thermal motion compo- 
nent, increases in activated and erupting prominences. 

Fi gure E] shows the NTV as function of temperature measured bv lParenti and Viall 
(|2007l ) using SUMER data for the quiet Sun (left) and a quiescent prominence (right) 
in 1999. These plots show a similar correlation of the NTV with temperature: an 
increase of velocity from chromospheric to TR temperatures, and a decrease towards 
coronal values. In both cases, the maximum NTV is found at logT = 5.5. The main 
differences between the two distributions are the smaller velocities in the prominence, 
and a different gradient in temperature. These two aspects have been reported by 
several authors. However, th e amplitude of the gradient difference and the absolute 
values of the velocities differ (jWiik et al.lll993l : ICirigliano et al.ll2004h . 

Somewhat differ ent be haviours of the NTV-T distribution can be found in the liter- 
ature. For example. IStellmacher et al.. (2003. ) report the variation with temperature of 
the UV line widths in the form A\/\, without separating the thermal and non-thermal 
components. They noticed a discontinuity in this distribution at around logT = 4.8 
with a sudden increase of the velocity. They attributed this to the different structures 
of the "cooler" and "hotter" plasma. However, the impossibility of resolving the promi- 
nence finest structure leaves ambiguity in this kind of interpretation. These authors also 
studied the widths of hydrogen Lyman lines. They estimated an NTV > 30 km s~^. 
However, these lines are optically thick, and this result should be considered with some 
precaution. 

It is clear that high-resolution measurements of the unresolved motions, as well 
as Doppler shifts, may bring import ant informa t ion o n the prominence fine structures 
and their dynamics. Using SUMER. Iwiik et all (|l999l ) performed a statistical analysis 
of the moments of the N V 1238 A line (log T* = 5.3) in a prom inence. They then used 
the results from the multi-thread model of IWiik et al.l 1 1993l l to infer the properties 
of the structure's fine geometry. They estimated the NTVs as a function of the line 
intensities, and found an averaged value of 23 km s~^. This value, together with the 
measured line shift, was used to infer the number of threads along the LOS (Fig. [9]). 
Their results were consistent with the multi-thread model, where the intensity along 
the LOS is proportional to the number of threads, and the velocity inside the thread 
has a Gaussian distribution given by the value found for the averaged NTV. 
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Fig. 9 Th e num ber of threads in the LOS (crosses) as a function of the line intensity. From 
IWiik et al.l lll999f) 



A variation of the NTV within the prominence (central region versus periphery) has 
been noticed in some cases jMariska et aLlll979l : lde Boer et al.lll99^ : lstellmacher et alJ 



20031 ). but not in others (|Vial et aLlll98d) . This piece of information may bring insight 



on the heating deposition in prominences. The velocities seem to increase towards the 
peripheral regions and to be more marked for hotter lines. 



3.5 Small-Scale Structure 

It is clear from numerous observations that prominences are made up of many fine 
structures, very likely smaller than the best resolution currently achievable (~ 150 km). 
In addition to images showing that prominences are made up of numerous small-scale 
structures, measurements of the filling factor in prominences indicate that these small- 
scale structures inhabit only a small percentage of the prominence volume, while efforts 
to model small-scale prominence structures indicate that there are numerous threads 
along the LOS. 

In addition to being relevant to the understanding of the structure of prominences 
in general, this question also relates to the question of the structure of the PCTR. 
As depicted in Fig. 1101 there are a number of models describing the possible spatial 
relationships between cold and TR temperature prominence material. 

Filling factors can be formulated in various ways. Here we will discuss measurements 
of the volume filling factor which measures the fraction of the volume filled by material, 
fv- Another filling factor discussed in the literature is the filling factor perpendicular 
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Fig. 10 Possible models of the PCTR. a) Each cold core (~ 8000 K) is surrounded by a 
PCTR. b) Individual threads at different temperatures, c) Cold cores in a larger enveloping 
PCTR. Here the properties vary with height (ho). The length of the threads is not specified 
— in some cas es the prominence is modelled as a series of short knots or blobs. Adapted from 
iPoiogal J 19941) . 



to the LOS, fs, which describes the degree t o which the entranc e slit of the instrument 
is filled by the fine structure elements (see lHeinzel et al]|l996l ). This factor depends 
on the instrumental resolution as well as the physical properties of the object. The 
filling factor calculations discussed here assume that the prominence consists of either 
material of uniform density or of empty s pace, which is, of course, a simplification 
which could affect spectral diagnostics (see ljudg3l200d ). 

For the PCTR, the basic method for measuring the filling factor is to divide the 
volume emission measure EM = fy UedV derived from optically thin emission, by the 
square of the density derived fro m a line ratio d i agnost ic (see Sect. l3.2TT|l . This method 
has been used, for instance, bv iMariska et all (|l979l ). who calculated vo l ume filling 
factors of 0.018 — 0.024, and in an erupting prominence bv IWiding et al.l (|l986l ) who 
found a volume filling factor of < 0.023. 

Lyman continuum absorption measurements (see Sect. 15.2")) have bee n used to cal- 
culate the filling factor of neutral hydrogen and helium in the prominence. iKucera et al] 
(|l998l ) used data from CDS to calculate fs <l 0.3. Using an estimate of the prominence 
dimension along the LOS and s tandard values for th e densi ty, this yields a volume 
filling factor fy of 0.001 — 0.1. IChiuderi-Drago et al.l l|200lh used observations of a 
filament by CDS and common estimates for the prominence density to find a volume 
filling factor of cold material of ^ 0.2. 

Given that the fraction of space taken up by prominence material is relatively 
small, how is this ma terial structured? Clearly the structures must be on a small-scale. 
Ide Boer et al. I (|l998h cite a lack of small-scale (2" ) spatial coherence in lines formed at 
temperatures from 10'* — 10^ K as evidence of multi-thermal threads on that scale or 
smaller. 

A number of researchers have approached this problem by modelling the distribu- 
tions of int ensity, shift , and w idth of spectral lines. Using such an analysis on HRTS 
EUV data, IWiik et all (|l993l ') calculated that there were 15 threads contributing to 
the emission from each pixel at 20000 K and about 30 threads at TR temperatures. 
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IWiik et alJ (|l999l ) found a similar result of 15 — 20 threads for N V emission observed 
wit h SUMER. Both works note spatial scales on the order of a few 100 km or less. 

ICirigliano et alJ (2004) based a multi-thread model on multi-line measurements of 
intensity and Doppler velocity. They find that the number of threads increases with 
temperature, with two to eight 50 km threads needed at logT « 4.8 — 5.6, and 35 
threads at logT ~ 5.8. Over the entire prominence, this would mean about 20 low 
temperature and 800 high temperature threads. These thread estim ates ar e consi stent 
with line moment analyses using Ha data performed by Mein and MeirJ (|l99l[ ) and 
IZirker and Koutchmvl l|l99lh who found ~ 10 threads in a resolution element. There 
are non-LTE models which also use multiple threads to model prominence hydrogen 
emission (see Sects. [8l and [9)l . iMein and Meiru y_991i) note, however, that there seems 
to be a discrepancy between velocity-based and density-based thread models, with the 
density-based models indicating that there are many more threads. They suggest that 
the velocity threads may i n fact be bundles of smaller density structures. 

IChiuderi-Drago et all l|200ll ) studied a filament using both EUV (including Lyman 
continuum absorbing coronal lines, as mentioned above) and radio data. They found 
that the data support a model in which cool threads surrounded by TRs are suspended 
in the hot coronal gas . 

iKucera and LandU ( 2008h compared times and locations in an erupting prominence 
in which there was substantial Lyman continuum absorption and in which the absorp- 
tion has disappeared. The change in TR emission between the two regions was about 
a factor of five. This would require at least 2 or 3 alternating layers of cool and TR 
temperature material along the LOS. 

Although the concepts of filling factor and characteristics of small-scale structures 
are connected, the details of how are still not clear. 



4 Bulk Motions 



Prominences have long been known to be dynamic structures, displaying internal mo- 
tions of various ki nds even wh e n glob ally at rest. A number of review s have been written 
on this topic, e.g. . ISchmiedeij (|l988l ) and lTandber g-HanssenI (119951 1 . In this section, we 



discuss internal bulk motions of non-erupting prominences. Oscillatory motions are 
described and discussed in Paper II, Sect. 4. 

A good understanding of fiows (their trajectories, velocities, and thermal proper- 
ties) is highly important as a test for models of prominence formation and stability, 
as they may help us to distinguish between different mechanisms of mass supply, and 
for the closely related issue of the prominence magnetic field structure. These are dis- 
cussed in more detail in Paper II (Sects. 2 and 3). For this, we need to know the detailed 
trajectories of prominence plasma blobs, including the plasma origin, and any change 
in temperature and velocity. How do their motions compare to various magnetic field 
models - especially in the prominence barbs (Paper II, Sect. 2.1.3)? What connections, 
if any, are there between fiows at different temperatures? 

In the last decade, new instrumentation has yielded more information concerning 
fiows. In the visible, the combination of high temporal and spatial resolutions, along 
with Hq spectral information, has allowed new investigations: counter-streaming fiows 
have been seen in filaments on the disk, using high-resolution instruments which can 
provide images with a spatial resolution of ~ 0.2", and a time resolution of ~ 1 minute 
or higher, along with some Doppler information. The SOT instrument on board Hinode 
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allows for equally good temporal and spatial resolution with excellent image stability 
over long periods of time, offering another probe of prominence motions in Ha and 
Ca II H (3968 A) lines. 

In the EUV range, progress has been made using spectrographs (SUMER, CDS, 
and EIS) which can give information in a range of lines formed at chromospheric, TR, 
and coronal temperatures. Higher cadence UV and EUV imaging information from 
TRACE and, at times, from EIT, has been important as well. 



4.1 Measurements of Velocities 

Measurements of velocities in filaments and prominences are not trivial. There are 
two main methods that have been used to infer motions. Motions in the plane of 
the sky are measured by tracking actual features, while LOS motions are detected 
using Doppler shifts or, in some cases, line profile distortions. Sometimes the plane- 
of-sky and line position methods can be c ombined to g o od eff e ct, as with the gr ound- 



oi-sky and Ime position metnods can be c ombined to g o od ett e ct, as witn tne gr ound- 
based Hq observa tions reported by, e.g., IZirker et al. ( 19981 '): iLin et al] (|2003l '): and 



Chae et al.l ll2006ll but this can only be done to a very limited extent in EUV [e.g., 



Kucera and Landill2006l l. With either method alone, we have no direct information 



on the 3D structure. However, in some cases we can make estimates based on the 
knowledge of the orientation of prominence features as viewed over several days, or 
from two points of view using the STEREO spacecrafts. 

Observations of motions in the plane of the sky have the advantage that it is possible 
to pick out actual moving features, although there are cases, especially in optically 
thin plasmas, in which multiple layers of plasma can make the data hard to interpret. 
Feature tracking is also insensitive to evenly mo ving flows . A co mmon variation of 



feature tracking is to use the time slice method ( Lin et al] |2003| . and Fig. Ilip . The 
intensity or velocity along a slice of the image of the filament is followed versus time, 
and the slope of the brightening or darkening in this 2D diagram gives a measurement 
of the velocity of the feature in the plane of the sky. This method can be easily applied 
to high spatial resolution images, e.g., with the Swedish Solar Telescope (SST) or the 
Dutch Open Telescope (DOT) at La Palma, and has also been used with EUV images. 
Another technique for tracking flows in the plane of the sky is local correlation tracking 
(LCT), although this method is not good at isolating fine features. In the case of the 
relatively low resolution EUV observations, the traceable moving features often make 
up the minority of the emission. There can also be problems in correctly tracing features 
in low cadence (~ 10 min) data sets. 

Doppler measurements more reliably give access to steady flows. However, espe- 
cially in optically thin plasmas, they represent an integration of the LOS velocity over 
many different features so that the measurement provides an average value along the 
LOS. We should also note that current space-based EUV spectrometers do not have 
an absolute wavelength calibration, and large uncertainties may therefore arise in es- 
tablishing the reference line position for Doppler shift measurements. For observations 
of filaments on the disk, the signal coming from the prominence must be disentan- 
gled from the chromospheric background. Various te chniques have been developed to 



gled rrom tne cnromospneric background. Various te cnniques nave been developed tc 
do this , principally based on cloud model methods (|BeckerslE963 : iMein et aLlEiii : 



iTziotzi ou 2007). In such models, four parameters are derived for the cloud: the LOS 
velocity V , the source function 5", the optical thickness r, and the line width AXj^. S 
and r are strongly coupled but V can be easily computed (see also Sect. 18.131) . 
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Fig. 11 Images illustrating the use of time slice diagrams to analyze flows in a Ha filament 
observed with the SST (19 June 1998). The left panel shows the locations of parallel slices used 
for analysing velocities in Ha fine structures. At right from top to bottom are the blue wing 
intensities, Doppler shift, and red wing inte nsities along a single time slice. Diagonal features 
show motions along the slice. From llTin et al.. t2003i) . 



4.2 Quiescent Prominences 

Even the most stable quiescent filaments exhibit flows of 5 — 20 km s~^. Recent research 
has focused on variations and motions in thread-like structures that seem to make up 

fila ments. 

IZirker et al.l 1 19981 ) reported fine-scale counter-streaming in Ha both along filament 



spines (Paper II, Sect. 2.1.3) and vertically in prominence barbs. Moving features were 
traced over 10000 to 100000 km at speeds of 5 to 20 km s~^ perpendicular to the 
LOS. The results were obtained by using movies of the filament taken in the red and 
blue wings of the Ha line, which revealed oppositely directed flows. This is consistent 
with an earlier study of an active region filament bv lSchmieder et all (|l99ll ) in which 
the authors report Doppler measurements of possible intermittent counter-streaming 
in filament footpoints with LOS velocities of ±5 km s~^. They used a cloud model to 

calculate velocities of ±15 km s~"'"±15%. 

Later studies bv lLin et al. 1 (|2003l . l2005h and lSchmieder et al. I l|2008l ) found counter- 



streaming in threads in Ha with average transverse velocities in the range 7— 15 km s~ 
and speeds in one instance as high as 30 km s~^. They also report motions of the 
threads relative to one another with typical velocities of 2 — 3 km s~^, and report a 
case of sideways motion in an isolated bundle of threads in a barb which, they suggest, 
is related to ph otospheric motions. Fig. [TT] shows the time-slice techniques used by 
iLin et a j (l20Q3|) to analyse the data. 



Chae et aLr(l2006l . I2OO7I ) used a cloud model to analyse quiescent filament obser- 
vations along the spine and in a barb at five points in the Ha line. They found the 
data consistent with a filamentary structure in which some filaments show streaming in 
different directions with LOS velocities of 15 ± 3 km s~^, while others contain plasmas 
at rest or moving only very slowly (0 ± 3 km s^"'^). 

These observations and others listed in Table [5] are crucially important in under- 
standing the magnetic structure of prominences, especially in the barbs. The material 
seems to be moving both up and down the barb structures in fairly close proxim- 
ity, in such a way that these opposite flows seem to be crossing when seen in the 
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Fig. 12 Dark upward moving feature s observed on the lim b with Hinode/SOT (30 November 
2006) in the Ca II H (3968 A) line bv lBerger et alJ ll2008l') . 



plane of the sky. The downward motions are slower than what would be expected if 
the material were falling unimpeded under the force of gravity. This indicates that 
some process is at work, which provides an upwards force to the cool material (e.g., 
iMackav and Galsgaard 2001). while still allowing it to move up- and downwards fairly 
readily. These observations are a challenge to most common models of the magnetic 
field in prominence barbs (see discussion in Paper II, Sect. 2.1.3). 

The Hinode SOT instrument has made possible images of prominences and promi- 
nence flows on the limb with unprecedented resolution and pointing stability. iBerger et al.l 
liooi) discuss very different types of motions observed in prominences on the limb in 
Hof and Ca II H (3968 A) with SOT. In observations of large, quiescent hedgerow 
prominences, they report two types of vertically moving features. One of these consists 
of bright downfiow streams with velocities of about 10 km s~^ and lifetimes of about 10 
minutes. These bright downflows are sometimes seen to interact with what appears to 
be large scale vortex rotations. They also report turbulent-looking, dark upflows which 
ascend with velocities of 20 km s^^ in the plane of the sky, and have a large- headed 
"mushroom cap" shape (Fig. [12]). 

IChae et al.l (|2008| ) analysed another hedgerow prominence observed with SOT in 
Ha. This prominence showed numerous vertically oriented structures. They were seen 
to move horizontally across the field of view at speeds of 10 — 30 km until a 
point at which they appear to shed downwards moving blobs (see Fig. [T5)) . The bright 
blobs move downwards with acceleration of 0.015 to 0.083 km s"'^, far lower than the 
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'^20km/s 18:27:57 UT- 18:28:21 UT 

Fig. 13 Vertically oriented feature s seen in a hedgero w prominence (17 August 2007) observed 
with Hinode/SOT in the Ho line bv lChae et al.l ((200^). They appear to move both horizontally 
and vertically in the plane of the sky. White arrows point out a newly formed feature and one 
which is falling downwards. 



gravitational acceleration of 0.27 km s~^. The change in direction from horizontal to 
vertical motion results in a vortical-appearing motion. 

These observations of prominences at the limb present a somewhat different picture 
than that presented by Ha observations on the disk, which seem to describe relatively 
direct flows along straight thread-like structures with inclined up- and downflows in the 
barbs, and more or less horizontal motions along the spine. How do these disk obser- 
vations correspond to the complex motions observed in hedgerow-type prominences on 
the limb? Perhaps the vertical prominence motions observed with SOT are providing 
a side view of the flows in barbs seen in filaments on the disk. 

4.2.1 Flows Observed m UV and EUV 

In many ways, observations in the UV and EUV are quite limited compared to Ha 
observations since they have not yet reached the combination of high temporal, spatial, 
and often spectral resolutions available in Ha. However, the insight they provide with 
regard to temperature information makes them important to our understanding of 
bcisic prominence properties. 

Although fine counter-streaming threads cannot be observed directly in the EUV, 
larger-scale counter-streaming flows in prominences have been observed in H, He I, and 
He n Lyman continuum absorption features in coronal lines by TRACE (|Schmieder et al] 
I2OO8I I. Also, Doppler observations com bined wit h modelling indicate that prominences 
do show such behaviour in EUV. jie B oer et al.l (Il9S|S) observed a quiescent hedgerow 
prominence in a number of chromospheric and TR lines observed by SUMER. Doppler 
shifts were generally correlated between many, although not all, of the EUV lines. 
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In general, however, they did not see detailed spatial coherence in the intensities of 
the different lines, leading them to suggest that the emission may come from threads 
formed at different temperatures. 



IWiik et~al](|l999h observed a polar crown prominence in the N V line at 1238 A with 
SUMER. They found the prominence to be quite dynamic, with structures changing on 
time scales of a few hours and modelled the prominence as a collection of tiny threads 
with a Gaussian velocity distribution with a width of 23 km s~^. 

A numbe r of studies o f EUV observations of prominences have described relatively 
fast motions. Imn3(|l999l ') tracked individual knots and brightening in quiescent promi- 



nences seen by EIT at 304 A with a 20 min cadence. He measured speeds in the plane of 
the sky of 10 — 70 km s~^ which increased as the moving features "surged" from one end 
of the prominence as seen above the limb to the other. These motions were most clear 
during the period of prominence formation. Motions continued later in prominence 
development, with the prom inence becoming increasingly complex and filamentary. In 
a later paper, IWan j (|200lh suggested that such jets are related to flux cancellation 
occurring near the footpoints of the prominence barbs. Such jet-like motions have also 
been observed in active region prominences (Sect. 331) • 

iKucera et al.l l|2003[) analysed prominences in quiet and plage regions using CDS 
wide band movies taken in He I, O V, and Mg IX. They found multi-thermal features 
with velocities across the line of sight in the 5 — 70 km s^^ range. Some of the motions 
appeared to be horizontal in nature. For others, the direction in three dimensions was 
not as clear. iKuc cra and Landi (2006) measured a velocity of about 40 km s~"'^in a jet 
like feature in an activated prominence. 



IWiik et all (|l993h observed two prominences using the HRTS instrument. One 



prominence, described as hotter and more active (including a footpoint eruption during 
the observations), was oriented so that the LOS was along its axis. It showed Doppler 
shifts of ~ 30 km s~^ along the edge of the prominence and LOS velocities between 
±10 km in the rest of the prominence. The second prominence in a plage region 
had a more north-south orientation. It appeared cooler and quieter, and showed much 
lower bulk velocities of ~ 3 km s~^. One footpoint, however, exhibited a LOS velocity 
of 50 km s~^ 

One question of interest is whether there is a variation in velocities of plasmas 
formed at different temperatures. A general trend of higher velocities detected in EUV 
observations as opposed to Ha observations can be seen in Tables [S] and [SI which 
summarise many of the results discussed in this section. However, it is not clear 
if this trend is real, or if it is the result of differences in observing capabilities and 
target selection. There have been some efforts to combine observations of Ha and EUV 
motions, but they have been somewhat limited by the lower resolution of the EUV 
data. 

As to compa risons of plasmas observed at different temperatures in the EUV, 
IWiik et al.l ( 19931 ) and iKuccra and L andi (2006 , 20 08 , ) found motions at tempe r atures 



below about 2.5 x 10^ K to be quite similar, and Poioga and Molownv-Horad (Il999l) 



found similar motions in the range lO'' to 10*^ K. On the other hand 



.p..^ o......^. ... .... .....6. ... V... .... ... .... ......... 

([1.998) report variations. Using Doppler measurements, iCirigliano et al 



de Boer et al.l 
\200i ) found 



an increase in velocities with temperature: from 16000 to 40000 K, velocities were 
22 - 28 km s"\ while for 60000 to 400000 K, the mean velocity was 35 km s"^ 
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Table 5 Velocity flows in quiescent filaments and prominences observed in Ha. Acronyms: 
BBSO - Big Bear Solar Observatory; SVST - Swedish Vacuum Solar Telescope, now SST. 



Observation 



Speed 
(km s"-"-) 



Notes 



Method 



Telescope Source 



On Disk 

Spine, barbs 5-20 ± LOS 



Polar crown, 
spine 



Spine, barbs 
Barb 



8 ± LOS 



15 ± 10 ± 
LOS 

5-13 ± LOS 



Counter- 
streaming, 
vertical in 
barbs, horizon- 
tal in spine 

Counter- 
streaming 
horizontal in 
spine 



± to thread ori- 
entation 



Doppler 
time slice 



Doppler 
time slice 



BBSO 



SVST 



Doppler SST 
time slice 



IZirker et al.l 
JlQQST) 



iLinet al.l 
l l2003f) 



iLin et al\ 
J2005f) 



Spine, barbs 10 ± LOS 



Spine, barbs 15 ± 3 LOS 



Barb 

Limb 

Hedgerow 



10 ± LOS 



< 35 ± LOS 

10-30 ± 
LOS 



Intermittent 
counter- 
streaming 

Counter- 
streaming 
+ stationary 

Upwards 



Vertical 
Horizontal 



Doppler Hida 
time slice 



Doppler 
shifts, cloud 
model 

Feature 
tracking 

Feature 
tracking 



BBSO 



SVST 



SOT 



ISchmieder et aP 
l|2008f) 



IChae et al. | 
t2006l. 120071^ 



Ivan BallegooiienI 
{2003 ) 



IChae et al.l 
1I2OO8D 



Hedgerow 



10 ± LOS 
20 ± LOS 



Bright down- 
flow 

Dark upward 
moving feature 



Feature SOT iBerger et al] 

tracking (2008) 



Hedgerow 



2-10 LOS 
2-25 ± LOS 



Horizontal 

Up- and Down- 
flows 



Doppler 
shifts 
Feature 
tracking 



MSDP ISchmieder et alj 
ll201Cri 

SOT 



4-2.2 Cross-Field Diffusion 

Most of the motions discussed up to this point are assumed to be along or with mag- 
netic field lines. Another p ossible type of bulk flow is cross-field diffusion of neutral 
atoms. [Gilbert et al.l (l2002 ll calculated the diffusion times for neutral atoms in a simple 
prominence model to be 22 days for hydrogen and 1 day for helium. This difference is 
used to explain a decrease of the ratio of He to H as a function of heigh t inside stable 
quiescent filaments derived from He I 10830 A and Ha observations by [Gilbert et al.l 
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Table 6 Velocity flows in quiescent filaments and prominences observed in EUV. LCT stands 
for Local Correlation Tracking. 



Observation 



Speed Notes 
(km s"-"-) 



Method 



On Disk 

Barb 

Limb 

Polar crown 
Prominence 
Barb 

Polar crown 
Hedgerow 
Prominence 
Three 

prominences 



Prominence 



Prominence 
jet 



< 6 LOS 

± 10 LOS 
30 ± 3 LOS 
50 LOS 
23 LOS 



±8, maxR:i25 
LOS 

10-40 ± 
LOS 

10-70 ± 
LOS 

5-70 ± LOS 



22±3-28±2 
35 ± 2 LOS 
38 ±4 



Vertical 



Prominence ori- 
ented II LOS 
Prominence ori- 
ented ± LOS 



Width of veloc- 
ity distribution 



vertical 



Mostly horizon- 
tal 



Feature seen in 
Ho and PCTR 



Doppler 
shift 

Doppler 
shift 



Doppler 
threads 

Doppler 
shift 

LCT 



Feature 
tracking 

Feature 
tracking 



logr = 4.2-4.6 Doppler 
logT = 4.8-5.6 



Feature 

tracking, 

Doppler 



Telescope Source 



SUMER iKucera et aD 



HRTS 



'Wi ik et al.l 
1 11993 



SUMER 'Wii k et al.l 
(il999li 



SUMER Ide Boer et"all 
lll998fl 

ATM iPoioga and Molownv-Horaj 
Jl999f) 

EIT IWand llT999l l 



CDS iKucera et all 

j2003f) 

SUMER 'Cir idiano etal] 
{20m 



TRACE, IKuc era and Landil 
SUMER (200l) 



( 20071 ). and thus an overabundance of neutral heliu m in the lower part of the fila- 
ments. However, for more active or erupting filaments. iKilper et al.l ( 2009l l report that 
the variation in the H/He ratio disappears and is accompanied by a general increase 
in absorption in both lines in filaments prior to eruption, suggesting that this injec- 
tion of mass prior to eruption happens quickly enough or mixes the material such that 
cross-field diffusion cannot cause the separation of neutral If and He. 

ISchmieder et al.l (|2010l ) measured the velocity vector in prominence flows observed 
by SOT (prominence shown in Fig. [TJ by combining measurements of transverse veloc- 
ities from SOT images and Dopplershifts derived from Ha profiles observed with the 
MSDP spectrograph operating on the solar tower in Meudon. Although the motions 
appear nearly vertical in the plane of the sky, the velocity vectors suggest that the 
flows are inclined by an angle of 30 to 90 degrees. Perhaps the vertical prominence 
motions (down flows) observed by SOT p r ovide a side view of the flows in barbs seen 
in filaments on the disk. Schmic der et al.l (|20ld l suggest that the dark bubbles rising 
up inside the hedgerow prominences may be due to higher magnetized r egions which 
could correspond to parasitic polarity emergence in the filament channel. ISerger et al.l 
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(|2008l ) mention cross-field diffusion as a possible cause for the motions they observe 
with SOT (Fig. [HI). 

4.3 Flows in Active Region Prominences 

One of the signatures of active region prominences i s their high level of activity, which 



includes easily seen flows along the prominence axis (|Tandberg-H ansscn 1995). A num- 
ber of studies report the formation of active region prominences via abrupt jet-like 
flo ws, oft e n asso ciated with observed activity in the magnetic field. 

IChael ( 20031 ) found fast (80 — 250 km s^^), jet-like and eruptive behaviour in a 



forming active region filament observed by TRACE in the 1600 and 171 A channels. The 
material was often visible in absorption at 171 A, but appeared to be heated with time 
so that it eventually showed coronal temperature emissio n in the EUV. Th ey found that 
the jets were associated with cancelling magnetic fiux. iLiu et al] ( 20051 ) reported Ha 



observations of two active region filaments formed by surges with measured velocities 
of 150 and 180 km s~^ perpendicular to the LOS. The filaments appeared relatively 
stable, and one lasted at least two days before going over the limb. Obs ervations of 



the formation of two active region filaments by Schmieder et al.l 1 2004131 ) revealed 



dynamic process with downward velocities of 20 km s ^, but less frequent H q up- 
flows . Other studies have reported shorter-lived fllaments also formed by surges I ZirinI 
ll976l:lLitvinenko and Martinll 19991 ). 



Chae et al. I ('2000') report motions in an active region prominence seen in Hq in- 
cluding a flow system c onnected to neighb ouri ng loops that in jected material with 
speeds of 10 - 40 km s"^ lOeng et al.l (I2OO2I ) and lLin et ahl (|2008l ) described bulk flows 



in Ha, including counter-streaming, in existing active region fllame nts. Velociti e s were 
generally in the range 5 to 15 km s^^, but as high as 25 km s' " MOeng et all (|2002l ) 
report especially active periods characterised by twist motions and downflows which 
they associate with episodes of nearby magne tic activity in parasitic polarity regions 
seen in magnetograms. Ok amoto et al.l l|2007l ) observed an active region prominence 
in Ca H with SOT. They found ubiquitous horizontal flows. Some exhibited steady 
velocities of about 40 km s~^, while others showed more complex accele ration. 

In addition to the quiescent prominences mentioned in Sect. 14.2.11 iKucera et al.l 
( 20031 ) also analysed a prominence in an active region observed in Ha and with TRACE 



at 1216 and 1600 A. In one case, a horizontally moving feature with a plane of the 
sky velocity between 5 and 40 km s~^ was observed in both Ha and the 1600 A band 
(thought to represent C IV emission), suggesting a source with a cool core and a TR. 



4.4 Summary 

As seen in Tables [5] and [6l there is a significant range of fiow velocities that have been 
measured. While comparing these values, one should keep in mind obvious differences 
such as disk vs limb observations, or cool part vs PCTR. The true orientation of the 
flows clearly needs to be carefully appreciated before attempting a sensible discussion 
of their possible causes and consequences {e.g., enthalpy flux). We come back on these 
issues in Secti onfTTI Observations from the forthcoming Interface Region Imaging Spec- 
trograph (IRIqj) may help in addressing some of the questions outlined above, as this 



ihttp : //iris ■ Imsal ■ com/ , 
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instrument will be able to measure both LOS and transverse velocity components of 
neutral and ionized plasma by combining spectroscopy and imaging at high spatial and 
temporal resolutions. 

5 Mass Determination 

Prominence density, total mass loading, and mass composition are critical values for 
theoretical models of the origin of the prominence plasma (Paper II, Sect. 3). An 
accurate determination of these quantities is not a trivial task and, typically, only order 
of magnitude estimates can be made. Efforts over the past decade have significantly 
improved these measurements to start testing models of prominence support, mass 
draining, mass loading, dynamics within the material, and eruption mechanisms. 

5.1 Mass Estimates From Neutral Hydrogen Density and Geometrical Considerations 

A rough estimate of the mass M of a cool prominence (low ionisation ratio) can be 
made using the following relation: 

M ^ iYHmnV , (18) 

where N-^ and mjj are the mean neutral hydrogen number density and the mass of 
the hydrogen atom, respectively, and V is the volume occupied by the prominence 
plasma. The mass obtained in this way is only a lower limit to the total mass of the 
prominence. Using typical values of ~ 3 x 10^^ cm~'^ and V ~ 10^® — 10^^ cm^ 
(using lower and upper limits for the typical dimensions given in Section[T]), one obtains 
M « 5 X 10^^ - 10^^ g. 

Note that the precise determination of the dimensions of prominences and filaments, 
and in particular their geometrical thickness, is far from trivial. The volume V of the 
prominence in (|18[) can be estimated by making assumptions about the morphology. 
For instance, the geometrical thickness can be determined by tracking an object over 
several days and comparing its appearance when on the disk and above the limb. If 
this is not possible, then the thickness can be set equal to the apparent width of the 
structure. 

More accurate techniques are based on the determination of the hydrogen column 
density at each spatial position within the prominence or filament. One solution is 
to use non-LTE calculations consistently solving the radiative transfer problem in a 
promi nence slab as described in Sect.[6l Let us mention here the work of lHeinzel et al.l 
(|l996l ) who, by combining polarimetric data in the Ha and He I D3 lines on one hand, 
and observed Ha intensities (compared to theoretical ones) on the other hand, ob- 
tained values for the geometrical thickness ranging between a few hundreds km up 
to a few 10** km for 18 measurements in different prominences. They then derived a 
corresponding total column mass of 10~^ g cm~^. A different kind of non-LTE calcula- 
tions can also be made using the so-called cloud model ( Mein ot al.ll 199(3 : iHeinzel et al.l 
1 19991 . and also Sect. I8.l3|) . Using the cloud model with the optical thickness, Doppler 
width, source function, and velocity as input parameters, the absorption in Ha can be 
calculated and compared with observations. However, uncertainties in the ionisation 
fraction of hydrogen result in a large range of values for the total hydrogen number 
density. 
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Photoionization for Prominences with Assumed Composition 

3.5 I 1 1 — \ 1 1 r 1 1 1 r 




1000 



Wavelerngth (Angstroms) 

Fig. 14 Plot of the average pfioto-ionization cross-section per atom/ion in a prominence 
with an assumed composition (45% H", 45% H+, 9% He^, and 1% He+). In this plot, the 
contributions to the total photo-ionization cross-section are separated by the absorbing species: 
those due to (below the blue line, wavelength shorter than 912 A), He'' (between the blue 
and the red line, which starts at 504 A), and He+ (between the red and the black line, which 
starts at 228 A). 



5.2 Determinations from EUV Continuum Absorption 

Prominences are observed in EUV coronal lines as dark structures embedded in the 
brighter background of coronal radiation (see Fig. 1161 right). Since the prominence 
material is much cooler than the corona (where such lines are formed), the usual 
complicated combination of atomic processes is reduced to the relatively simple pro- 
cess of absorption, which is the absorption by the photo-ionization continuum. For 
prominences in EUV, the main absorbers are neutral hydrogen, neutral helium, and 
singly-ionized helium, for which the cut-off wavelengths are 912 A, 504 A, and 228 A, 
respectively. The photo-i onization cross-section is inversely proportional to so me power 
of the photon frequency dKeadv and Kilcreasellioool : I Anzer and Heinzefcoosh : a plot of 
the average cross-section per atom/ion for an assumed prominence composition is pro- 
vided in Fig. 1141 The used composition is based on assumed heliu m abundance and on 
the model-dependent ionization fraction for hydrogen and heliu m i Heaslev and MilkevI 



Il978l : lAnzer and Heinzelll999l : iLabrosse and Gouttebroz3 2004|). Since the vast major- 



ity of electrons are in the g round state in H and He (jMilkev et alll973bDl : lHeinzel et al.l 
Il987l : ILabrosse and Goutte brozc 20011 ) , only the ground state photo-ionization cross- 
sections need to be used in an analysis of the continuum absorption (in general, much 
weaker subordinate continua overlap the resonance ones). Much of the research de- 
scribed below has been conducted to calculate the number density and total mass, and 
to reduce errors related to the foreground radiation, unknown depth of the material, 
and emissivity blocking due to the cavity. 
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Fig. 15 Figures from I Gilbert et a! of a prominence on 20 March 1999 at around 14:30 

UT in EIT 195 A (left), the mass map calculated from CDS 625 A absorption due only to H" 
(center), and the mass map calculated from CDS 368 A absorption (right), which is due to a 
combination of H" and He". In the mass maps, a brighter pixel corresponds to greater mass, 
and both wavelengths have the same scaling in brightness, ranging from zero to 3.0 X 10^"^ g 
arcsec"^. Note the additional mass detected in the lower part of the prominence at 368 A, 
presumably due to absorption from He". 



Th e earliest calculations of the continuum absorption in EUV were done bv lOrrall and Schmahll 
who used the filament observations taken by the Harvard EUV spectrometer on 
ATM. Since then, a large number of cases has been observed by SO HO, TRACE, 
STER EO, and Hinode. Mass measurements were first attempted by iKucera et al.l 
(| 19981 ) using several spectral lines observed with CDS to estimate a column density 
of neutr al hydrogen o f the o rder of 10^® cm~^. Follow-up mas s calcu lations have been 
done by IColub et all (| 19991 ) using TRACE observations , by | Pennl ( 2000) for an ac- 
tive region filament eruption, and bv lCilbert et al.l ( 20051 . l200a ) in a large study of 23 
prominences that took into account the coronal radiation in front of the prominence 
material, correcting an underestimation of the prominence mass, to yield total mass 
values ranging from 1 x 10^^ K to 2 x 10^^ g . Similar value s have been obtained for a 
quiescent filament by iHeinzel et al.l ( 2003lJ ). ICilbert et al.l ( 2010l ) extend the method 
by interpolating the coronal background emission for the entire prominence area, yield- 
ing a mass calculation at each pixel to form a 2D map (Fig. I15|) . Total mass estimates 
of the relatively small prominences analyzed are on the order of 10^^ g, with margins 
of error ~ 20%. 

A related line of study has used multi-wavelength EUV observations to determine 
the overall 3D structure of prominences, with an emphasis placed on the geometry and 
extended area of prominences and their related cavities in the EUV. The early work on 
non-LTE radiative transfer diagnostics found a much larger opacity (sometimes 1 — 2 
orders of magnitude) for the hydrogen Lyman continuum below 912 A as compared 
to the Ha line opacity, which can explain observations of extended filament structures 
in EUV co-aligned to much narrower Ha counterparts ( Heinzel et al.l [ioOlbl . [2003al : 
ISchmieder eraklboosl . l2004al ). This indicates the presence of more filament material 
than can be visible in Ha. Further research utilized these findings to determine the 
vertical extent of material in filaments, their 3D structure, and the amount of emis- 
sivity blocking, which is due either to t he lack of coronal emission in the cavity or to 
the presence of the absorbing material ('Schwartz_et_alJ l200 4': 'Anzer and Heinzelll2005l : 
[s.chwartz et al. 2006; Anzcr et al. 2007). Hcinz cl et al. (12008 ') used observations from 



the X-Ray Telescope on Hinode (XRT. IColub et al.ll2007l )" of the prominence shown in 



Fig.[T]to determine the cavity blocking and isolate the continuum absorption (Fig. ll6|) 
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Fig. 16 Multi-wavelength images obtained with XRT (left) and with TRACE at 195 A (right) 
of a prominence observed on 2007 April 25 at 13:19 UT. The prominence is visible in absorption 
at 195 A but is transparent to X-rays, allowing a determination of the emissivity blocking due 
to the cavity and a more accurate calculation of the continuum absorption. Adapted from 
iHeinzel et al] ||2008|) . 



The column density can be related to electron density and effective thickness using 
non-LTE models (see, e.g., Sect. [87lT2|) . They calculated the optical thickness at vari- 
ous wavelengths and a neutral hydrogen column density of 1 — 5 x 10^^ cm~^. 

Previous d e termi nations of the prominence mass using radio measurements by 
llrimaiiri et al.l J 19951) we r e exp anded upon and compared with EUV measurements 
by Chiuderi-Drago et al.l ( 200 li ). who came to similar conclu s ions about the emissiv- 
ity blocking as the studies described above. [Pel Zanna et al.l ( 20041 ) went further and 
attempted to distinguish between a prominence with a transition region and corona 
below it, or a thread-like PCTR without any quiet-Sun tr ansition region bel ow (w hich 
is discussed more in iKucera and Landil I2OO i). They and Ichiuderi-Dragol also 
calculated the neutral helium-to-hydrogen ratio to be A'^(He'^)/A'^(H'') = 0.1 — 0.2. 

Although the various methods do not lead to large variations in the determination 
of the mass, one should be aware of the uncertainties due to the adopted hypothesis 
and also be aware of the unknown nature of the extended filaments (possibly leading to 
some blocking or absorption) which could play a role in the total mass determination. 



6 Basics of Radiative Transfer 

Observations show that there are central cool parts of prominences (or fine prominence 
structures) surrounded by the gradually hotter PCTR (which is optically thin in most 
transitions). Densities are generally higher in the central parts. In these regions, sev- 
eral spectral lines and also some continua {e.g., the hydrogen Lyman continuum) are 
optically thick and this requires a proper treatment of the radiative transfer. In this 
section, we introduce the basics of the radiative transfer relevant to prominence physics 
and outline its most important aspects to understand how prominence diagnostics can 
be done from the analysis of optically thick lines. 

We will start with the relatively simple ID slab models. These models have proven 
to be very useful for the basic understanding of radiation processes in prominences and, 
moreover, they describe relatively well the global prominence radiation properties as 
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Fig. 17 Sketch of the ID geometry of the radiative transfer problem in the case of a vertical 
slab (prominence) above the limb. Dotted arrows indicate the incident radiation. 



observed at lower spatial resolution. Indeed, the prominences look like vertical plasma 
slabs when observed above the limb with low resolution. On the other hand, high- 
resolution images show that prominences are highly heterogeneous, containing many 
fine-structure threads, blobs or fibrils. This is discussed in Sect. (5] 

The radiative transfer equation determines the modification of the specific inten- 
sity of radiation I{i>) along an elementary geometrical path ds due to absorption and 
emission processes {v is the frequency): 



(19) 



The absorption and emission coefficients are denoted as Xi^ and rji, , respectively. Using 
X as the reference coordinate in a simple ID prominence slab (see Fig. I17|) : 



da; n - 

— — — cos o = a . 
ds 

Defining r,/ as the optical depth at frequency u with: 



(20) 



dTi, — —Xudx 



(21) 



we can express the transfer equation in its standard ID plane-parallel form: 



The source function is then: 



dTu 



Xv 



(22) 



(23) 



The radiative transfer equation must be solved numerically because the source 
function usually strongly depends on the radiation field. Only after the source function 
is specified, the integration of the radiative transfer equation is relatively simple (the 
so-called formal solution). In Sect. [T] we discuss various modern approaches used for 
solving the non-LTE transfer problem under prominence conditions. 

A critical issue in all prominence non-LTE modelling is a realistic determination 
of the boundary conditions for the radiative transfer equation. Since prominences are 
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typically rather low-density objects, the scattering of the incident radiation plays a 
dominant role in determining the source functions. The boundary conditions at the 
surfaces of modelled structures are specified by the radiation coming from the sur- 
rounding solar atmosphere (or possibly from other structures around) for all directions 
and considered frequencies. The most important irradiation comes from the underly- 
ing photosphere and chromosphere, and also in UV or EUV lines and continua from 
the chromosphere-corona TR and the corona. The actual conditions depend on the 
geometry of the problem and on the line or continuum transitions under consideration. 
However, once the source function is determined and we want to compute the synthetic 
spectrum, the formal solution of the radiative transfer equation has to be performed 
along a prescribed LOS, typically the LOS along which we observed the structure. In 
this case, only the incident radiation in the direction of the LOS is considered in the 
boundary conditions of the radiative transfer equation. Observing the prominence on 
the limb, the LOS directed towards the observer contains the coronal radiation emitted 
by material located behind and in front of the prominence. This may be relevant for 
some UV or EUV transitions, but not for 'cool' lines such as the hydrogen La line. 
In fact, the La line shows emission from the corona, but this is quite negligible in 
comparison to values of the prominence source function determined by the incident 
chromospheric and TR radiation. On the disk, when observing the filaments, the situ- 
ation is different. In such case one has to consider the radiation passing along the LOS 
from the solar surface below the filament. This poses a great challenge because we don't 
see directly this region, and thus it is difficult to estimate this incident radiation. This 
is a typical pro blem of the so-called c loud model. For Lyman lines, a possible solution 
is discussed bv I Schwartz et aU (|2006l l. 

Depending on the atomic model considered, the incident radiation for various line 
and continuum transitions were specified by various authors. In most cases, the ob- 
served specific radi ation intensities are used. Com piled data from variou s sources have 
been given e.g., by Gouttebroze et aD (| 19931 ) and iHeinzel et al ] (l2005l) for hydrogen, 
iRudawv and Hcinzcl ( 1992) for hydrogen subordinate continua.'Lab rosse and Gouttebrozel 
(|2001 ) and .Labrosse ct al. (2007a,) for helium. These and complementary data files can 
be obtained from the authors. 



6.1 Opacity and Emissivity in Spectral Lines 

We introduce the notion of the line absorption profile for a transition between the lower 
and upper atomic levels i and j, respectively: 

<P, = —J—H{a,x). (24) 

H{a, x) is the Voigt function, a is the damping parameter a = ai + aj pertinent to the 
respective atomic levels i and j, and x = Au/Ai/j^ is the frequency displacement from 
line centre expressed in units of the Doppler width. The Voigt function is a convolution 
between a Gaussian and a Lorentzian profile. This refiects the two dominant line broad- 
ening mechanisms: Doppler broadening in the line core, and coUisional broadening in 
the wings of the line. The absorption profile is normalised: 



I <j)vdu = I 




(25) 
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The energy absorbed in the line is Ej — Ei — huij, where Eij is the excitation 
energy of a given atomic level, and i^ij is the line centre frequency. 

Using the Einstein coefficients for absorption (Bij), spontaneous emission (Aji), 
and stimulated emission (Bji), the absorption coefficient corrected for stimulated emis- 
sion is written as: 

huij hvij 
Xu = UiBij^-^fjji, - njBji^-^'tp„ , (26) 

and the emission coefficient as: 

rji, = njAji^-^-il}„ . (27) 

In these relations, rii and nj are the atomic level populations. The frequency depen- 
dence of the absorption and emission processes is given by the profiles 0^ and tpi, , 
respectively. In most practical cases one assumes that ^l>i, = and this approximation 
is called complete redistribution (see Sect. I6.3"1l . 



6.2 Line Source Function 

As we have seen, the line source function is defined as the ratio of emission and ab- 
sorption coefficients, i.e., 

= ^ (28) 



niBij<j)u — UjBjiipu 



(29) 



^ P.,(-) (30) 



^iBij TLjBjj^ 



with pijiu) = il)v/(j)v. The term of stimulated emission in the denominator has been 
simplified because of its small contribution. This corresponds to the so-called non- 
LTE situation typical for prominences and filaments. The atomic level populations 
and pijiy) depend on the radiation intensity and thus are coupled to the radiative 
transfer equation. Contrary to that, in LTE the level populations are given by the 
Boltzmann distribution and the source function is simply equal to the Planck function. 
Unfortunately, LTE is of no use in the physics of prominences because of rather low 
plasma den sities. Basic co ncepts of the non-LTE physics are well explained in the 
textbook of lMihalad (|l978l '). 



6.3 Partially-Coherent Scattering in Prominence Plasmas 

In the general case of photon scattering in spectral lines, the emission profile differs 
from the absorption one, and p[v) 7^ 1. The emission profile ipu is expressed in terms 
of the so-called scattering integral J: 

/•oo 

J = / J^'(j)y'du' , (31) 
•/o 

1 f°° 

V"!^ = ^ / Ru'.uJu'Av' , (32) 
•J Jo 
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where J^i is the mean intensity at frequency v' , and R^i is the redistribution function, 
i.e. the probability that the radiation absorbed at frequency u' will be reemitted at 
frequency v. As we have mentioned above, the complete redistribution (CRD) assumes 
that p = 1. In such a case, the redistribution function has the simple form R^i ^i, = 
(piiicfiv and, when inserted into (|32|l. one gets ^})v = The photon frequencies are 
completely uncorrelated in this case. The other extreme is purely coherent scattering 
in which the scattered photon has the same frequency as the absorbed one. 

A realistic situation is partially-coherent scattering, somewhere between purely 
coherent scattering and CRD, and we thus call it partial redistribution (PRD). For 
resonance lines (transitions between the ground state g and an upper level j), one has 
in PRD: 



Ru',u^iRii + (J--i)Riii (33) 



The function Rjj follows from purely coherent scattering in the atom's frame, while 
RlII reflects the complete redistribution in the atom's frame due to elastic collisions 
having the rate Qe- 7 is the branching ratio, i.e. the probability that the coherence 
in the atom's frame is destroyed by coUisional perturbation of the upper atomic state. 
Ryi ,1, is the velocity-averaged (i.e. in the observer's or laboratory frame) redistribution 
function, here also averaged over all directions. The critical importance of PRD for 
resonance lines li ke hydrogen La emit ted by quiescent prominences was first clearly 
demonstrated by iHeinzel et al although the nece ssity to consider PRD for 

prominences was mentioned by several authors before fsee lCram and Vardavaj|l978l : 
lHeinzellll983l l. However, the subordinate lines like the hydrogen Ha line which arise 
between two excited atomic levels can be well described by the CRD approximation 
- the coher ence is partia lly destroyed by the lower-level broadening (this is further 
discussed in iHeinze 

ID and 2D non-LTE models of prominences and filaments, including their fine struc- 
tures (threads), have been constructed using the angle- averaged PRD for the first two 
Lyman lines (for higher members of the Lyman series the coherence effects become less 
important and one can use CRD). This is described in Sects. [8] and [9] Within ID mod- 
els, P RD was also us ed for strong resonance lines of helium (Labrosse and Gouttebroz3 



20011 ) and of Ca H l|Gouttebroze and Hein^l2002l '). The emergent profiles computed 



with PRD may significantly differ from those computed assuming CRD and this sub- 
stantially affects the resonance- line diagnostics. However, the resulting line profiles do 
refiect both the PRD scattering physics and the actual shape of the incident line pro- 
files to be scattered. This can be easily understood by inspecting the scattering integral 
in (|32p . If the frequency distribution of the radiation field is flat enough over the line 
absorption profile, then = (j)^ and we get the CRD case (we call this special case 
'natural excitation'). This shows how critical is the real shape of Ju in combination 
with the redistribution function. Finally, let us mention that the line scattering in 
real heterogeneous prominences is highly anisotropic, and thus the angle-dependent 
redistribution functions should be considered for even more realistic diagnostics. To 
our knowledge, no realistic optically thick line transfer was performed for prominences 
using the angle-dependent PRD. 
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6.4 Formal Solution of the Radiative Transfer Equation in a Finite ID Slab 

Here we consider two kinds of schematic ID slabs of a finite geometrical thickness, 
oriented either vertically above the solar surface and irradiated symmetrically on both 
sides (the case of prominences seen on the limb, Fie:. ll7|l . or oriented horizontally to the 
solar surface and irradiated mainly from below (the case of filaments). More realistic 
geometries will be discussed in Sect. (9] 

The formal solution of the radiative transfer equation gives us the outgoing radia- 
tion intensity at the slab surface (i.e. for r = 0) and in direction fi: 

7(0, ^i) = loir, fi) exp(-r/M) + / S{t) exp(-t/^)dt/^, (35) 

Jo 

where /o(t, /i) is the incident radiation on the opposite side of the slab. 
Assuming a constant source function, we get analytically: 

1(0, m) = loir, /i) exp(-r/^i) + S[l - exp(-r//i)]. (36) 

Two limiting cases are important: 

T < 1 ^ ^[1 - exp(-r/^)] ~ St/i^l (37) 
T> 1 ^ S'[l-exp(-r/^)] iiS. (38) 

Now we can discuss the distinction between prominences and filaments. 

For prominences on the limb, the spectral line is in emission (/q = 0, no background 
radiation considered). A* = 1, and so: 

1(0) = ^[l - exp(-r)] . (39) 

In a special case of an optically thin slab we get: 

1(0) ~ = T]D, (40) 

where r ^ 1 and D is the geometrical thickness of the slab. We use the fact that S 
and 77 are related by (|28p . and r = xD- 

However, we have to remember that Iq from other directions drive the source 
function or, equivalently, -q. This is the case in central cool parts, where the radiation 
scattering is the dominant process determining the source function. On the other hand, 
inside the PCTR the temperature is steeply increasing, emission lines become optically 
thin and the coUisional excitation starts to dominate over the scattering. 

For filaments on the disk, the spectral lines are typically in absorption, and for 
/i = 1 we get: 

7(0) = Iq exp(-r) + S[l - exp(-r)] . (41) 

Since the line source function in central cooler parts is generally controlled by the 
photon scattering, we can approximately write: 

S ^ \h. (42) 

where 1/2 is the dilution factor by which the incident solar-disk radiation 7o has to 
be multiplied because there is roughly only one half of the prominence or filament 
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surrounding from which the incident radiation illuminates it (assuming no radiation 
from the corona). For filaments seen against the disk we compute their contrast as: 

M= l[i + exp(-r)] (43) 

and in two limiting situations we get: 

r « 1 ^ ~ 1 (44) 

(45) 

For an optically thin filament, the line centre contrast approaches unity, and in the 
case of a large optical thickness, it becomes 1/2. This is why we can see the filaments 
as dark structures relative to the background chromosphere. Since r decreases from 
the line centre towards the line wings, we don't see the filaments when shifting the 
narrow-band filter out of the line centre. 

To summarize, the reason why we observe prominences on the limb in emission and 
filaments on the disk in absorption is the following. The cool prominence plasma ab- 
sorbs the radiation coming from the solar disk and scatters it in all directions. Because 
there is no coronal background in 'cool' optical lines, we see on the limb only the scat- 
tered radiation and the line is thus in emission. On the other hand, the chromospheric 
background of the filament is the absorption line which becomes even darker due to 
filament absorption. The radiation scattered in the direction toward the observer rep- 
resents only a small fraction of the absorbed one and thus cannot compensate for the 
absorption. We thus see filaments darker than the background solar surface. 

There is still another specific aspect of the radiation absorption by prominences or 
filaments, which is now attracting more and more attention with respect to prominence 
mass loading: this is the absorption of EUV line radiation by the resonance continua 
of H I, He I and He H (see Sect. [5]). 



6.5 Statistical Equilibrium Equations 

The non-LTE source function is not known in advance and must be computed by 
solving the transfer problem. The atomic level populations, on which the absorption 
and emission coefficients depend, must be computed using the equations of statistical 
equilibrium which replace the Boltzmann equation valid in LTE. A general form of the 
equations of statistical equilibrium is: 

^ = ^ (i?j, + C,, ) - ^ {R,j + ) (46) 

drij ^ drn dn^V 

dt ^ dt dx ' ^ ' 

Rij are the radiative rates, those for absorption and stimulated emission depend on the 
line and continuum radiation field. Cij — ricSlij (T) are the coUisional rates proportional 
to the electron density ric and dependent on temperature T through the function 
flijiT). Note that the coefficients Cij are defined here in such a way that they include 
Tio, unlike the definition ^ in Sect. 13. Il The time-derivative on the left hand side of 
(|47l) splits into the local temporal variations of [e.g., due to time-dependent heating 
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processes) and the divergence of the flux of atoms in state i {V is the macroscopic 
flow velocity of the prominence plasma). Other equations to be used are the charge- 
conservation equation "^Nf^Z^. — ric, and the total particle number (N) evaluation 
which comes from the equation of state for the gas pressure p: 

p = iVfcr . (48) 

Here Z^. is the ionization degree of k-th species and = ^ Nf^ + Uc [N^. is the total 
density of atoms in a given ionization state k). Finally, knowing A^, the electron density, 
and the atomic abundances together with the atomic masses, one can compute the gas 
density p . 

To be more speciflc, we write the line radiative rates in the form R^j — B^jJij 
for absorption, Rji{spont) — Aji for spontaneous emission and Rji{stim) — BjiJij for 
stimulated emission. Then Rji = Rji[spont) + Rji{stim) , and Jij is the integrated spa- 
tially averaged mean intensity weighted by the absorption proflle, deflned by This 
quantity tells us how many line photons are actually absorbed from the mean radiation 
field, owing to the frequency dependence of the absorption coefficient represented by 
the line profile function 0^. 

In the special case of the so-called 'two-level atom model', the equations of statis- 
tical equilibrium can be written simply as: 

"-I-B12J12 + "■1C12 = ^2^121 + "-2-621 Jl2 + ^2(721 , (49) 

where 1 and 2 refer to the lower and upper levels of the transition, respectively. Com- 
bined with the expression for the line source function, we get, after some algebra and 
dropping the subscripts relative to the energy levels, the well-known formula: 

S = {l-e)J + eB^a . (50) 

In the typical case of a UV resonance line {hv/kT ^1), such as the La line of hydrogen, 
t can be expressed as: 

" " (C21 +A21) ■ ^^^^ 

It represents the probability of a photon destruction. At high densities, e reaches unity, 
and the LTE conditions are achieved. On the other hand, at low densities, e is very 
small, and thus large departures from LTE take place. This is the case for solar promi- 
nences. Taking the typical temperature and electron density values as T = 8000 K and 
He = 10^^ cm^'^, we get e ~ 10~^ for the hydrogen La line. This also means that 
roughly 10^ scatterings are needed before the photon is destroyed (thermalized) . This 
example demonstrates the necessity of using the non-LTE approach for prominence 
radiative transfer modelling. A ge neral formulatio n for both lines and continua within 
a multilevel atom can be found in iMihalad (| 19781 ). 

The above formulation assumes a sta tionary state. N on-stationary radiative state 
of prominences was first considered by lEngyoldl (|l980l l. who estimated the relax- 
ation time needed to reach the statistical equilibrium for Lyman continuum photo- 
ionization/recombination. For typical prominence densities, this time is of the order 
of 1 — 2 minutes or even longer for low-pressure structures. Ife finally concludes that 
these radiative relaxation times are comparable with the ti me scales of th e observed 
variations of the fine structure in quiescent prominences (e.g.. lEngvold|[l97^ . However, 
detailed time-dependent non-LTE modelling will be required to determine departures 
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of the prominence plasma from stationary state. In this respect, the temp erature relax- 
ation towards the radiative-equilibrium state was recently considered bv iGouttebrozj 
( 20071 ') who also found rather long time scales (larger than 10^ s). 



7 Multilevel Non-LTE Problems 

The solution of the mutually coupled radiative transfer and statistical equilibrium 
equations (non-LTE problem) for multilevel atoms can be achieved only numerically, 
and various techniques have been developed for this task. They are briefly reviewed in 
this subsection. 

For prominence non-LTE modelling, two types of methods have been extensively 
used in the seventies and eighties, namely the complete linearization method - CL of 
lAucr and Mihalas (1969), and the equivalent two-level-atom method - ETLA described 
bv lAvrett and Loeserl |l987|). The solution of the radiative transfer equation is achieved 
by the Feautrier method (differential form of the radiative transfer equation), or using 
its integral form. All these and auxiliary techniques are thoroughly described and 
discussed in the textbook of iMihala^ ( 19781 ). 



CL means that all equations which have to be solved simultaneously (namely the 
radiative transfer and statistical equilibrium equations, but also others like e.g., mo- 
mentum balance and energy balance equations) are linearized and then solved as a 
system of linear algebraic equations. Because of this, a Jacobi iterative solution is re- 
quired which converges quadratically. The great advantage of this approach is that 
all equations are solved simultaneously, and the change in one variable automatically 
affects other relevant variables. On the other hand, ETLA uses a generalized form of 
the two-level-atom source function and solves the radiative transfer equation for it, 
transition by transition. Then the radiation fields thus obtained are used to iterate 
the statistical equilibrium equations. In this case, the radiative transfer and statistical 
equilibrium equations are not solved simultaneously for all transitions, but only for one 
at a time. In the case of a strong coupling of atomic transitions, this method can fail 
or may have convergence difficulties. 

A CL-based code, developed for stellar atmospheres, was modified for prominences 
i n the form of ID vertical slabs, and ext ensively used in a serie s of pap ers bv lHeaslev et al.l 
(|1974 hydrogen plus helium plasma), iHeaslev and Mihalad (Il976 | . hydrogen plus he- 
lium, magneto- hydrostatic equilibrium and radiative equilibrium) . IHeaslev and MilkevI 
1^78, other species included) , Hcaslcv and Milkcv ( 1983, jxydrogen Lyman continuum 
modelling). CL was also used bv lFontenla and Roviral l|l983r ) who modelled the hydro- 
gen La line. All these studies used the com plete redistribution approach for hydrogen 
Lyma n lines, althou gh Milkev et al.l ( 19791 ) noticed that PRD may play a role for La. 



Later, iHeinzel et al ] (| 19871 ) used their version of CL to study the hydrogen line forma- 
tion in an isothermal isobaric ID slab. In their work, PRD in Lyman lines was first 
used in a consistent way (with the help of intermediate ETLA iterations), which led 
to an excellent agreement with OSO-8 observations of the La line. 

An other prominence code, based on ETLA approach, was developed by Goutte- 
broze (|Gouttebroze and Labross^ bOOCl . and references therein) , and after extensive 
tests against the results from CL, it was used to compute a large grid of 140 isother- 
mal isobaric ID slab models which still represents an important benchmark for other 
pr ominence modellers or for observers. A subset of these 'GHV models was published 
by iGouttebroze et aL I l|l993h . and the whole set is available from the authors upon 
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request. iHeinzel et alJ (Il994l l then used all of those models to generate various correla- 
tions between prominence radiation and plasma parameters, which have proved to be 
very useful for prominence diagnostics. 

An extensive prominence non-LTE modelling was also done by two groups in Kiev, 
namely by Morozhenko, Zharkova, and Yakovkin and Zel'dina. They studied the forma- 
tion of lines of different elements, by using the in tegral method for mult i level a toms. 
Two pioneering attempts are worth mentioning. lYakovkin and Zel'dinal ( 19681 ) have 



first noticed that PRD might be of importance for hyd rogen La modelling in disk fila- 



ments, but their approach was rather schematic. Then, orozhenkg (|l970l ) considered 
a multi-sla b model to accou nt for the prominence fine structures. This was later de- 
veloped by IZharkoval l|l989h who took into account the mutual radiative interaction 



of individual slabs (see also Heinze]|[l989l ). Multi-slab models (sometimes called multi- 



thread models, according to the fact that we frequently observe many fine-structure 
threads inside the prominence body) have also been used, without mutual radiative 



interaction, by iFontenla et al.l ( 19961 ) and by others. Without radiative interaction. 



one slab is first modelled in detail, and then several (usually identical) slabs are put 
together. The emergent synthetic spectrum is finally computed by a formal solution 
along the LOS through all slabs. This is discussed in more details in Sect. 18. Il 

A principally difi'erent approach to solve the multilevel non-LTE problem is based 
on the so-called Accelerated Lambda Iteration or ALI technique, widely used in stellar 
atmospheric modelling since the end of eighties and beginning of nineties. ALI and 
related methods have proven to be extremely efficient, namely for multi-lev el atoms 
with c omplex line and continuum transitions. We recommend the review bv iHubenvl 
(|2003l ) . A specific approach which is now used in some prominence non-LTE codes is 
called the preconditioning o{ the statistical equilibrium equations, and is thoroughly dis- 
cussed by Rybick i a,nd jrlum m cr (199]; .; , 1992) . It was first in corporated into prominence 
codes by iHeinzell (|l995l ) and IPaletorJ (| 19951 ) , and later by ICouttebroze and Labross3 
(|2000l ). In st ellar atmospheric codes, it is used in connection with CL (hybrid code 
TLUSTY of iHubenv and Lanz|[l995l ). or in connection with ETLA (PANDORA of 



lAvrett and LoeseJ 19921 ). The preconditioning of Rybicki and Hummer was named as 



MALI (Multilevel ALI) method. Combined with the necessary linearization to account 
for non-linear terms in the statistical equilibrium equations, it repr esents a robust tool 



tor non-lmear terms m the statistical equilibrium equations, it repr esents a robust tool 
for prominence non- LTE modelling (jHeinzell [l995t IPaletoul [l995l ) . Using purely the 



MALI technique, no problems have been met with the convergence, even when some 
transitions like the hydrogen La line are extremely optically thick (say up to 10^) 
and others are optically thin. However, by adding other constraint equations (such as 
MHD or energy balance), it is possible that one should have to solve the radiative 
transfer equation in the strongest transitions (which are most coupled to the plasma 
thermodynamic conditions) by CL or ETLA, and then combine this with MALI for 
other less critical transitions. This 'hybrid' approach is now used in some stellar atmo- 
spheric codes. Within MALI, further sophistications like var ious types of acceleration s 
or multigrids are now being applied (see the recent work bv lLeeer and Paletoulboosl ). 

Several sophisticated non-LTE codes have been developed and are used for various 
types of prominence modelling and spectral diagnostics. The most important ones 
which are currently used or are of potential use in prominence physics are the following: 
lAS-code of Gouttebroze (ID, 2D, ETLA, MALI, PRD), Ondfejov-code of Heinzel (ID, 
2D, MALI, linearization, PRD, radiation magneto- hydrostatic), and the code of Paletou 
(ID, 2D, MALI, PRD). Of potential use are also the codes PANDORA of Avrett and 
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Loeser (ID, ETLA, MALI, PRD, multielement) or TLUSTY of Hubeny (ID, hybrid, 
PRD, multielement). These codes are generally available to potential user^. 

8 Results from ID Non-LTE Modelling 

One can ask oneself why it is necessary to build and use a heavy computing machinery 
such as non-LTE radiative transfer calculations in order to get some atmospheric in- 
formation through comparison with radiative observables. The answer lies in the very 
complex processes of line and continua formation, and in particular when the solar 
feature (such as a prominence in our case) is optically thick in some lines and continua. 
In this case, the atomic level populations are non-linearly coupled to the radiation 
field which has a non-local character extending to the whole medium. A well known 
example is the Ha line which may be considered as optically thin in some cases, but 
is inevitably linked to the ha and L/3 transitions which are definitely optically thick 
in prominences. Another answer lies in the importance of e.g., the Lyman lines in the 
radiative budget of prominences. Consequently, the direct inversion of physical param- 
eters from the observed radiation (as described in Sects. [5] and [3]) is an impossible task, 
which leads to the necessity of extensive forward modelling through the solution of the 
non-LTE problem (as described in Sects. [6] and [7|). 

In this section, we review the main results obtained from ID non-LTE computa- 
tions. In spite of its intrinsic limitations, so obvious from the point of view of fine- 
structure observations, we put some emphasis on ID modelling for the following rea- 
sons. On one hand, ID modelling is an historical milestone since it was performed at 
a time when observations were limited in terms of spatial and spectral resolution. On 
the other hand, the progress made in the frame of ID allows to shed some light on 
the continuous improvement of modelling in relation with the physics of the radiation 
{e.g., PRD vs CRD) and the availability of new UV line profiles. We start by present- 
ing the non-LTE modelling of the hydrogen lines and continua (Sect. (ST}, followed by 
the helium lines and continua (Sect. [5^ . and then the spectra from ionized calcium 
(Sect. 18. 3|) . We also devote some space to the case of moving prominences (such as 
active and eruptive prominences), as the radial component of the prominence plasma 
will cause a Doppler effect which has significant consequences on the radiation emitted 
by the prominence (Sect. 18. 4[) . Results from 2D non-LTE modelling are reviewed in 
Sect. [3 

In the following, we will refer to the hydrogen Lyman series as the Lyman series, 
unless otherwise stated. 

8.1 Hydrogen Lines and Continua 
8.1.1 Early Observations and Modelling 

The Lyman series was discovered in the sola r chromosphere as early as 1960 with the 
fir st La spectrum ob tained from a rocket bv lPurcell and TousevI ([l96u) and analyzed 
bv lMorton and Widina (1961 '). The spectrum of the L/3 line (more difficult to obtain 
because located below the 1150 A cut-off of magnesium fiuoride coatings) was recorded 

* The IAS prominence code was described in lGouttebroze and Labrossd 1I2OOOI ') and is avail- 
able at http://www.las.u-psud.fr/pperso/pgoutteb/RTC/rtc.html. 
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Fig. 18 La prominence (3 July 1979) observed with the Transition Region Camera (TRC). 
From lBonnet et al.l ||1980D . 



bv lTousevet ID ^M)- No prominence spectrum in the Lyman series was obt ained. 
The attempts to measure it with Skylab and HRTS ( Kieldseth-Moe et ahlEoT^ l were 
hampered by the fact that the La line was too bright for the exposure time used. This 
is no sur prise when one look s at the beautiful La p ictures of the Transitio n Region 
Camera (|Bonnet et al.lll980l . Fig. [TSll, and VAULT dKorendvke et al.ll200ll . Fig. [28)) . 
since the prominence La intensity is about the quiet Sun intensity multiplied by a 
factor of 0.3 — 0.4, roughly the value of the dilution factor. 

At the end of the 1970s, the non-LTE modelling was ahead of the o bservations 
( Heaslev and Mihalasll 19761: iHeaslev and Milke"vlll976l . Tl97i 'l. EarUer work bv lHiravamal 
( 19631 ) and [Poland and Anzeij (|l97lF " had shown the importance of the incident (UV 
and EUV) radiation for ionizing the prominence plasma. The same was valid for the 
excitation of the atomic levels, especially in the Lyman series. It was realized that 
with reasonable ID models, the Lyman series was optically thick, with r up to 10^ 
for La. Such an important opacity combined with the strength of the incident radi- 
ation as compared to the local (thermal) radiation field in these lines was a natural 
explanation for the (predicted) bright appearance of p rominences in the Lyman series. 
Moreover, numerical efforts were already performed ( Morozhenkol [l973 ) to treat the 
multi-slab radiative transfer within a set of tiny structures (such as threads) in the 
frame of a two-level atom wi th ap plication to the Ha lin e . As for the Lyman series, 
lYakovkin and Zel'dinal l| 19641 ) and lYakovkin and Zel'dinal (|l968l ) considered the exci- 



51 



CYCLE 5020 Lyman-a LAST PROFILE 




-1.00 -0.76 



-0.26 0. 0.26 
WWELENOTH [Afii 



CYCLE 5020 Lyman-p LAST PROFILE 




-1.00 -0.76 



-0.26 0. 0.25 

VWWELENOTH [A^ 



Fig. 19 Observed ha and L/3 prominence profiles. From lViall lll982bl'l . 



tation and ionization from levels 1 and 2, derived electron densities (as proportional to 
the square root of the population of level 2) and finally produced emergent La profiles 
for which the frequency variation of the scattering was properly taken into account in 
the core and the wings of the line. 

The first (photoelectric) prominence spec tral observa tions in La along with L/3, 
Ca II, and Mg II fines came from OSO-8 (jViall Il982al . and Fig. [HI. Earlier on, 



iMihalas et al.l ( 19781 ) had published a seminal pap er on 2D non -LTE transfer in an ex- 
ternally illuminated structure. Their code allowed fviall 1 1982bh to compare La, Mg II 
h and k, and Ca II H and K computed profiles to observed ones (see Sect. [9]). A two- 
level atom was assumed and the ionization was treated with the assumption of photo- 
ionization only but the OSO-8 incident radiation was properly taken into account. This 
allowed the author to build a reasonable model and evaluate the 2D (boundary) effects. 

At this time, the issue of PRD h ad already been addressed for the case of stel- 
lar atmospheres I Hummeilll962l . 19691) . but no clear way for handling the problem of 
collisions had been found before Omont et al.l ([l972) . A fir st trea tment in the frame 
of a ID prominence modelling was applied bv lMilkev et all (Il979l') who evaluated the 
infiuence of PRD on the ratio of Ha to La intensities. Heinzel and Via] l|l983ll also 
discussed the infiuence of PRD on in tegrated inten sities of the La line, and compared 
with OSO-8 observations. Later on, iHeinzel et al. 1 (|l987. . HGVl) performed the first 
realistic modelling with emphasis on the La profile whose detailed shape is very sen- 
sitive to PRD (Fig. I20|l . They showed that PRD not only lowered the far wings of the 
La line but also had a visible influence on the near wings, namely a reproduction of 
the incident La profile. 



8.1.2 Computed vs. Observed Line Profiles and Continuum Intensities 
Since then, a series of modelling effor ts was produced by Gou ttebroze, Heinzel, and 



Dy uou tt 

Vial, epitomized by the GHV paper (|Gouttebroze et all Il993l ) which, among other 



observables, provided emergent line profiles. It should be pointed out that the basic 
work of Heasley, Mihalas and Milkey mentioned above built ID models on the basis 
of prominence integrated intensities only (with the assumption of a fiat incident spec- 
trum), and did not provide the spectral signatures in terms of UV line profiles of the 
various models used. There are a few reasons for that: the incident (chromospheric and 
coronal) line profiles were not well known, and the proper treatment of the line radi- 
ation scattering (see Sect.O was not easy. Moreover, there were only a few observed 
prominence spectra to compare model products with. 
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Fig. 20 L« line profiles emergent from a ID prominence model at three angles. PRD in full 
line and CRD in dashed line. Prom HGVl. 




Fig. 21 Lo, L/3 and L7 line half-profiles emergent from a ID prominence model characterised 
by its temperature (6000 K), its pressure (0.05 dyn cm~^), its geometrical thickness (5000 km) 
and microturbulence (5 km s~^). Abscissae are in A. The line profiles are normalised to the 
value of the maximum intensity, given by GHV in a table reproduced here in Table [7] From 
GHV. 



Working on a set of 140 ID models, GHV were able to produce the correspond- 
ing set of observables (especially the Lyman line profiles and continua, see Fig. 1211 
and Table [T]) comparable with new pro files provided by OSO- 8 and by the Ultraviolet 
Spectrometer and Polarimeter (UVSP. IWoodeate et al.lll98 0') on the Solar Maximum 
Mission (SMM) spacecraft. This was done with a proper account of the incident radi- 
ation profiles (tak en from OSO-8 for La an d L/3, and properly scaled for the rest of the 
Lyman series, see ICouttebroze et ahlEgTSh . and a rigorous treatment of the diffu sion 
within t he line p r ofiles on th e basis of works bylH ummor (1969. ') ; Milkev and Mih alasI 
(llOTSal ): iHubenvl l|l985h : and iHeinzel et al.l (|l987n . An important feature of this mod- 
elling was the use of a 20-level plus continuum atom, which allowed to predict a set 
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Table 7 Table of radiative quantities computed from the same model as in Fig. 1211 The six 
columns give: lower level of the transition; upper level; optical thickness of the slab at line 
centre; maximum intensity of the profile in erg cm~^ s~^ sr~^ Hz"'^; integrated intensity in 
erg cm~^ sr~^; full width at half-maximum (in A). From GHV. 
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Fig. 22 Correlation between the Ha intensity and the emission measure (in CGS units). From 
GHV. 



of Lyman profiles more complete than what could be actually observed with the then 
available spectral instrumentation. They also derived some correlations connecting the 
observables to the thermodynamic parameters {e.g., the Ha line intensity to the emis- 
sion measure, see Fig. I22p . Another example is the color temperature of the Lyman 
continuum, which was found representative (at not too high temperatures) of the elec- 
tron temperature, a confirmation of the result of Heaslcy and Milkcy (1983). Of course, 
one should also keep in mind the limits and assum ptions behind the G HV models (ID, 
static, isothermal, isobaric). But these correlations ( Heinzel et al.|[l994 ) provided a uni- 
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Fig. 23 Correlati on between the squa re of the electron density and the density of level 2 (in 
CGS units). From lHeinzel et alj ||1994| '). 



fying thread for combined well-chosen multi-wavelength observations and a valuable 
tool for their interpretation. For instance, the above-mentioned correlation between 
Ha intensity and the emission measure (see Sect. 13.2. ip is simply explained by the 
quasi-constancy of the term "n^jn^ (Fig. I23[) . The two quantities are related to the 
photo-ionization and recombination terms from and to the second energy level of the 
hydrogen atom, respectively. 

These basic ID models describe the central cool parts of the prominence or its fine 

structures reasonably well. A prominence whose physical propert ies vary with altitude 

may b e simulated by a sequence of individual models. For instance. iGouttebroze and Labrossd 
(|2000l ) computed 20 models to simulate a prominence with a vertical extension of 
10^ km, a temperature increasing from 5000 K at the bottom to 10000 K at the top, 
and a pressure decreasing from 0.2 to 0.02 dyn cm~^ . Other parameters were kept con- 
stant. These authors showed how the hydrogen La line profile varies from the bottom 
to the top of the prominence. Near the bottom, the prominence is very optically thick 
in La and the emission process is dominated by the scattering of the incident radi- 
ation. As the temperature increases and the pressure decreases, the optical thickness 
decreases. The incident radiation also decreases as the effect of altitude. Therefore, 
the scattering decreases and the thermal emission of the slab becomes more important 
at the top of the prominence. As a consequence, the line profile is relatively broad 
and flat at the bottom of the prominence (similar to the incident profile). As altitude 
increases, the profile is transformed into a thermal emission profile with sharp peaks 
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and low wings. It is worth stressing that as the altitude of the prominence is changed, 
the boundary conditions of radiative transfer are changed, since the incident radiation 
is diluted with height. 

Using a simple non-LTE radiative transfer code. lwiik et al derived electron 

densities 1 — 5 x 10^*^ cm~'^ from Ha observations of a prominence observed at Pic du 
Midi with the MSD P. This corresponds to a gas pressure of ~ 0.01 — 0.1 dyn cm 
iHeinzel et al.l (| 19961 ) compared Ha observations with results from G HV and derived 
electro n densities in very good agreement with measurements from iBommier et al.l 
( 1994 2.5 X 10^ - 6.3 x 10^° cm~^). This enabled them to obtain values for the geo- 
metrical thickness ranging between a few hundreds km up to a few 10^ km in different 
prominences, corresp onding to a fairly constant total column mass of 10~^ g cm~^. 
iHeinzel et all (|2001al '). using a more elaborate ID code and including a PCTR, derived 
a constant gas pressure around 0.2 dyn cm~^ for reversed profiles and a variation from 
0.12 to 0.04 dyn cm~^ for non reversed profiles from Lyman lines spectra obtained by 
SUME R. This corresponds to slightly higher electron densities than in the lWiik et al.l 
( 19921 ) study. An other approach to derive physical parameters of prominences has 
been developed recently based on the absorption mechanism of coronal radiation by 
H, He I, and He II continua (see Sect. I5.2"|l . The opacity derived from observations of 
EUV lines (TRACE, Hinode/EIS) was consistent with models of prominences having 
an electron density of the order of 10^^ cm^'^ with typical parameters for the temper- 
ature (6000 — 8000 K) and a thickness of 1 — 5 x 10'^ k m. These val ues yield electron 
densities in prominences in the range 10^° - lO" cm ffleinzel et al. 2008). 

Let us now discuss the observational material concerning lines and continua. In 
the eig hties, the only Lyman spectroscopic results came from OSO-8 (La and L/3) and 
UVSP ( Fontenla et al.iri988l . La). There was a debate about the La reversal, a feature 
which is about 0.4 A wide (peak-to-peak) but is perturbed by the geocorona absorption 
when the observing spacecraft is in low orbit. Moreover, contrary to the La profile where 
the geocorona absorption, about 30 mA wide, is consequently resolved with the OSO-8 
20 mA spectral resolution, the L/3 profile, as observed with OSO-8, was obtained with a 
lower resolution (60 mA) which did not resolve the geocorona absorption. This was not 
of minor importance, since non-LTE modelling always provided a strong disagreement 
between La and L/3 intensities, the L/3 computed intensities being much lower than 
the observed ones (HGVl), a feature also present in the solar chromosphere (see VAL 
and FAL modelsf]- 



8.1.3 Fine-Structure of Modelled Prominences 

In order to overcome such a discrepancy between La and L/3 intensities, several so- 
lutions, apart from the graft of a PCTR (see below), were proposed, most dealing 
with the concept o f fine structu ring of prominences and well supported by observations 
as early as 1960 (|Dunnlll96ll l. Nowadays, SST observations show str uctures down 



to O.l " resolution. Various approaches were considered as early as 1978: iMorozhenkol 
(|l978l l introduced the concept of filamentatton degree and considered a ID multi-slab 
prominence where the mutual interaction between (isothermal) threads was take n into 
accou nt. On the other hand, Fontenla and Rovira (1983, 1985), followed by Vial et al.l 
( 1990l ). considered ID threads with a thermodynamic structure determined by a sim- 



^ The most sophisticated computations of iHubenv and Litej JlQQST) led to some (minor) 
increase in the L/3 wings only. 



56 



Table 8 Intensities of La, L/3, and Ho (in erg s ^ cm ^ sr ^) for various numbers of tiireads 
and different models compared with observed OSO-8 integrated intensities. From FRVG. 



Number of Threads 


Ly«(104) 


Ly/3 Lya/Ly/3 


Ha(10'') 


H«/Lyo 




Model A: 


p = 0.02 (dyn cm"^) 






1 

10 
100 


1.08 
1.47 
2.04 


587 18 
805 18 
1100 18.5 


0.02 
0.21 
1.99 


0.02 
0.14 
0.97 




Model B: 


p = 0.05 (dyn cm~-^) 






1 

10 
100 


1.35 
1.8 
2.47 


1780 7.5 
2419 7.5 
3290 7.5 


0.06 
0.56 
4.81 


0.04 
0.31 
1.95 




Model C: 


p = 0.1 (dyn cm~^) 






1 

10 
100 


1.76 
2.34 
3.12 


4352 4 
5887 4 
8003 4 


0.13 
1.28 
9.16 


0.08 
0.55 
2.94 




Model D: 


p = 0.2 (dyn cm~^) 






1 

10 
100 


2.62 
3.43 
4.42 


10560 2.5 
14240 2.4 
19350 2.3 


0.33 
2.99 
14.9 


0.12 
0.87 
3.4 


Model with p 


= 0.1 (dyn cm 


and Cold Isothermal Core (35 km) 


1 

10 
100 


1.71 
2.94 
3.86 


4127 4.1 
5670 5.2 
7554 5.1 


0.33 
6.29 
18.3 


0.19 
2.14 
4.74 


Model with p 


= 0.1 (dyn cm 


^) and Cold Isothermal Core (45 km) 


1 

10 
100 


1.97 
3.01 
3.86 


4368 4.5 
5563 5.4 
7362 5.2 


1.01 
7.5 
19.2 


0.51 
2.49 
4.97 


GHV Model with p = 


0.1 (dyn cm-2) and T 


= 6000 K 




D = 200 km 
D = 1000 km 


1.63 
2.48 


50 326 
66.5 373 


1.77 
7.15 


1.1 
2.9 


Observed 


OSO 8 


2.8 - 3.6 


440 - 550 65 


1-10 


0.25 - 2 



pie energy equation, including conductive energy transport and radiative l osses. For 
a discussion of the pros and cons of the two approaches, see iHeinzell (|l99(]| ) who de- 
termined heuristically the effect of the mutual radiative interaction upon the exact 
shape of the La, L/3, and Ha profiles. This then led to a systematic modelling of 
threads where ambipolar diffusion (i.e., the diffusion of neutral hydrogen atoms with 
respect t o the protons) betwe en the cool core and the ambient corona was taken into 
account dFontenla et all 1 19961 . FRVG). Also note that PRD was taken into account. 
Agreement between some models and OSO-8 observations was obtained but the high 
necessary number of threads (about one hundred along any LOS) could be questioned 
(see Table [8]). 



57 



PROMINENCE 5 June 1997 file:1 pix=eO (L 4) 
■ 0.3 




923.0 923.5 
Wovslongth (A) 



920.5 921.0 
Wfl^elen^th (A) 



Fig. 24 L4 to L9 prominence p rofiles at a given positi on (solid line) and averaged over the 
SUMER slit (dashed line). From lSchmieder et al.l Il999l) . 



Along the classical approach of ID monolithic (slab) models, another improvement 
was performed with the graft of a PCTR on top of the isothermal isobaric models 
( Heinzel et al.|[T988l ). This was an important step towards the modelling of realistic 
structures while keeping the full complexity of atomic physics and radiative transfer. 
Even though new data in the Lyman series was not expected until SOHO, the modelling 
(and its Lyman radiative signatures) made impressive progress (see Sect. [5)). 

Another path was followed bv lAnzer and Heinzell (|l998l ) who noticed a strong dis- 
crepancy between the physical parameters of the Kippenhahn-Schliiter model and the 
parameters derived from non-LTE analysis. Later on, this led them to develop a 2D 
multi-component model (Sect.[9ll whose spectral outputs were compared to incoming 
SUMER observations of the Lyman series. Rather rapidly, the Lyman series profiles 
available in prominences a nd fil aments were confronted to model predictions. From 
SUMER, ISchmieder etal] l|l999l ) compared the observed L4 to L9 profiles in a quies- 
cent prominence fFig. l24p with the profiles provided by three different classes of models: 
the (isobaric, isothermal) GHV models, the filamentary FRVG models, and finally a 
superposition of GHV-type models with very small thickness which represent an ele- 
mentary filamentary structure. One can repeat two of the conclusions of the paper: 1) 
all Lyman lines seem to be formed at the base of the PCTR, and 2) some temperature 
gradient corresponding to a PCTR is needed to explain the behaviour of several higher 
Lyman lines. It should be pointed out that the modelling used the up-to-date values 
of incident Lyman intensities provided by SUMER. 
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Similar observations had also been obtained in a filament bv lSchmieder et all (|l998l ) 
for L4 to L7 lines. A ID modelling including various types of PCTR was performed, 
which validated the concept that a thin PCTR is necessary in order to obtain Lyman 
profiles which are not too reversed. However, the L4 profile was not well reproduced. 
L4 to L8 line pr ofiles were observed in the pre-eruption and disparition brusque phase 
of a filament bv lSchmieder et al.l ( 2000l '). and the authors found an agreement between 
the velocities derived from the Lyman profiles and the velocity of the front of the 
associated CME bright loop. In the frame of the Joint Observing Program 1070, three 
quiescent prominences were observed with SUMER from LI to L9 together with o ther 
space and ground-based observatories, including the Ha line taeinzel et al.| [2001al '). As 
shown in Fig. 1251 two classes of profiles were obtained: deeply reversed and unreversed 
ones (a reminder of the differences between OSO-8 and UVSP La profiles). In both 
cases, the integrated intensities are similar. Note that because the La profile is obtained 
on the attenuator part of the SUMER detector, its detailed multiple-peaked profile is 
not reliable - but the total intensity is correct. Also note that lStellmacher et al.l (|2003l ) 
observed up to L19. 

The modelling performed along the GHV type of models (isothermal and isobaric) 
proved to be unsuccessful in reproducing the two classes of observed profiles, especially 
the unreversed ones. The addition of a PCTR which nat urally followed, combined with 
a magneto- hydrostatic Kippenhahn-Schhiter-type model ( Heinzel and Anzeijl999l l. proved 



to be decisive in predicting the two types of profiles, provided the angle between the 
LOS and the field lines were taken into account. The results from lHeinzel et aL I (|2001ah 



can be summarized as follows: the PCTR across the magnetic field is very thin because 
of the strongly reduced perpendicular conduction, while the PCTR along the magnetic 
field, being governed by strong parallel conduction, is rather extended. Consequently 
unreversed profiles correspond to observations along the field lines (where the PCTR 
provides line centre photons) and, on the contrary, deeply reversed profiles correspond 
to observations across the field lines (where the cool prominence core i.e., the dips, is 
the only contribution). The whole piece of work which then followed essentially dis- 
cussed the various signatures of such models including the possibility of adding some 
of them in order to reproduce the observed filamentary fine structure. This led to 
the development of a fu ll 2D magneto-hydrostatic model of prominence fine struc- 
ture l|Heinzel et al.ll2005l ). which confirmed the sensitivity of the line profile shapes to 



the direction of the LOS and the magnetic field (see Sect. [5)). This idea was tested 
throu gh the continuous (3 days) observations of a filament/prominence close to the 
limb (jSchmieder et al.ll2007h from L3 to L7 (SUMER) and in Ha. With the filament 



being circular in shape, the authors managed to interpret the daily differences of pro- 
files as the result of the changing angle of the LOS with the magnetic field direction. 
Another use of the Lyman series (at least from L3 to L7) allowed Schwartz ct al. ( 2 0061 ) 
to derive various models fitting different portions of a quiescent filament, to calculate 
the optical thickness at the head of t he Lyman continuum (rgi2), and address the 
issue of the EUV filament extensions ( Heinzel et aLllioOllJ ). The results raised some 



questions concerning the central temperature of the Ha filament (found in the range 
4000 — 6000 K), the ratio of opacities Tgi2/TRa (larger than 50, a surprising high value) 
and the conclusion that EUV filament extensions are, most probably, due to emissivity 
blocking (i.e. a lack of emitting material at mild or coronal temperatures). 



http : //sohowww.nascom.nasa.gov/soc/ J0Ps/jopl07 . txt. 
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Fig. 25 Lq, L/3, and L7 profiles observe d by SUMER in 3 di fferent prominences. Note the 
large asymmetries in the third case. From lHeinzel et al 



Note that iHeinzell l|l995h developed a new code for vertical o r horizontal ID slabs 
based on the MALI method of lRvbicki and Hummed (|l99ll . ll992l ). This code was sub- 
sequently used in all computations made by Heinzel and collaborators. The MALI 
method has proven to be very efficient and stable for all transitions w ithin mu l tileve l 
atoms. PRD is routinely used for hydrogen La and L/? lines, following IPaletoul (|l995l ) 
who first applied MALI to 2D models (see Sect.^. 

Without perturbation by an attenuator on the La channel, IVial et al.l (|2007l ) de- 
rived average ha and L/3 profiles (Fig. I26|l . From the different (spatial) behaviours of 
La and L/3 profiles, these authors worked on a large range of observed La and L/3 
profiles and intensities, and compared this set of values with predictions of ID and 
multi-thread modelling (their Table 3). On one hand, the L/3 intensities seem to be 
better matched by ID isothermal models, and on the other hand, the observed La/L/3 
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Fig. 26 Av erage La and L/3 profiles obtained by SUMER on 3 different regions of a promi- 
nence. From I Vial et al.l j2007l 'l. 



ratio never goes higher than 180 while the thread modelling leads to much lower val- 
ue s. Let us men t ion th at the issue of the number of threads has also been addressed 
by IPoioga et al.l ( 19981 ) who used a different technique (namely the dispersion of line 



intensities) and applied them to La, L/3, L7, and Lyman continuum Skylab data. They 
essentially derived a thread number along the LOS of the order of 20 for (Lyman) 
optically thick lines. 

8.1.4 Electron Temperature From Hydrogen Lyman Continuum 

The hydrogen Lyman continuum in solar prominences is mostly produced through 
photo-ionization due to the chromospheric emission, followed by radiative recombina- 
tion. This emission is related to regions affected by the penetration of the incident 
Lyman continuum radiation. Under the assumption that the electron temperature is 
constant in the emitting region, the source function (Sx) may be approximated by a 



black-body ra diation function (Bx), div ided by the LTE departure coefficient bi of the 
ground state i Noves and Kalkofenlll97ol i: 



Sx = ^ , (52) 

where Tc is the color temperature of the emitt ing plasma. Under promi nence conditions, 

we 



wnere J c is tne color temperature 01 tne emitt mg plasma. U nder promm ence conditu 
Tc can represent the electron temperature T { Gouttebroze et al]ll993l ) . In addition 



may roughly assume that the source function is constant with depth and position, so 
that the emitted intensity can be written as: 



2/ic^ , he 



Here, h, c, and k are the Planck constant, the speed of light, and the Boltzmann 
constant, respectively. In this case T and bi are obtained by fitting (I53|l to the observed 
continuum intens ity. Note that a variation o f T with depth can be derived with the 
technique used bv lVernazza and Noved (|l972l ) i n the chromospheric context. Only a few 



measurements by SUMER have been reported. l Ofman et al.l ( 19981 ) found a large range 
of temperatures (from 5000 K to 15000 K), the largest values bei ng probably re l ated to 
the disappearance (and associated heating) of the prominence. iParenti et al. 1 l|2005al ) 



found different temperatures (8280 K and 7560 K) in different parts of the observed 
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Fig. 27 Lyman co ntinuum specific inten sity obtained by SUMER in two different regions of 
a prominence. From IParenti et al.l i2005al l. 



prominence. The authors note that the uncertainties are of the order of 250 K, mostly 
due to the difficult photometric calibration. They also notice s ome variation of the 
slope with wavelength fFig. I27p. which can be of solar origin (see lVernazza and Novej 
I1972I '). It is obvious that such measurements should be repeated on various parts of 
prominences and on a large set of prominences. It also should be reminded that this 
technique is only valid well below 15000 K, on one hand, and that it also underestimates 
T when the layer is too thin (i.e., optically thin in the Lyman continuum). 



8.1.5 Filament ModelUng 



There are various ways for modelling filaments. Filament modelling in 2D geometry 
is discussed in Sect. [51 In ID, either one takes the usual geometry of a prominence 
and tilts the layer by 90 degrees: its lower boundary receives the whole chromospheric 
(and part of coronal) radiation. Or, one considers the filament as part of a full ID 
atmospheric model starting from the photosphere and ending in the corona. In the 
latter case, the boundary condition for the incident coronal radiation is more difficult 
to define. A nother approac h uses the so-called cloud model, proposed by Beckers in his 
Ph.D thesis (|Beckerslll96i ) in order to deal with structures such as spicules, surges and 
filaments lying above the photosphere and the lower chromosphere, and seen against 
the solar disk. It has been the object of many develop ments, the latest ones being 
summarized in Mein et al. I (|l996l ') - see also the review bv lTziotzioiil (|2007l ). According 
to iMein et~al ] U996'), once the source function is known with respect to depth, three 
important parameters can be derived: the total optical thickness, the Doppler width 
and the Doppler shi ft. The source func t ion c an be derived from non-LTE modelling 
such as per f ormed bv lGouttebroze et al.l ( 19931 ). using the geometry of horizontal slabs. 
iMein et all d 19961 1 successfully applied their method to an Arch Filament System ob- 
served in Ha with the Multichannel Subtractive Double Pass (MSDP) of the Vacuum 
Tower Telescope (VTT). Such a method is especially efficient for observations with a 
wide field of view and a high temporal cadence. 



8.1.6 Summary 

Both observations and non-LTE modelling of Lyman lines have made a tremendous 
progress during the last few decades. The synergy between observations and modelling 
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Fig. 28 Composite pi cture of prominence and filament observed in La with VAULT. From 
iKorendvke et a.1.1 ||2001| 1. 



has come from challenging observations in a variety of structures, the completeness 
of ID modelling shedding light on radiative signatures, and the move from ID to 2D 
and multi-thread modelling. The proof of a necessary PCTR where the various Lyman 
lines are (partly) formed at various temperatures is not a minor feature, even if the 
issue had been raised more than twenty years ago. So, the picture of a bundle of 
magnetostatic threads with a low central temperature (< 7000 K) and an overall low 
pressure (< 0.02 dyn cm~^) seems to be the right one. 

However, one should realize that the range of line profiles and consequently physical 
parameters remains very large. Moreover, the number of threads assumed in the various 
modellings (larger than 10 along the LOS) implies that observations are insufhcient for 
resolving them, in spite of progress in spatial resolution, as shown by the VAULT 
pictures (Fig. I28|) . In the case of filament observations, there could also be a limit due 
to the very high opacity, especially in the La line. This means that for a given set of 
observations performed with the present instrumentation, there might be more than 
one solution. One way to go around this problem is to combine the analysis of hydrogen 
lines with lines emitted by other elements {e.g., He - see Sect. I8.2")l . It should also be 
realized that multithread modelling does not so far take into account the radiative 
interaction between threads. 

On the observational side, one should be aware that the first ten Lyman lines can 
only be obtained through, at least, four consecutive exposure times with SUMER, which 
means that observations are never strictly simultaneous. The total SUMER sequence 
for the Lyman lines takes about a couple of minutes, a duration which may appear very 
small in comparison with the overall lifetime of prominences (at least a few days), but 
which is comparable to the lifetime of fine structures, as evidenced by the SOT images. 
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Having in mind the potential offered by high spatial, spectral and temporal resolution 
of Lyman line profiles, further progress can only be obtained with a dedicated fast UV 
spectrograph combined with a filament /prominence magnetograph. 

On the modelling side, a great step has been made with the use of a non-empirical 
(magneto- hydrostatic) 2D model of dips (see Sect.|9]), but fully 3D MHD models are 
certainly necessary for taking into account the complexity and dynamics of prominences 
and their fine structures. 



8.2 Ifelium Lines and Continua 
8.2.1 Early Studies 

After solving the integral diffusion equation for the radiation, lYakovkin et al.l (|l982l ') 
calculated the populations of selected neutral helium levels with quantum number 
n < 4 under low-temperature plasma conditions. They found that the excitation and 
ionization of helium are maximum near the boundaries. Other authors have considered 
different heliu m lines in diff erent geometries in order to compare with specific obser- 
vations. iMorozhenk^ ( 19841 ) investigated the excita tion of singlet helium levels in the 



frame of homogeneous and filamentary prominences. iLi et al.l (|2000l ) computed the He I 
5876 A (D3) emission line profile using the two-cloud model and taking into account 
the multiplet of He L 



However, iHeaslev et al.l ( 19741 ') were the first to deal with the non-LTE hydrogen 



and helium spectra in model prominences self-consistently. Using simple ID homoge- 
neous and static models, they were able to make the first satisfactory comparisons 
with observed intensities of helium lines in quiescent solar prominences. These authors 
removed several restrictive assumptions that were made before (see references in their 
paper). For example, they used a multi-ion helium model atom with 20 energy lev- 
els and treated the radiative transfer problems in detail in the resonance lines and 
continua for neutral and singly ionized helium. Their prominence models consisted of 
isobaric and isothermal slabs. Then, a series of papers following the same approach 



was published in the following years jHeaslev and Milkevl[l97i . Il978l . 1 19831 ) . However 



there were some limitations, such as the use of CRD (see Sect. [6^ in the formation 
of the resonance lines, or the use of frequency-independent incident intensities for each 
line. They presented integrated intensities, line ratios and optical depths for various 
hydrogen, helium, and ionized calcium lines. However, no emergent line profiles were 
shown (see discussion in Sect. 18. 0]) . 

8.2.2 Isothermal Isobaric Models 

Using a helium system with 34 energy levels (29 fo r He I , 4 for He II, and one level 
for fully- ionized helium). iLabrosse and Gouttebroz3 1 200ll ) computed the neutral and 



ionized helium spectra emitted by a model quiescent solar prominence. The geome- 
try was that of a ID plane-parallel slab standing vertically above the solar surface. 
Their computations allowed departure from LTE for all atomic levels, and the radia- 
tive transfer equations were solved for all optically thick lines and continua. PRD was 
included in the calculations of resonance lines. This represented a new step compared 
to the above-cited papers towards a more physical approach. They obtained several 
radiative quantities related to the emergent spectrum. An important point is that they 
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also presented the emergent line profiles. This is particularly useful for the interpreta- 
tion of prominence spectra, especially when one observes optically thick lines such as 
the resonance lines of neutral (He I 584 and 537 A) or ionized helium (He 11 304 and 
256 A) from which the derivation of the plasma parameters is not straightforward. 

These authors presented detailed comparisons with the pioneering work of Heasley 
and co-workers. This showed the importance of considering the PRD effects in the for- 
mation of the resonance lines of H and He, a n d of u sing detailed, frequency-dependent 
incident line profiles. In fact, iHeinzel et al.l l|l987ll had already shown that PRD is 
im portant in calculating hyd r ogen line profiles (see also Sect. I8.1.ip . but the study 
bv lLabrosse and Gouttebrozel 1 200 ll ) also revealed that PRD for the radiative transfer 
calculations in hydrogen lines and continua is important for the subsequent modelling 
of helium lines and continua, since the electron density stru cture inside the prominence 
slab changes between CRD and PRD _|Heinzel et al."l987^. 

Using a grid of 480 models. iLabross c and GouttcbrozQ ( 200 ll ) investigated the ef- 
fects of the slab temperature, pressure, and width, as well as of the He abundance, on 
the mean population densities and on the line profiles and integrated intensities. This 
investigation illustrated the complexity of the coupling existing between the different 
line transitions which results from the coupling between the radiative transfer and sta- 
tistical equilibrium equations. Their results can be briefly summarised as follows: at low 
temperatures/pressures, the main mechanism of formation of the He resonance lines is 
the scattering of the incident radiation. A substantial amount of ionized helium can be 
found at the boundary between the prominence and the corona. This ionized helium 
is produced via the photo-ionization of neutral helium by the incident radiation, and 
it can scatter the incident radiation at 304 A. Consequently, even cool prominences 
can in principle emit in the He II line at 304 A. Fig. 1291 shows population densities of 
several energy levels of He I and He II averaged across the prominence slab as a func- 
tion of temperature. It is interesting to note that the population of the He II ground 
state (level 30 in this figure) can be larger than the population of the ground state 
of He I (level 1) for temperatures just over 10^ K and at a low pressure of 0.02 dyn 
cm~^. At higher temperatures/pressures, coUisional effects become non negligible and 
contribute to the thermal component of the neutral helium lines. However the temper- 
atures and pressures considered in these isothermal and isobaric models were not high 
enough to produce a significant coUisional contribution in the formation of the ionized 
helium lines. The effect of the slab width is usually negligible on the optically thick 
resonance lines: when the line centre is saturated, only the optically thin wings of the 
line p rofile may become brighter for larger slab thickness. iLabrosse and Gouttebro"z3 
( 200 ih also examined the correlations between various line intensities as a function of 
the prominence plasma parameters. They notably revisited the ratio between the He I 
D3 line and the H/3 line of hydrogen (see Fig. I30|) . This study also showed that the 
different sensitivities of line intensities (optically thick vs. optically thin, singlet vs. 
triplet) could be used together with hydrogen lines, to improve the diagnostics of the 
hel ium abundance in promin ences . 

iGouttebroze et al. 1 l|2002l ) considered the superposition of several slabs along the 
LOS (without radiative interaction) and found that it improves the agreement between 
observed and computed integrated intensity ratios in some cases (i5(HQ)/i5(H/3) and 
_E(D3)/_E;(H;g)), but not for E{Ca II 8542)/i;(H/3). 

A co mparison betwee n obse rved and computed H and He I line profiles was pre- 
sented in ILabrosse et al.l ( 20061 ) . Profiles of the L/3, Le and He I 584 A lines were 
obtained by SUMER in a prominence observed in June 2004. The agreement was satis- 
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Fig. 29 Mean population densities (in cm~'^) as a function of temperature for two differ- 
ent pressures (top: 0.02 and bottom: 0.2 dyn cm~^) and one slab width (1000 km). Solid 
lines: ground states of He I (1), of He H (30), and He HI continuum (34). The population 
densities for these levels are divided by 10*. Singlet excited levels represented: ls2p (5, 
dashes) and ls3d (10, long dashes). Triplet levels are: ls2p '^P (4, dots) and ls3d 
(9, long dashes/dots). Level 31 is the n = 2 level of He II (short dashes/dots). From 



iLabrosse and Gouttebrozd ||2001| ) 
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factory, but the non-simultaneity of the acquisition of the profiles of each line prevented 
the authors to make an in-depth study. They concluded that the temperature of the 
prominence central part was ~ 8600 K and the pressure 0.03 dyn cm^^. 



8.2.3 Models with PCTR 

It is necessary to include temperature and pressure variations in the prominence slab 
if we want to understand how these cool structures can exist in the hot corona. Promi- 
nence models with a PCTR usually represent some sort of semi-empirical atmospheric 
model. This means that the temperature structure is determined empirically in order to 
achieve agreement between the synthetic spectra and observations. By observations, we 
mean either the data for a specific prominence or filament or, more generally, some kind 
of canonical data set. However, problems arise with the variations of the gas pressure 
or density. In fact, all non-LTE prominence models so far used are either isobaric, or 
co nsider certain magne to- hydrostatic equilibrium (MHS - see also Paper II, Sect. 2.2). 

lAnzer and Heinzel (fl999l addressed the question of energy balance in quiescent 
solar prominences using ID slab models in MHS equilibrium and with a PCTR. Their 
prominence models included two parts, namely the core of the prominence body with 
the base (cool part) of the PCTR, where the plasma is still optically thick for some 
lines, and the outer part of the PCTR where the plasma becomes optically thin. They 
pr esent pressure and temperature pr ofiles inside the slab. These profiles were then used 



bv ILabrosse and Gouttebrozd ( 200j) and in subsequent pa pers for the modelling of He 
lines (for the inner part of the prominence only). Later on. lHeinzel et all 1 2001a ) used 



the same kind of models and showed that they were able to reproduce different types 
of hydrogen Lyman line profiles and integrated intensities that were observed by the 
SUMER spectrometer. In this respect, the presence of the PCTR in the modelling was 
found to be critical, e specially for the low members of the Lyman series (see Sect. 18. l1) . 
ILabrosse et al. 1 (|2002l ) presented a study of the PCTR structure and its signatures on 



the emergent spectra of hydrogen, helium, and ionized calcium. They computed the 
spectr a emitted by a ID prominence model in MHS equilibrium as in Anzcr and H einzell 
( 19991 ). The radiative transfer equations were solved for all optically thick transitions 



(when the optical thickness is > 0.1) out of LTE. They presented line intensity ratios 
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Fig. 31 i?(10830)/-B(D3) as a function of the optical thickness at the head of the He I ioniza- 
tion continuum edge (504 A). Left; isothermal isobaric models; right: models with PCTR. The 
different curves distinguish the models according to their mean temperature. Solid line: mean 
temperature T < 6000 K; dotted line: 6000 < T < 9000 K; dashed line: 9000 < T < 12000 K; 
dash-dotted line: 12000 < T < 16000 K; dash-dot-dotted line: 16000 < T < 20000 K; 
long-dashed line: T > 20000 K. Note the different opacity scales in the left- and righ- 
hand figures. Here, all isotherma l models have a mean temperature less than 20000 K. From 
iLabrosse and Gouttebrozd ll2004l ) . 



as well as line profiles. It is worth noting here that it is the first time that predicted 
spectra are computed simuhaneously for H, He and Ca II in prominence models with 
a PCTR. They have shown that the presence of the PCTR affects the line profiles 
in different ways, depending on the optical depth and the region of formation of the 
spectral lines. For example, the inclusion of the PCTR has a dramatic impact on the 
resulting emergent profiles of hydrogen and helium resonance lines. However, it is also 
evident in their study that through the radiative coupling between optically thick and 
optically thin lines, the intensities of the latter are also affected by the inclusion of the 
PCTR in the models. 



ILabrosse and Gouttebro^ (|2004l l investigated the He I triplet line intensities in 
prominence models with and without a PCTR. This study focused on the emergent 
intensities of the neutral helium triplet lines (such as He I 10830 A and D3 at 5876 A) 
and the level populations of the relevant helium states. The authors found that the 
presence of the PCTR affects the emitted intensities of the triplet lines by reducing the 
impact of coUisional excitations at high temperatures in comparison with the isother- 
mal and isobaric case. A simple study of helium energy level populations demonstrated 
how statistical equilibrium is changed when a transition region is present (Fig. I3ip . 
This points to the necessity of including an interface between the prominence body 
and the corona to predict all emergent intensities whatever the region of formation 
of the radiation is, due to the non-local natur e of the coupling between t h e rad iative 
transfer and statistical equilibrium equations. ILabrosse and Gouttebro"z3 1 2004 ) also 
found a correlation between most of the He I triplet line ratios and the altitude of the 
model prominence. This allowed the authors to solve some long-standing discrepan- 
cies in comparisons between predicted triplet line intensity ratios and observations by 
extrapolating the computations to higher altitudes. 



68 



8.2-4 Summary 



After the initial attempt bv lLabrosse et al] ( 20061 ) ■ there should be more attempts 



to 



compare simultaneously observed H and He line profiles with com puted ones. In this 
respect, the He II Balmer lines in the wings of the H Lyman lines ( Ebadi et al.lboo'gl ') 



provide an excellent opportunity to simultaneously observe H and He spectra (repre- 
sentative of cool and hot material, respectively), and put new additional constraints to 
be taken into account in future modelling efforts. The modelling of the helium spec- 
trum should also be extended to disk filaments. Much progress can be expected from 
the modelling of the helium spectrum in 2D geometry (discussed in Sect. O, such as 



ag 



presented bv lLeger and Paletou {2009, ) and .Gouttebroze and Labrossa (|2009i ). 



8.3 Ca II Spectra 

Apart from hydrogen and helium lines, the optical and IR spectral region contains, 
among others, also five lines of ionized calcium which are frequently observed in promi- 
nences. These five lines are the UV Ca II resonance lines H and K plus three IR lines 
(triplet). The high spatial and temporal resolution images p roduced in the Ca II H 
line by SOT (Fig. I12|) were presented in lBerger et al.l l|2008l ). In SOT movies, highly 
dynamical fine-structure features are recorded, sometimes showing dark bubbles mov- 
ing upward. The question thus arises whether these brightness changes are due to the 
ionization state of calcium or due to other effects. This may be related to prominence 

energy and m omentum balance. 

Early on, iHeaslev and MilkevI ( 19781 ') computed the ratio r of integrated intensi- 



ties of the Ca II IR line at 8542 A to the hydrogen H/3 line. They found that in the 
temperature range 6500 — 9000 K, this ratio is relatively insensitive to temperature, 
practically independent of column mass, and decreasing with gas pressure. This theo- 
retical relation was then used to derive the gas pressure from observations. However, the 
model curve was restricted only to pressures below 0.2 dyn cm~^ . Calcium spectrum 
formation was discussed also by iMorozhenko and Zha rkoval (1 1 9821). Hydrog en and cal- 
cium semi-empirical models were then produced bv lzhan^ and Fane! ( 1987 ) who found 



a reasonable agreement with observed, uncalibrated (i.e. , in relative units) line pro - 
files. In order to extend the rang e of gas pressures used bv lHeaslev and MilkevI (1978 



ICouttebroze and Heinzel l|2002h computed a larger grid of models of the GHV type, 
up to p = 1 dyn cm~^. Their principal result is that at pressures higher than 0.1 dyn 
cm~^, the ratio r is strongly dependent on the temperature: r starts to increase for 
temperatures below 8000 K, and decreases at higher temperatures. This means that 
there exists a temperature bifurcation o f the ratio r for high pressures. Contrary to 
the results of iHeaslev and MilkevI ('l978'). this can explain higher measured values of 



say up to 0.8 l|SteUmacher and Wiehr.,200(i ) which do indicate lower temperatures 
and gas pressures higher than 0.1 dyn cm~^. This has important consequences for 
prominence momentum balance since higher gas pressures mean larger values of the 
plasma-/3 parameter. 

Concerning the non-LTE mod elling, the Ca II UV resonance continua are partially 
driven by hydrogen Lyman lines (|Gouttebroze and Heinzell2002h . This presents a sig- 
nificant complication for the calcium transfer modelling. First, the local radiation fields 
in the Lyman lines must be carefully computed and then used to evaluate the photo - 
ionization rates of calcium UV continua. Note that ICouttebroze and Heinzell ( 20021 ) 
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used PRD for H and K lines, but did not treat details of the cross- redistribution with 
IR tr i plet lines. PCTR and multi-t hreads were also treated later on (|Gouttebroze et aP 
I2OO2I : iLabrosse et allbool l2007bh . Fin ally, it is worth mentioning that the measure- 



ment s of ISteUmacherandWiel^ led to a pretty constant value of r over the 

observed prominences. This indicated a rather constant value of the gas pressure, as- 
su ming a uniform t e mper ature in the cool parts. 

iTziotziou et al.l l|200ll ) used a two-step method to interpret Ca II (8542 A) obser- 
vations obtained with the MSDP at the VTT telescope. The first step involved the 
computation of a large grid of models by a multi-level non-LTE transfer code which 
yields the Ca II line depth-dependent mean intensity inside an isolated, isothermal 
cloud lying above the chromosphere. The second step involves the inversion of the ob- 
served profiles with the grid of computed Ca II profiles. Maps of the temperature, LOS 
velocity, and microturbulence were built inside a quiescent filament. The temperature 
was found to peak at 8500 K, the velocity indicated an excess of blueshift (material 
moving upwards), and the microturbulence peaked at 5 km s~^. This method has not 
been applied to SOHO data because it does not work if the medium is too optically 
thick in the observed line (then the incident radiation component of the cloud model 
is not visible). 



8.4 Moving Prominences: Doppler Dimming and Doppler Brightening 

It is still not clear how solar prominences reach a state of equilibrium. However it 
is known that during the formation stage, significant large-scale plasma motions are 
taking place (see also the discussion in Paper II, Sect. 3). Similarly, mass motions are 
important during the disappearance of prominences. Between these two stages of the life 
of prominences, one may also observe periods of activity with internal fine-structure 
plasma motions (see Sect. |4|. Hydrogen and helium resonance lines are difficult to 
interpret as the prominence plasma is optically thick in these transitions. When the 
prominence is moving radially outwards, the incident radiation coming from the solar 
disk and illuminating the structure is shifted to lower frequencies. Therefore, if mass 
motions are taking place, the diagnostic of the prominence is even more complex. In 
order to improve our understanding of dynamic structures, it is necessary to study the 
effect of motions in pro minences on the relevant lines. Most of this modelling has been 
done in ID (apart from lViallE982bh . and we therefore describe it here. 

Let us start with a simple description of the effect of radial motions on the emitted 
spectrum. Consider a simple two-level atom whose upper level is excited by the radi- 
ation coming from the Sun: the absorption profile of the radiative transition between 
the two levels gets out of resonance with the incident radiation when this atom is 
moving away from (or towards) the Sun, due to the Doppler effect. We namely have a 
Doppler dimming effect when the incident line is in emission, or a Doppler brightening 
effect if the incident line is in absorption (see also iHvder and Liteslll97Gl ). In a more 
realistic situation, an atom has more than two energy levels. Consequently, coupling 
effects take place between the atomic levels. For that reason, a combination of Doppler 
dimming and brightening can occur, just as it happens wh en the coupling between the 



first two excited levels of hydrogen is taken into account (iHeinzel and RompoltlEgSTl : 
[Contikakis ct al. 1997b). Note that this is the same Doppler effect which is also used for 
the diagnostic of the radial component of the solar wind (see e.g. . iKohl and Withbroel 
[19821 ). 
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The main factors determining the effects of the radial motions on the emitted 
prominence spectrum are the line formation mechanisms, including the relative contri- 
butions of coUisional and radiative excitation, and the characteristics of the incident 
radiation {e.g., whether it is in absorption or in emission, or the strength of the line). 
The maximum effect is achieved whe n radiative excitation d omin at es the collisional 
one. M o re discussion ca i i be found in Heinzel and Romp old (|l987f l ; iGontikakis et al.l 
(ll997al):lLabrosse et al' ('2006b', '2007a''^ 



iHeinzel and Rompolt (1987) and Gontikakis et al] (Il997al lbl) computed the hydro- 
gen spectrum emitted by an eruptive prominence. The prominence is modelled as a 
ID plane-parallel slab standing vertically above the solar limb. Its geometrical thick- 
ness along the LOS and its altitude above the limb are free parameters. The other 
free parameters, which define the prominence atmosphere, are the electron tempera- 
ture, the gas pressure, and the microturbulent velocity. The incident radiation field is 
identical on both sides of the prominence. This velocity-dependent incident radiation 
illuminating the prominence slab can be represented by the mean intensity ,Ju: 



Jy = 



_}_ I T A. , '^0 



/o ^V-n',n') dn' , (54) 



where 10(^^1 n) is the specific intensity of the incident radiation, n' is the direction of the 
incident photon, and V is the prominence velocity. Ju is calculated at a given height, 

taking into account the ccntrc-to-limb varia tion s (if any) of the in cident radiat ion. 

As illustrated in Hcinzcl and Rompolt (,19871 ') . [Gontikakis et al.l ( 1997a ). and lLabrosse et ahl 



( 2007al '). he Doppler effect induces a shift in frequency of the incident line profile, and 



the variation of the Doppler shift with direction induces a distortion of the incident 
profile. Figure [321 shows this effect for the hydrogen La and Ha lines, and the helium 
line at 584 A. As the outward velocity increases, the central peak of an incident emis- 
sion resonance line profile is less intense, it is moved towards the red, and the line 
width becomes larger. In order to study this Doppler effect, it is of course necessary to 
use a realistic incident radiation that varies with frequency in the calculations. 

I Heinzel and Rompoltl (|l987h did their calculations in CRD and were able to show 



the main effects of the velocity-dependent radiation field on the integrated intensities 
of the hydrogen Lyman and Balmer lines. However, it is important to compute the 
line spectrum in PRD if we want to compare the calculated line profiles with observa- 
tions, and infer the plasma parameters in moving prominences. PRD differs from CRD 
mainly in the line wings (in the line core, CRD dominates). The line wings themselves 
are greatly affected by the frequency shi ft of the incident radiation, and this yields 
potentially large line profile asymmetries. I Gontikakis et al. 1 (|l997bl lal) were the first to 



inc lude PRD in th e ir mod elling of the hydrogen lines emitted by moving prominences. 
iLabrosse et al.l ( 2007al ') presented the first computations of the helium line profiles 



emitted by a moving prominence. The authors focused on the line profile properties 
and the resulting integrated intensities, and their sensitivities to the prominence plasma 
parameters. They also studied the effect of frequency redistribution in the line forma- 
tion mechanisms. Figure [33] shows the differences between CRD and PRD for the 
profiles of three helium lines (584, 537, and 304 A) at different velocities. Note the 
wing asymmetries at higher velocities which are due to PRD effects. These authors 
o nly co nsidered isothermal and isobaric prominence models similar to those presented 
in lLabr osse and Gouttcbroze (2001). 

The effect of the PCTR on t he H and He lines in a ctive and eruptive prominences 
was presented for the first time in lLabrosse et all 1 20081 '). The pressure and temperature 
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Fig. 32 Mean intensities of the incident La (top left), Ha (top right), and He I 584 A (bottom) 
emission profiles as seen by a moving prominence at the height of 5000 km. Top panels: the 
profi les are displayed for the velocities ranging from — 280 km s~^ ffieinzel and Rompoltl 
Il987l ). Bot tom panel: velocities range from (solid line) to 240 km s~^with a step of 
80 km s- ^Labrosse et al.|[2007ah . 



profiles were taken from lAnzer and Heinzell ( 1999h . [Labrosse et al. I (Eooi) showed that 
the He II 304 A line is strongly dependent on the radial velocity, which can be explained 
by the dominant role of scattering in the line formation. It is thus necessary to carefully 
consider the role of the radial motion of the plasma in the modelling of this line. Other 
lines such as He I 584 A and the hydrogen Ha line are also sensitive to the Doppler 
dimming / brightening due to the radial motion of the plasma at velocities up to ~ 
100 km s~^ (Fig. I34|) . At higher velocities, the absorption profile of the transition gets 
significantly out of resonance with the incident radiation, and the resulting variation of 
the emergent intensity with the radial velocity then only depends on parameters such 
as the strength of coUisional excitations, or the coupling with other transitions. 



8.4-1 Summary 



The study of the variation of the resonance line intensities with the radial velocity of 
the plasma indicates that the Doppler dimming effects are essentially present when the 
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AX/aX„ M/M„ Oy/t^K 

Fig. 33 Differences between CRD (dashed line) and PRD (solid line) for the profiles of three 
helium lines at different velocities: 0, 80, 160, and 240 km s~^from top to bottom. Abscissa: AX 
(distance to line centre) in Doppler units and limited to 10 Doppler widths around line centre; 
vertical axis: specific intensity in ergs~^ cm"'^ sr~^ on a log scale. Model parameters are 
altitude H = 50 000 km, width L = 650 km, temperature T = 6500 K, and pressure p = 0.1 dyn 
cm~^. The thin solid line shows the incident line profile, gradual ly shifting towards the red as 
the prominence is moving radially at increasing velocities. From iLabrosse et al.l ll2007al 'l. 



relative contribution of the thermal emission compared to the scattering of incident 
radiation in the lines studied is low. 

Most of the information on the radial velocity can be obtained through the variation 
of the relative integrated intensity, i.e., by measuring the amount of Doppler dimming 
on the EUV resonance lines of hydrogen and helium. Although this measurement is 
difficult to perform (mainly because one has to be observing at the right place and at 
the right time), space-based instruments such as EIT, SUM ER, and CDS on SOHO, 
EUVI (the Extreme Ultraviolet Imager. [Howard et al.lboosi ') on STEREO, SOT and 
EIS on Hinode, as well as TRACE, are able to obtain a temporal series of intensity 
measurements covering an entire prominence eruption, starting before the radial mo- 
tions take place. In this case, it is possible to measure the dimming in intensity, and 
to compare it with predicted relative intensity curves. 
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Fig. 34 Variation of the He I 584 A (top), He H 304 A (middle), and hydrogen Ha (bottom) 
line intensities (in erg s~^ cm~^ sr~^) as a function of the radial velocity of the prominence 
plasma (in km s~^) for 4 different types of promi nence atmospheres. The p arameter 7 refers 
to the pressure and temperature profiles derived bv lAnzer and Heinzell lliggSTl and is related to 
the temperature gradient in the PCTR: the higher 7, the higher the gradient, and the narrower 
the PCTR. From lLabrosse et al.l ll2008l'l . 
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Detectors with a spectral pixel size around 0.05 A or less {e.g., SUMER) are ideal as 
they allow line asymmetries to be observed. A lower spectral resolution still allows good 
comparison with theoretical line profiles, however it makes the task of finding a unique 
model that fits the observed optically thick profiles more diflicult. In combination with 
the study of the apparent motion of the prominence material brought by imagers on 
SOHO, Hinode, and STEREO, the full velocity vector can be inferred. 

Further modelling is necessary to take into account the fact that the prominence 
plasma is heated during an eruption. In this situation, Doppler dimming will be af- 
fected by an enhanced contribution of coUisional excitation. This is one of the major 
challenges in this type of study: it is essential to make a distinction between a change in 
the thermodynamic state of the plasma {e.g., an increase in temperature and density) 
and a change of the radial velocity. Does a decrease in the intensity of a line formed 
at low temperatures (say below 10000 K) correspond to the Doppler dimming effect 
due to the radial motion of the plasma, to variations of the thermodynamic plasma 
parameters, or to a combination of both? To answer this, other lines are needed to 
perform an independent diagnostic of the prominence plasma. For instance, the inten- 
sities of the He I 10830 A and D3 lines are not very sensitive to the radial velocity 
( Labrosse erai]|2007bh . essentially because of a fiat incident spectrum. Therefore these 
lines can be used to determine the plasma parameters. Another mean of distinguishing 
between variations in intensity due to a change of the plasma state or to the radial 
velocity is to use the full line profiles. The radial motion of the prominence induces 
asymm etries in the line profile which cannot be attrib uted to temperature or pressure 
effects l|Gontikakis et al.lll997al : lLabrosse et al.ll2007al ). 

Future studies will need to address active and eruptive filaments. As these struc- 
tures are seen on the disk, the radial motion of the plasma cannot be inferred from 
imaging measurements, and Doppler shifts cannot always be interpreted in a straight- 
forward manner. The case of eruptive disk filaments is particularly relevant as they 
can trigger Earth-directed CMEs. Again, the question of the contribution of density 
variations or heating of the plasma versus the Doppler dimming effect has yet to be 
resolved. This is a complex issue, since many properties of the prominence plasma will 
vary during the eruption. As we progress towards more sophisticated modelling, we 
will be able to disentangle this. Efforts should also be directed towards 2D modelling, 
as discussed in the following section. 



9 Results from 2D Non-LTE Modelling 



The previous section is focused on ID modelling of solar prominences. Although these 
ID models still represent a useful and computationally efficient approach for a number 
of situations, they do not allow us to consistently study the variation of the radiation 
and the plasma parameters in two dimensions, e.g., along and across the magnetic field 
lines. However, prominence observations carried out both from space and from ground 
reveal the structure of prominences at various angles with respect to the magnetic field 
orientation. Proper interpretation of such observations requires the use of more general 
2D or 3D prominence models. 
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Fig. 35 Sketch of the geometry of 2D prominence models. An infinite dimension is oriented 
tangentially (a), or perpendicularly (b), to the solar surface. Arrows oriented towards the slab 
indicate the incident radiation. 



9.1 2D Models of the Whole Prominence 

The geometry of 2D prominence models is characterized by two finite and one infinite 
dimension which can be oriented tangentially, or perpendicularly, to the solar surface 
(Fig. I35|) . The variation of all quantities takes place within the plane of the cross- 
section perpendicular to the infinite dimension along which the medium is assumed to 
be homogeneous. 

A robust method for solv i ng the radiative transfer in 2D planar geometries was first 
developed bv iMihalas et all (|l978h . It makes use of a Hermitian integration formula 
on ray segments through grid points (short characteristics) and of CRD, and is also 
well suited for solving the velocity-dependent problems. This method was used by 
Ivi5l l|l982b[l to construct a 2D horizontally infinite vertical slab model representing 



the whole prominence structure. In this model, the cross-section dimension along the 
axis perpendicular to the solar surface (height) dominates over the other cross-section 
dimension (thickness). Synthetic profiles of H I La, Ca II H and K, and Mg II h and 
k lines, obtained using the two-level atom approximation, were compared with OSO-8 
observatio ns and a good agreem ent was found for La and Ca II lines. However, the 
method of iMihalas et al.l l|l978l ) is computationally prohibitive an d was designed to 
treat CRD and the two-level atom approximation only. Therefore, lAuer and Paletoul 
(|l994') developed a new method for the solution of the non-LTE radiative transfer 
problem in 2D pro minence geometry wh ich can treat also PRD, important f or stro ng 
resonance lines (see iHeinzel et al.l Il987l ') . The method of lAuer and Paletoul ( 19941 ) is 



based on the ALI technique (s ee Sect.[7l), with u se of the OAB (Olson-Auer-Buchler) 
approximate lambda ope rator ( Olson et al.lll98^ ) and the short characteristics method 



I Kunasz and Aueij|l988l ) :o obtain the f ormal solution along individual ray segments 



This method was used bv lPaletou et al.l (|l993l ') to compute the synthetic profiles of H I 



(La), Ca II (H & K), and Mg II (h & k) lines in a 2D horizontally infinite slab model 
using the two-level atom approximation. The authors confirmed a strong influence of 
PRD on synthetic La profiles. 

The models described above treat the 2D radiative transfer problem with the sim- 
plified two-level atom approx imation. An ALI method for mult i-level radiative trans- 
fer was further developed by iRvbicki and Hummed (|l99ll . Il992l ) (the so called MALI 



method - Multilevel Accelerated Lambda Iterations - see Sect.[7|. It is based on cou- 
pled iterative solution of radiative transfer using ALI and preconditioned equations of 
statistical equilibrium for the level populations. Thi s method was suc cessfull y imple- 
mente d into 2D horizontally infinite slab ge ometry by Auer et al. lll994h and bv lPaletoul 



( 19951 ). both for CRD and PRD. Recentlv. [lleger and Paletoul (|2009l ) computed 
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gent neutral helium line profiles (including the He I 10830 A and D3 multiplets) using 
2D horizontal prominence fine structure models in both single-thread and multi-thread 
configurations. They considered isothermal and isobaric horizontal slabs irradiated by 
the solar surface on the sides and on the bottom, and also by the corona on the sides 
and on the top. Radiative transfer computations were performed for hydrogen and 
He I transitions. The authors obtained a good agreement with the results of previous 
models, and confirmed the importance of the multi-thread approach for the modelling 
of the prominence fine structures. 

In the case of the horizontally infinite slab geometry, the incident radiation from 
the solar surface irradiates the bottom and both side surfaces of the prominence slab. 
The top of the slab is not irradiated if one does not consider radiation coming from 
the solar corona (Fig. 135b .). The situation is somewhat different if the prominence 
model assumes a vertical geometry, with the infinite axis pointing perpendicularly 
to the solar surface. Then th e incident radiation irradiates the prominence structure 
from all sides (see Fig. [35b') . ICouttebrozel (|2005l ) developed 2D vertical models with 
cylindrical cross-section. In this case, the plasma parameters and the radiation field 
var y with the di s tance to the cylinder axis and with azimuth. The ALI method used 
by ICouttebrozj (|2005l ) for the solution of the non-LTE radiative transfer in cylin- 
drical geometry with assumption of the two-level atom and CRD was successfully 
modified to t reat al so the multi- l evel a tom with anisotropic incident radiation by 
iGouttebroz i l|2006l ). iGouttebrozj (|2007l) studied the temperature relaxation of 2D 
cylindrical threads to the radiative equilibrium. Gouttcbrozc (200§) developed 2D 
cylindrical models which take into account the Doppler effects produced by 3D ve- 
locity fields. The 2D prominence models with cyl indrical geometry can be used also for 
other solar features such as coronal loops (see e.o.. |Patsourakos et al.|[2"007l : lGouttebrozel 
bOOSi ). The radiative transfer calculations performed under cylindrical geometry could 
be used to provide additional constraints to the models described in Sect. 4.3.3 of 
Paper II to study the effect of partial ionization on the damping of MHD waves. 



9.2 2D Models of Prominence Fine- Structure Threads 

High-resolution prominence observations reveal that solar prominences are not ho- 
mogeneous bodies, but exhibit a variety of fine structures, in particular elongated 
thread-like features and knots of plasma. A considerable number of prominence fine- 
structure threads have nearly vertical orientation, although observations of quiescent 
prominen ces composed m ainly of horizontal threads also exist (see Sect. [J] and also the 
review by iHeinzell |2007^ . Vertical threads which are perpendicular to the mean mag- 
netic field in the prominence could be formed as plasma condensations propagating 
ve rtically due to the int eraction of the neighbouring magnetic field dips, as suggested 
by [Poland and Mariskal jl988) . 



Heinzel and Anzeij I2OOII) described MHS equilibrium of such vertical 2D fine- 



structure threads hanging in magnetic dips of the horizonta l magnetic field. This was 
obtained by a generalization of the ID prominence models of lAnzer and Heinzell (|l999l ) 
in MHS equilibrium of Kippenhahn-Schliiter type. The multi-l evel non-LTE r adiative 
transfer is solved in this 2D structure using MALI technique (lHeinzel|ll995l : IPaletorJ 
with short characteristics method t o obtain the formal solution along individ- 
ual ray segments l Auer and Paletoulll994 ). The 2D temperature variation fFig. I36() . 



prescribed semi-empirically, accounts for two different PCTRs. The PCTR across the 
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Fig. 36 Semi-empirical variation of the temperature within the cross-section of the 2D vertical 
prominence fine-structure thread. The x and y axis, representing geometrical dimensions of 
the thread, are not drawn to the same scale. Iso-contours of the temperature are shown. This 
plot clearly shows two different shapes of the PCTR, with A) a gentle rise of the temperature 
along the x-axis (along the field lines), and B) a steep gradient of the temperature along the 
j/-axis (across the magnetic field lines). 



magnetic field lines is very narrow (Fig. I36B ). with steep gradient of the temperature 
from the cool central part of the thread towards its boundaries. The PCTR along the 
magnetic field hnes is much more extended (Fig. I36I A_). with a shallow rise of the tem- 
perature. Synthetic profiles of the hydrogen Lyman lines obtained along and across 
the magnetic field lines show significant differences with considerably reversed profiles 
obtained across the magnetic field and unreversed profiles obtained along the field. 
Such a b ehaviour is in agreeme nt with findings based on the SUMER observations pre- 
sented bv lHeinzel et al.l ll2001al). A modified version of th e 2D prominence fine- structure 



thread model of iHeinzel and Anzerl 1 200 j ) was used bv iHeinzel et al.l 1 20051 ) to com- 
pute a grid of 18 2D models. The 12-level plus continuum hydrog en model atom and 
adap tive geometrical grid for 2D MHS equilibrium computations ( Heinzel and Anzerl 
I2OO3I ) was used. This study confirmed previous findings about the beh a viour of pro- 
files of the hydrogen Lyman series - see also Schmieder et all (|2007l ). [g unar et al.l 
( 2007al ) supplemented the previous study with the investigation of the Lyman contin- 
uum behaviour and showed that a comparison between synthetic and observed Lyman 
continuum profiles in the wavelength range 800-911 A produces new useful constraints 
on the temperature variation of the prominence fine structures (see the discussion on 
the determination of the electron temperature from the slope of the hydrogen Lyman 
continuum in Sect. I8.ITI]) . 

The same 2D thread model was used bv lGunar et al. 1 l|2007bl ) to derive the prop- 
erties of the prominence fine-structure threads. A trial-and-error method was used to 
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find a model with the best agreement between synthetic profiles of the hydrogen Ly- 
man series and the observed ones. The observed profiles were obtained by SUMER 
on May 25, 2005 and contain the full hydrogen Lyman line series, including unique 
La line observations outside the attenuator. The prominence fine-structure properties 
ob tained in this way are in a good agreement with values given in the textbooks, such 
as iTandberg-HanssenI (199^). However, the authors concluded that the observed Ly- 
man line profiles can be better reproduced by using multi-thread fine-structure models 
consisting of a set of identical 2D threads hanging on separate magnetic field lines with 

the LOS perpendicular to the magnetic field. 

Such multi-thread models were used also bv lCunar et al.l ( 2008l l to study asymme- 
tries of the hydrogen Lyman line profiles. For this reason each thread of the multi-thread 
model has a randomly assigned LOS velocity. A set-up of such multi-thread models 
is shown in Fig. 1371 This model consists of 10 identical 2D fine-structure threads ran- 
domly shifted with respect to the foremost thread (thread #10). LOS velocities are also 
generated randomly, in this case from the interval (—10, 10) km s~^, and their orienta- 
tion and magnitudes are indicated by arrows. Synthetic Lyman line profiles obtained 
by such multi-thread modelling exhibit substantial asymmetries, even with the LOS 
velocities being only of the order of 10 km s~"'^(note that no mutual radiative interac- 
tion is taken into account). The La profiles can also exhibit an opposite asymmetry to 
that of the higher Lyman lines. This is in agreement with the (so f ar) ra t her pu zzling 
behaviour of observed La and L/3 line profiles noticed already bv IViall ( 1982al l. and 
confirmed bv iGunar et al.i ()2007b ') and iVial et al.i (120071 ). 



10 Energy Balance Considerations 

After the work of iHeaslev and MihalasI (|l976h . who demonstrated that the radiative 
equilibrium temperature inside ID prominence slabs is much lower than what is typi- 
cally deduced from spectral observations (see temperature diagnostics in Sects. [5] and 
Is}, various authors have attempted to investigate possible sources of prominence heat- 
ing. 

The radiative equilibrium means that the radiative fiux Fr integrated over all fre- 
quencies is conserved inside the prominence slab: 

dFr 



dx 



Fr 



/ Ft^dp 
Jo 



This also means that the net radiation energy emitted at a given point must be exactly 
equal to that absorbed: 

L^A-K inv - XuJv)du ^ A-K xASu - Jv)du = Q . (55) 
Jo Jo 

In a general case there must exist sources of heating and cooling which will establish 
the energy equilibrium. The energy balance equation then reads: 
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35 




Fig. 37 Multi-thread 2D fine-structure model with randomly assigned LOS velocities, in- 
dicated by the arrows at the top. The grey-scale represents the variation of the tempera- 
ture from the minimum central temperature (black) to the maximum boundary temperature 
(white). White contours represent the density from c entral 5.0 X lO"'^^ g cm~^ to external 
0.05 X 10-13 g cm-3 values. From lGunar et a.1.1 ||2008| '). 
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where F — Fc + Fv represents the conductive flux and the flux associated with the 
mass flows (enthalpy flux), and L are the so-called radiative losses defined in (|55p . 
Other possible sources of the heating or cooling are included in the term H {e.g., wave 
heating, magnetic and mechanical heating, shock dissipation). This term is usually 
considered to be proportional to gas density. This energy balance equation constrains 
the temperature structure of the prominence (ID slab in the above case). 

The evaluation of the radiative losses inside the cool prominence represents a dif- 
ficult task because as we see from (I55|l . L depends on the source function and on the 
radiation field. Moreover, the integration over frequencies includes various lines and 
continua of species which are supposed to be important coolants of the prominence 
plasmas. In the chromosphere, where the thermodynamic conditions are similar to 
prominences, the mo st important coolants are hydrogen, calcium and magnesium (see 
IVernazza et aLlFigSll l. In the case of prominences, no comprehensive stu d y was devoted 
to radiative losses of elements other than hvdrogen. IPoland and Anzeij ( 197lh made a 
crude estimate of Ca II H and K line losses and claimed the i r rela tive importance. 
However, more detailed transfer modelling of IZhang and Faii3 (|l987l l has shown that 
ca lcium losses are negligible for prominence slabs and a similar result was also obtained 
bv lAnzer and Heinzell ( 19991 ). On the other hand, some authors have considered the op- 
tically thin losses in addition to hydro gen ones computed in detail from non-LTE leve l 
populations and radiation intensities jFontenl^_et^ 19961: Anzer and Heinzel Il999l l , 



or added also the helium component ( KuhranT^olandlll99ll )." 

Concerning the heating, the classical conductive-heating term in (|56p can be gener- 
alised by assuming that the P CTR layers are a l so hea ted due to the ambipolar diffusion. 
This was studied in detail bv iFontenla et aD ( 1996ll who have demonstrated that in- 
clusion of the ambipolar diffusion dramatically increases e.g., the L/3 line intensity in 
multi-slab models. However, such large i ntensities are inc onsistent with SUMER ob- 
servations of hydrogen Lyman lines (e.o.. lvial et al.lboOTi '). It is not clear whether the 
problem lies in the method used to compute the Lyman lines with PRD under the 
energy balance with ambipolar diffusion, or whether the ambip olar diffusion itself i s 
improperly treated in case of prominences (see the discussion in lHeinzel et al.ll2005l ). 
This certainly deserves further attention. 

Prominence heating associated with the mass infiow was estimated for various ID 
models bv i Anzer and Heinzel ( 2000. ) who considered the inflow of enthalpy and ioniza- 
tion energy. The resulting energy gains were c o mpar ed with integrated radiative losses 
obtained for such slabs bv lAnzer and Heinzell ( 19991 ). For reasonable inflow velocities, 
many of the considered models can be in energy equilibrium. On the observational side, 
we have a very large range of velocity values and orientations (see Tables [S] and O . 
The magnitude of velocities is in the bracket 2 — 35 km s~^ with some indication of 
increased magnitude for increased temperature (compare Table [5] with Table [G]) . But 
the orientation of these velocities with respect to the prominence axis is unclear. More- 
over, one would suspect that in the frame of a symmetrical prominence where most 
UV lines are formed in the PCTR, the inflow velocities, which would be of opposite 
signs, cancel each other. In fact, there is no direct evidence of infiow of hot material 
which supports the enthalpy contribution. We also have no values of the radiative losses 
corresponding to the observed prominences. Consequently, it is difficult to draw any 
conclusion concerning the role of enthalpy in the energy budget. However, having in 
mind the values of radiative losses in the modelling ofiAnzcr and Heinzel (2000) (which 
represent some observations) on one hand, and noting that the velocity magnitude is 
relatively high in transition region lines (more than 10 km s~^ in Table |6] compared 



81 



to 1.7 km s ^ in the work of lAnzer and Heinzelh on the other ha nd, one can derive 
that t he enthalpy flux can match the radiative losses of most of the Unzer a.nd Heinzell 
( 20001 ) models. The issue of comparing enthalpy (derived from observed mass flows in 
their whole complexity) with radiative losses also derived from observations is certainly 

an important one to be addressed. 

Returning back to radiative equilibrium models, iHeaslev and Mihalad (|l976l ') have 
also demonstrated that the penetration of the hydrogen Lyman continuum radiation 
into a filamentary prominence causes an increase of the central temperature by about 
10^ K. The reason is that the radiation temperature of the incident Lyman contin- 
uum is higher as compared to that of Balmer continuum. The latter one is optically 
thin in all prominence models, and thus fully penetrates into the central parts. This 
result suggests that for geometrically thin fine-structure threads, the central radiative 
equilibrium temperature will be highe r compared to more compa ct and thick slabs 
where it can reach 46 00 K according to Heaslev and Mihalad ( 19761 ) . Such a behaviour 
was demonstrated by iGouttebrozd ( 20071 ) for vertical cvlinders. By increasing the gas 
pressure, the equilibrium temperature decreases beca use the radiative losses are h igher 
for higher density. However, contrary to the resu lts of Heasle v and MihalasI ( 1976h . the 
radiative equilibrium temperatures obtained bv iGouttebrozg i 2007 ) are much higher, 
re aching more than 9000 K for low-pressure models. A similar result was also obtained 
bv lAnzer and He inzelj |l 999|). see their F igs. 7 and 8. By adding losses due to other ele- 
ments (as e.g.. iri lKuhiand Polandlll99lh , the central temperatures may again decrease 
unless some non-radiative heating is considered. This complex behaviour deserves fur- 
ther study based on detailed non-LTE treatment of radiative losses. 

Energy balance in hotter parts of the prominence weis studi ed by v arious authors us- 
ing the DEM (see Sect. 13.2"!^ . Older w ork was summarized by Engvol d (1989 ). Amon, 



Anzer and Heinze" 



recent studies let us mention those of IParenti and Viall (2007) or 
(|2008l ). Within the hotter prominence plasmas, for temperatures above 30000 K, the 
energy balance was treated in dynamical models of prominence condensations - see 
Paper II (Sect. 3). 



11 Open Issues 

This review focuses on results on the thermodynamic parameters brought by spectral 
inversion techniques and non-LTE modelling of the prominence plasma. While most 
of the thermodynamic paramete rs are found to be in overall agreement with the Hvar 
Reference Model Atmosphere of lEngvold et ak I (|l990l ) - see Table [B we still notice a 



large range of values concerning densities, ionization degree, filling factor, etc. Such 
a large range implies uncertainties on the thermodynamic parameters, a variety of 
conditions within an observed prominence and also from prominence to prominence. 
As we have seen with the case of the Lyman series, it may also imply that the properties 
which are derived depend on the viewing angle, the spatial resolution, or other factors. 
It is true that we are now starting to obtain high resolution ground-based and space- 
based observations, and to be able to perform sophisticated non-LTE modelling of 
prominences with fine structures and complex geometries. As a result, today's picture 
of a prominence is that of a dynamic bundle of small-scale structures in constant 
evolution. The density, temperature, and velocity of the plasma all change with time 
and location, and can vary greatly. This is vitally connected to the key questions of 
the magnetic structure and the formation and evolution of prominences, topics that 
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are discussed in Paper II. It is tempting to ascribe specific properties to any observed 
prominence and perform modelling adapted to each observational situation. However, 
we must remember that there are some constraints that all prominences are subjected 
to, e.g., mechanical support, energetic balance, or long lifetime. In order to be able 
to impose these constraints, it is necessary to constrain the diagnostic and to obtain 
typical (if not canonical) values of quiescent and active prominences. 

With these considerations in mind, below we attempt to describe some of the 
outstanding issues that need to be addressed in the future, both observationally and 
theoretically. 

1. The issue of the determination of the thermodynamic parameters is completely 
interrelated with the issue of the fine structure. 

(a) Geometry. What can imaging observations from space tell us about the 3D 
structure of prominences? It is difficult to pick out small features against the 
disk in the He II 304 A channel of EIT and EUVI. Nevertheless, the EUVI 
observations from STEREO A and B, when A and B are in quadrature, should 
certainly help to derive the 3D fine-scale structure of filaments and prominences. 
Meanwhile, non-LTE modelling should progress towards 3D radiative transfer 
including radiative interactions between threads. Impor tant steps have been 
taken recently (e.o.. lGunar et al.ll20oi : lGouttebrozell2008l '). Non-LTE modelling 
should also be coupled with MHD. 

(b) Filling factors. Even if future instruments lack the spatial resolution neces- 
sary for identifying the prominence fine structure (cool core and PCTR - 
prominence-to-corona transition region), efficient and multi- wavelength spec- 
trometers should allow the derivation of local densities and average squared 
densities (from DEM studies), from which sound filling factor values could be 
derived. 

(c) Ionization degree. It is a critical issue because of its impact on the derivation 
of the mass of the prominence. It also is critical if one tries to understand how 
material motion can take place perpendicularly to the magnetic field lines,, 
or what the effects of partial ionization on MHD waves and oscillations are. 
There is now a need for numerical coupling between oscillations modelling (as 
described in Paper II, Sect. 4) and non-LTE ionization balance modelling. 

2. The prominence-to-corona transition region is a critical region for the transfer 
of conductive energy (or its absence) from the corona to the prominence. 

(a) General properties derived from EM and DEM analyses. In spite of many efforts 
aimed at a careful analysis of the data and intensive modelling development, the 
properties of the PCTR remain uncertain. The sampling of the thermal struc- 
ture through the PCTR is still not detailed enough. For example, the different 
results concerning the slope of the temperature gradient for logT < 5 need 
to be understood to engage a correct investigation on the heating mechanisms 
in prominences. This also includes identifying the location of the temperature 
where the DEM has its minimum. How to invert this information to constrain 
modelling? Some attempts have been made, but a more systematic interaction 
between theory and observations is needed. 

(b) The PCTR velocities. The reduced spatial resolution of the UV-EUV compared 
to optical observations maintains some open questions on PCTR small-scale 
structure and magnetic configuration. The NTV (and DEM) may give access 
to unresolved spatial scales properties. It is still not clear if the NTV at a 
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given temperature changes with the location inside the prominence. This can 
indicate the location of more intense activity inside the prominence (including 
heating deposition). Then, much effort needs to be put in the interpretation 
of this component of the plasma velocity, either in terms of waves and micro- 
turbulence, or flows in a small-s cale geometry. I fere a gain, the role of the spatial 
resolution should be reminded. Engvold et al.l (|l980h found an anti-correlation 
between the resolution and the Doppler velocities, and a correlation between 
the spatial resolution and the derived NTV. 
(c) PCTR mass flows compared to cool mass flows. A pattern has become appar- 
ent, especially in filaments observed against the disk, in which roughly vertical 
motions appear in barbs, and horizontal ones appear in the prominence spine. 
How does this observed pattern connect with observations of prominences at 
the limb? What is the connection between observed Ha and EUV motions? 
Are they from separate elements of the prominence, or aspects of the same 
motions? Current observations suggest that such motions tend to be faster in 
the EUV than in Ha. Is this a correct picture, or is it simply due to limitations 
in the EUV observations? How far can individual moving features be traced? 
What ramifications are there in these comparisons for our understanding of the 
PCTR? 

3. The long-standing issue of energy balance still needs careful investigation. 

(a) Spatial and temporal considerations. On which spatial and time scales does 
the prominence heating take place? Energy balance issues including the role of 
radiative equilibrium, the importance of ambipolar diffusion, and the heating 
by waves or shocks, need to be carefully studied using sophisticated non-LTE 
models. 

(b) Energy losses due to radiation emitted by the prominence plasma. What would 
be the modelled time-dependent evolution of prominence condensations (see 
Paper H Sect. 3) if these simulations were considering temperatures lower than 
30000 K at which hydrogen lines form? Efforts should be made to couple the 
non-LTE radiative transfer calculations to other types of prominence plasma 
numerical modelling. This modelling will then lead to radiative signatures and 
yield new observational scenario. 

(c) Enthalpy. As we have seen in Sect. [TOl we do not have precise information on 
the directions of the observed flows with respect to the cool core of prominence 
structures. A determination of the direction and magnitude of velocities, along 
with a precise measurement of densities in a large range of temperatures, is 
necessary in order to evaluate the enthalpy contribution. 

4. Modelling of filaments. 

What do observations of filaments on the disk in EUV optically thick lines, in 
Hq, and in He I 10830 A tell us about their physical state? How can we relate 
the observations to the filaments thermal and magnetic structure? What would 
be a monochromatic appearance of fine-structures corresponding to various mag- 
netic models {e.g., dips)? In the past decade, much effort has been devoted to the 
modelling of the radiation emitted by prominences on the limb. There should be a 
similar focus on filaments. The fine-scale structure and motions of filaments should 
be consistently taken into account in the modelling. 

5. Modelling vs observations. 

How can we explain the multi-wavelength and multi-temperature appearance of 
prominences and filaments and of their surrounding cavities? How would synthetic 
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spectra from non-LTE prominence models compare with a full reference spectrum 
obtained by an EUV spectrometer? To answer this, we need to extend the non-LTE 
modelling to other species. This will give an unprecedented insight on the thermo- 
dynamic structure of prominences, from their cool core to the boundary with the 
corona. Much effort needs to be put in deriving prominence abundances, an impor- 
tant parameter for the choice of the best scenario of prominence formation. Such 
a multi-species non-LTE modelling would add constraints and help interpreting 
observations. 

6. Some basic physical processes may need to be revisited. 

Are prominences in statistical equilibrium? Is ionization equilibrium always a valid 
Eissumption? The recombination time scales seem to be comparable to the fine 
structure lifetime, as shown by the SOT movies. It is necessary to investigate the 
consequences of these new observations for modern non-LTE prominence models, 
and to develop time-dependent numerical codes to properly evaluate ionization 
and recombination times, and to bridge the gap between simplified models and 
high-resolution observations. 

Large multi-instrument observing campaigns have highlighted the need to reconcile 
observations performed by various instruments operating at different wavelciigtlis and 
at a large range of spatial and temporal resolutions. The physical processes responsible 
for the emitted radiation in these various wavelength regimes are different. This has 
to be kept in mind when interpreting such observations. A mission offering two simul- 
taneous vantage points, equipped with instruments of high spectral resolution capable 
of simultaneous observations over a large set of lines, would likely help us to answer 
central questions related to the prominence-to-corona transition region. Ambitious mis- 
sions {e.g., in the frame of Cosmic Vision) should aim at sub-one tenth of arcsecond 
resolution. Imaging spectroscopy is an important aspect of many solar observations 
and indeed of prominences and filaments. Having imaging and spectrometry capabili- 
ties on board future spacecraft with a temperature coverage in the range 10^ - 10^ K 
(including key chromospheric and coronal lines in addition to La, Ho, He II 304 A), 
together with density sensitive lines in the range 10^ - lO" cm , would be a major 
asset, as this provides the crucial plasma parameters, along with the necessary spatial 
information, to address the questions outlined above. 
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